A closed-form numerical algorithm (CFNA) is analyzed in detail. CFNA iswidely used in mechanical dynamics for periodic solution of second-order original differentialequations (SODE) with periodic time-variant coeffici...A closed-form numerical algorithm (CFNA) is analyzed in detail. CFNA iswidely used in mechanical dynamics for periodic solution of second-order original differentialequations (SODE) with periodic time-variant coefficients. The principle of the algorithm is todiscretize the motion period into many short time intervals, so the coefficient matrices of theequation set are regarded as constant in a time interval. Defects are found in the originalalgorithm in treating the modal coordinates at the two end-nodes and important modifications to thedefects is made for the algorithm. The modified algorithm is finally used to solve the dynamicproblem of a three-ring planetary gear transmission.展开更多
In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two th...In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.展开更多
Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and sha...Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and shape(k),is crucial in describing the actual wind speed data and evaluating the wind energy potential.Therefore,this study compares the most common conventional numerical(CN)estimation methods and the recent intelligent optimization algorithms(IOA)to show how precise estimation of c and k affects the wind energy resource assessments.In addition,this study conducts technical and economic feasibility studies for five sites in the northern part of Saudi Arabia,namely Aljouf,Rafha,Tabuk,Turaif,and Yanbo.Results exhibit that IOAs have better performance in attaining optimal Weibull parameters and provided an adequate description of the observed wind speed data.Also,with six wind turbine technologies rating between 1 and 3MW,the technical and economic assessment results reveal that the CN methods tend to overestimate the energy output and underestimate the cost of energy($/kWh)compared to the assessments by IOAs.The energy cost analyses show that Turaif is the windiest site,with an electricity cost of$0.016906/kWh.The highest wind energy output is obtained with the wind turbine having a rated power of 2.5 MW at all considered sites with electricity costs not exceeding$0.02739/kWh.Finally,the outcomes of this study exhibit the potential of wind energy in Saudi Arabia,and its environmental goals can be acquired by harvesting wind energy.展开更多
Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-fron...Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-front instability problem in fluid-saturated porous rocks is no exception.Since this kind of instability problem has both the conventional(i.e.trivial)and the unconventional(i.e.nontrivial)solutions,it is necessary to examine the effects of different numerical algorithms,which are used to solve chemical dissolution-front instability problems in fluid-saturated porous rocks.Toward this goal,two different numerical algorithms associated with the commonly-used finite element method are considered in this paper.In the first numerical algorithm,the porosity,pore-fluid pressure and acid/solute concentration are selected as basic variables,while in the second numerical algorithm,the porosity,velocity of pore-fluid flow and acid/solute concentration are selected as basic variables.The particular attention is paid to the effects of these two numerical algorithms on the computational simulation results of unstable chemical dissolution-front propagation in fluid-saturated porous rocks.The related computational simulation results have demonstrated that:1)the first numerical algorithm associated with the porosity-pressure-concentration approach can realistically simulate the evolution processes of unstable chemical dissolution-front propagation in chemical dissolution systems.2)The second numerical algorithm associated with the porosity-velocity-concentration approach fails to simulate the evolution processes of unstable chemical dissolution-front propagation.3)The extra differential operation is the main source to result in the failure of the second numerical algorithm.展开更多
The backtracking search optimization algorithm（BSA） is one of the most recently proposed population-based evolutionary algorithms for global optimization. Due to its memory ability and simple structure, BSA has powe...The backtracking search optimization algorithm（BSA） is one of the most recently proposed population-based evolutionary algorithms for global optimization. Due to its memory ability and simple structure, BSA has powerful capability to find global optimal solutions. However, the algorithm is still insufficient in balancing the exploration and the exploitation. Therefore, an improved adaptive backtracking search optimization algorithm combined with modified Hooke-Jeeves pattern search is proposed for numerical global optimization. It has two main parts： the BSA is used for the exploration phase and the modified pattern search method completes the exploitation phase. In particular, a simple but effective strategy of adapting one of BSA＇s important control parameters is introduced. The proposed algorithm is compared with standard BSA, three state-of-the-art evolutionary algorithms and three superior algorithms in IEEE Congress on Evolutionary Computation 2014（IEEE CEC2014） over six widely-used benchmarks and 22 real-parameter single objective numerical optimization benchmarks in IEEE CEC2014. The results of experiment and statistical analysis demonstrate the effectiveness and efficiency of the proposed algorithm.展开更多
There are many population-based stochastic search algorithms for solving optimization problems. However, the universality and robustness of these algorithms are still unsatisfactory. This paper proposes an enhanced se...There are many population-based stochastic search algorithms for solving optimization problems. However, the universality and robustness of these algorithms are still unsatisfactory. This paper proposes an enhanced self-adaptiveevolutionary algorithm （ESEA） to overcome the demerits above. In the ESEA, four evolutionary operators are designed to enhance the evolutionary structure. Besides, the ESEA employs four effective search strategies under the framework of the self-adaptive learning. Four groups of the experiments are done to find out the most suitable parameter values for the ESEA. In order to verify the performance of the proposed algorithm, 26 state-of-the-art test functions are solved by the ESEA and its competitors. The experimental results demonstrate that the universality and robustness of the ESEA out-perform its competitors.展开更多
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.展开更多
We give a study result to analyze a rather different, semi-analytical numerical algorithms based on splitting-step methods with their applications to mathematical finance. As certain subsistent numerical schemes may f...We give a study result to analyze a rather different, semi-analytical numerical algorithms based on splitting-step methods with their applications to mathematical finance. As certain subsistent numerical schemes may fail due to producing negative values for financial variables which require non-negativity preserving. These algorithms which we are analyzing preserve not only the non-negativity, but also the character of boundaries (natural, reflecting, absorbing, etc.). The derivatives of the CIR process and the Heston model are being extensively studied. Beyond plain vanilla European options, we creatively apply our splitting-step methods to a path-dependent option valuation. We compare our algorithms to a class of numerical schemes based on Euler discretization which are prevalent currently. The comparisons are given with respect to both accuracy and computational time for the European call option under the CIR model whereas with respect to convergence rate for the path-dependent option under the CIR model and the European call option under the Heston model.展开更多
The information of seismic response spectra is key to many problems concerned with aseismic structure and is also helpful for earthquake disaster relief if it is generated in time when earthquake happens. While curren...The information of seismic response spectra is key to many problems concerned with aseismic structure and is also helpful for earthquake disaster relief if it is generated in time when earthquake happens. While current numerical calculation methods suffer from poor precision, especially in frequency band near Nyquist frequency, we present a set of improved parameters for precision improvement. It is shown that precision of displacement and velocity response spectra are both further improved compared to current numerical algorithms. A uniform fitting formula is given for computing these parameters for damping ratio range of 0.01-0.9, quite convenient for practical application.展开更多
A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can ...A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.展开更多
A non-incremental time-space algorithm is proposed for numerical. analysis of forming process with the inclusion of geometrical, material, contact-frictional nonlinearities. Unlike the widely used Newton-Raphso...A non-incremental time-space algorithm is proposed for numerical. analysis of forming process with the inclusion of geometrical, material, contact-frictional nonlinearities. Unlike the widely used Newton-Raphson counterpart, the present scheme features an iterative solution procedure on entire time and space domain. Validity and feasibility of foe present scheme are further justiced by the numerical investigation herewith presented.展开更多
A novel heuristic technique has been developed for solving Ordinary Differential Equation (ODE) numerically under the framework of Genetic Algorithm (GA). The method incorporates a sniffer procedure that helps carry o...A novel heuristic technique has been developed for solving Ordinary Differential Equation (ODE) numerically under the framework of Genetic Algorithm (GA). The method incorporates a sniffer procedure that helps carry out a memetic search within the solution domain in the vicinity of the currently found best chromosome. The technique has been successfully applied to the Korteweg- de Vries (KdV) equation, a well-known nonlinear Partial Differential Equation (PDE). In the present study we consider its solution in the regime of solitary waves, or solitons that is first used to convert the PDE into an ODE. It is then shown that using the sniffer technique assisted GA procedure, numerical solution has successfully been generated quite efficiently for the one-dimensional ODE version of the KdV equation in space variable (x). The technique is quite promising for its applications to systems involving ODE equations where analytical solutions are not directly available.展开更多
If we use Littlewood-Paley decomposition, there is no pseudo-orthogonality for Ho¨rmander symbol operators OpS m 0 , 0 , which is different to the case S m ρ,δ （0 ≤δ 〈 ρ≤ 1）. In this paper, we use a spec...If we use Littlewood-Paley decomposition, there is no pseudo-orthogonality for Ho¨rmander symbol operators OpS m 0 , 0 , which is different to the case S m ρ,δ （0 ≤δ 〈 ρ≤ 1）. In this paper, we use a special numerical algorithm based on wavelets to study the L p continuity of non infinite smooth operators OpS m 0 , 0 ; in fact, we apply first special wavelets to symbol to get special basic operators, then we regroup all the special basic operators at given scale and prove that such scale operator’s continuity decreases very fast, we sum such scale operators and a symbol operator can be approached by very good compact operators. By correlation of basic operators, we get very exact pseudo-orthogonality and also L 2 → L 2 continuity for scale operators. By considering the influence region of scale operator, we get H 1 （= F 0 , 2 1 ） → L 1 continuity and L ∞→ BMO continuity. By interpolation theorem, we get also L p （= F 0 , 2 p ） → L p continuity for 1 〈 p 〈 ∞ . Our results are sharp for F 0 , 2 p → L p continuity when 1 ≤ p ≤ 2, that is to say, we find out the exact order of derivations for which the symbols can ensure the resulting operators to be bounded on these spaces.展开更多
The study on the application of Genetic Algorithms(GA) to numerical simulation has been carried out. The simulation with GA is aimed at to realize the operation optimization of the coal fired MHD generator channel. Th...The study on the application of Genetic Algorithms(GA) to numerical simulation has been carried out. The simulation with GA is aimed at to realize the operation optimization of the coal fired MHD generator channel. The computer program for this purpose has been developed. By simulating numerically the operation optimization of IEE’s 25MWt coal fired experimental MHD generator, the feasibility of the application of GA procedure to the MHD power generation field has been verified.展开更多
A new method to reduce the numerical dispersion of the three-dimensional Alternating Di-rection Implicit Finite-Difference Time-Domain (3-D ADI-FDTD) method is proposed. Firstly,the numerical formulations of the 3-D A...A new method to reduce the numerical dispersion of the three-dimensional Alternating Di-rection Implicit Finite-Difference Time-Domain (3-D ADI-FDTD) method is proposed. Firstly,the numerical formulations of the 3-D ADI-FDTD method are modified with the artificial anisotropy,and the new numerical dispersion relation is derived. Secondly,the relative permittivity tensor of the artificial anisotropy can be obtained by the Adaptive Genetic Algorithm (AGA). In order to demon-strate the accuracy and efficiency of this new method,a monopole antenna is simulated as an exam-ple. And the numerical results and the computational requirements of the proposed method are com-pared with those of the conventional ADI-FDTD method and the measured data. In addition the re-duction of the numerical dispersion is investigated as the objective function of the AGA. It is found that this new method is accurate and efficient by choosing proper objective function.展开更多
A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat tr...A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat transfer and chemical reaction of the resin on the flow and temperature are considered. The numerical algorithm of the resin flow and temperature variation in the process of RTM are studied. Its accuracy and convergence are analyzed. The comparison of temperature variations between experimental results and model predictions is carried out for two RTM cases. Result shows that the model is efficient for evaluating the flow and temperature variation in the filling stage of RTM and there is a good coincidence between theory and experiment.展开更多
Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,...Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.展开更多
In almost all frozen soil models used currently, three variables of temperature, ice content and moisture content are used as prognostic variables and the rate term, accounting for the contribution of the phase change...In almost all frozen soil models used currently, three variables of temperature, ice content and moisture content are used as prognostic variables and the rate term, accounting for the contribution of the phase change between water and ice, is shown explicitly in both the energy and mass balance equations. The models must be solved by a numerical method with an iterative process, and the rate term of the phase change needs to be pre-estimated at the beginning in each iteration step. Since the rate term of the phase change in the energy equation is closely related to the release or absorption of the great amount of fusion heat, a small error in the rate term estimation will introduce greater error in the energy balance, which will amplify the error in the temperature calculation and in turn, cause problems for the numerical solution convergence. In this work, in order to first reduce the trouble, the methodology of the variable transformation is applied to a simplified frozen soil model used currently, which leads to new frozen soil scheme used in this work. In the new scheme, the enthalpy and the total water equivalent are used as predictive variables in the governing equations to replace temperature, volumetric soil moisture and ice content used in many current models. By doing so, the rate terms of the phase change are not shown explicitly in both the mass and energy equations and its pre-estimation is avoided. Secondly, in order to solve this new scheme more functionally, the development of the numerical scheme to the new scheme is described and a numerical algorithm appropriate to the numerical scheme is developed. In order to evaluate the new scheme of the frozen soil model and its relevant algorithm, a series of model evaluations are conducted by comparing numerical results from the new model scheme with three observational data sets. The comparisons show that the results from the model are in good agreement with these data sets in both the change trend of variables and their magnitude values, and the new scheme, together with the algorithm, is more efficient and saves more computer time.展开更多
A three-dimensional,two-phase,five-component mathematical model has been developed to describe flow characteristics of clay particles and flocs in the profile control process,in which the clay particle suspension is i...A three-dimensional,two-phase,five-component mathematical model has been developed to describe flow characteristics of clay particles and flocs in the profile control process,in which the clay particle suspension is injected into the formation to react with residual polymer.This model considers the reaction of clay particles with residual polymer,apparent viscosity of the mixture,retention of clay particles and flocs,as well as the decline in porosity and permeability caused by the retention of clay particles and flocs.A finite difference method is used to discretize the equation for each component in the model.The Runge-Kutta method is used to solve the polymer flow equation,and operator splitting algorithms are used to split the flow equation for clay particles into a hyperbolic equation for convection and a parabolic equation for diffusion,which effectively ensures excellent precision,high speed and good stability.The numerical simulation had been applied successfully in the 4-P1920 unit of the Lamadian Oilfield to forecast the blocking capacity of clay particle suspension and to optimize the injection parameters.展开更多
A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and...A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.展开更多
基金This project is supported by National Natural Science Foundation of China (No.50205019) Development Foundation of Shanghai Municipal Commission of Education, China (No.04EB03).
文摘A closed-form numerical algorithm (CFNA) is analyzed in detail. CFNA iswidely used in mechanical dynamics for periodic solution of second-order original differentialequations (SODE) with periodic time-variant coefficients. The principle of the algorithm is todiscretize the motion period into many short time intervals, so the coefficient matrices of theequation set are regarded as constant in a time interval. Defects are found in the originalalgorithm in treating the modal coordinates at the two end-nodes and important modifications to thedefects is made for the algorithm. The modified algorithm is finally used to solve the dynamicproblem of a three-ring planetary gear transmission.
基金supported by project XJZ2023050044,A2309002 and XJZ2023070052.
文摘In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.
基金The author extends his appreciation to theDeputyship forResearch&Innovation,Ministry of Education,Saudi Arabia for funding this research work through the Project Number(QUIF-4-3-3-33891)。
文摘Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and shape(k),is crucial in describing the actual wind speed data and evaluating the wind energy potential.Therefore,this study compares the most common conventional numerical(CN)estimation methods and the recent intelligent optimization algorithms(IOA)to show how precise estimation of c and k affects the wind energy resource assessments.In addition,this study conducts technical and economic feasibility studies for five sites in the northern part of Saudi Arabia,namely Aljouf,Rafha,Tabuk,Turaif,and Yanbo.Results exhibit that IOAs have better performance in attaining optimal Weibull parameters and provided an adequate description of the observed wind speed data.Also,with six wind turbine technologies rating between 1 and 3MW,the technical and economic assessment results reveal that the CN methods tend to overestimate the energy output and underestimate the cost of energy($/kWh)compared to the assessments by IOAs.The energy cost analyses show that Turaif is the windiest site,with an electricity cost of$0.016906/kWh.The highest wind energy output is obtained with the wind turbine having a rated power of 2.5 MW at all considered sites with electricity costs not exceeding$0.02739/kWh.Finally,the outcomes of this study exhibit the potential of wind energy in Saudi Arabia,and its environmental goals can be acquired by harvesting wind energy.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-front instability problem in fluid-saturated porous rocks is no exception.Since this kind of instability problem has both the conventional(i.e.trivial)and the unconventional(i.e.nontrivial)solutions,it is necessary to examine the effects of different numerical algorithms,which are used to solve chemical dissolution-front instability problems in fluid-saturated porous rocks.Toward this goal,two different numerical algorithms associated with the commonly-used finite element method are considered in this paper.In the first numerical algorithm,the porosity,pore-fluid pressure and acid/solute concentration are selected as basic variables,while in the second numerical algorithm,the porosity,velocity of pore-fluid flow and acid/solute concentration are selected as basic variables.The particular attention is paid to the effects of these two numerical algorithms on the computational simulation results of unstable chemical dissolution-front propagation in fluid-saturated porous rocks.The related computational simulation results have demonstrated that:1)the first numerical algorithm associated with the porosity-pressure-concentration approach can realistically simulate the evolution processes of unstable chemical dissolution-front propagation in chemical dissolution systems.2)The second numerical algorithm associated with the porosity-velocity-concentration approach fails to simulate the evolution processes of unstable chemical dissolution-front propagation.3)The extra differential operation is the main source to result in the failure of the second numerical algorithm.
基金supported by the National Natural Science Foundation of China(61271250)
文摘The backtracking search optimization algorithm（BSA） is one of the most recently proposed population-based evolutionary algorithms for global optimization. Due to its memory ability and simple structure, BSA has powerful capability to find global optimal solutions. However, the algorithm is still insufficient in balancing the exploration and the exploitation. Therefore, an improved adaptive backtracking search optimization algorithm combined with modified Hooke-Jeeves pattern search is proposed for numerical global optimization. It has two main parts： the BSA is used for the exploration phase and the modified pattern search method completes the exploitation phase. In particular, a simple but effective strategy of adapting one of BSA＇s important control parameters is introduced. The proposed algorithm is compared with standard BSA, three state-of-the-art evolutionary algorithms and three superior algorithms in IEEE Congress on Evolutionary Computation 2014（IEEE CEC2014） over six widely-used benchmarks and 22 real-parameter single objective numerical optimization benchmarks in IEEE CEC2014. The results of experiment and statistical analysis demonstrate the effectiveness and efficiency of the proposed algorithm.
基金supported by the Aviation Science Funds of China(2010ZC13012)the Fund of Jiangsu Innovation Program for Graduate Education (CXLX11 0203)
文摘There are many population-based stochastic search algorithms for solving optimization problems. However, the universality and robustness of these algorithms are still unsatisfactory. This paper proposes an enhanced self-adaptiveevolutionary algorithm （ESEA） to overcome the demerits above. In the ESEA, four evolutionary operators are designed to enhance the evolutionary structure. Besides, the ESEA employs four effective search strategies under the framework of the self-adaptive learning. Four groups of the experiments are done to find out the most suitable parameter values for the ESEA. In order to verify the performance of the proposed algorithm, 26 state-of-the-art test functions are solved by the ESEA and its competitors. The experimental results demonstrate that the universality and robustness of the ESEA out-perform its competitors.
文摘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.
文摘We give a study result to analyze a rather different, semi-analytical numerical algorithms based on splitting-step methods with their applications to mathematical finance. As certain subsistent numerical schemes may fail due to producing negative values for financial variables which require non-negativity preserving. These algorithms which we are analyzing preserve not only the non-negativity, but also the character of boundaries (natural, reflecting, absorbing, etc.). The derivatives of the CIR process and the Heston model are being extensively studied. Beyond plain vanilla European options, we creatively apply our splitting-step methods to a path-dependent option valuation. We compare our algorithms to a class of numerical schemes based on Euler discretization which are prevalent currently. The comparisons are given with respect to both accuracy and computational time for the European call option under the CIR model whereas with respect to convergence rate for the path-dependent option under the CIR model and the European call option under the Heston model.
基金supported by Science for Earthquake Resilience (XH12032)
文摘The information of seismic response spectra is key to many problems concerned with aseismic structure and is also helpful for earthquake disaster relief if it is generated in time when earthquake happens. While current numerical calculation methods suffer from poor precision, especially in frequency band near Nyquist frequency, we present a set of improved parameters for precision improvement. It is shown that precision of displacement and velocity response spectra are both further improved compared to current numerical algorithms. A uniform fitting formula is given for computing these parameters for damping ratio range of 0.01-0.9, quite convenient for practical application.
基金This project is supported by National Natural Science Foundation of China (No.59805001)
文摘A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.
文摘A non-incremental time-space algorithm is proposed for numerical. analysis of forming process with the inclusion of geometrical, material, contact-frictional nonlinearities. Unlike the widely used Newton-Raphson counterpart, the present scheme features an iterative solution procedure on entire time and space domain. Validity and feasibility of foe present scheme are further justiced by the numerical investigation herewith presented.
文摘A novel heuristic technique has been developed for solving Ordinary Differential Equation (ODE) numerically under the framework of Genetic Algorithm (GA). The method incorporates a sniffer procedure that helps carry out a memetic search within the solution domain in the vicinity of the currently found best chromosome. The technique has been successfully applied to the Korteweg- de Vries (KdV) equation, a well-known nonlinear Partial Differential Equation (PDE). In the present study we consider its solution in the regime of solitary waves, or solitons that is first used to convert the PDE into an ODE. It is then shown that using the sniffer technique assisted GA procedure, numerical solution has successfully been generated quite efficiently for the one-dimensional ODE version of the KdV equation in space variable (x). The technique is quite promising for its applications to systems involving ODE equations where analytical solutions are not directly available.
基金Supported by the Doctoral programme foundation of National Education Ministry of China
文摘If we use Littlewood-Paley decomposition, there is no pseudo-orthogonality for Ho¨rmander symbol operators OpS m 0 , 0 , which is different to the case S m ρ,δ （0 ≤δ 〈 ρ≤ 1）. In this paper, we use a special numerical algorithm based on wavelets to study the L p continuity of non infinite smooth operators OpS m 0 , 0 ; in fact, we apply first special wavelets to symbol to get special basic operators, then we regroup all the special basic operators at given scale and prove that such scale operator’s continuity decreases very fast, we sum such scale operators and a symbol operator can be approached by very good compact operators. By correlation of basic operators, we get very exact pseudo-orthogonality and also L 2 → L 2 continuity for scale operators. By considering the influence region of scale operator, we get H 1 （= F 0 , 2 1 ） → L 1 continuity and L ∞→ BMO continuity. By interpolation theorem, we get also L p （= F 0 , 2 p ） → L p continuity for 1 〈 p 〈 ∞ . Our results are sharp for F 0 , 2 p → L p continuity when 1 ≤ p ≤ 2, that is to say, we find out the exact order of derivations for which the symbols can ensure the resulting operators to be bounded on these spaces.
文摘The study on the application of Genetic Algorithms(GA) to numerical simulation has been carried out. The simulation with GA is aimed at to realize the operation optimization of the coal fired MHD generator channel. The computer program for this purpose has been developed. By simulating numerically the operation optimization of IEE’s 25MWt coal fired experimental MHD generator, the feasibility of the application of GA procedure to the MHD power generation field has been verified.
基金the National Natural Science Foundation of China (No. 60271012)Research Foundation of ZTE Corporation.
文摘A new method to reduce the numerical dispersion of the three-dimensional Alternating Di-rection Implicit Finite-Difference Time-Domain (3-D ADI-FDTD) method is proposed. Firstly,the numerical formulations of the 3-D ADI-FDTD method are modified with the artificial anisotropy,and the new numerical dispersion relation is derived. Secondly,the relative permittivity tensor of the artificial anisotropy can be obtained by the Adaptive Genetic Algorithm (AGA). In order to demon-strate the accuracy and efficiency of this new method,a monopole antenna is simulated as an exam-ple. And the numerical results and the computational requirements of the proposed method are com-pared with those of the conventional ADI-FDTD method and the measured data. In addition the re-duction of the numerical dispersion is investigated as the objective function of the AGA. It is found that this new method is accurate and efficient by choosing proper objective function.
文摘A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat transfer and chemical reaction of the resin on the flow and temperature are considered. The numerical algorithm of the resin flow and temperature variation in the process of RTM are studied. Its accuracy and convergence are analyzed. The comparison of temperature variations between experimental results and model predictions is carried out for two RTM cases. Result shows that the model is efficient for evaluating the flow and temperature variation in the filling stage of RTM and there is a good coincidence between theory and experiment.
文摘Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.
基金the National Natural Science Foun-dation of China under Grant Nos. 40575043 and 40605024as well as 40730952the National Basic Research Program of China under Grant No. 2009CB421405The Innovation Project of the Chinese Academy of Sci-ences (Grant No. KZCX2-YW-220)
文摘In almost all frozen soil models used currently, three variables of temperature, ice content and moisture content are used as prognostic variables and the rate term, accounting for the contribution of the phase change between water and ice, is shown explicitly in both the energy and mass balance equations. The models must be solved by a numerical method with an iterative process, and the rate term of the phase change needs to be pre-estimated at the beginning in each iteration step. Since the rate term of the phase change in the energy equation is closely related to the release or absorption of the great amount of fusion heat, a small error in the rate term estimation will introduce greater error in the energy balance, which will amplify the error in the temperature calculation and in turn, cause problems for the numerical solution convergence. In this work, in order to first reduce the trouble, the methodology of the variable transformation is applied to a simplified frozen soil model used currently, which leads to new frozen soil scheme used in this work. In the new scheme, the enthalpy and the total water equivalent are used as predictive variables in the governing equations to replace temperature, volumetric soil moisture and ice content used in many current models. By doing so, the rate terms of the phase change are not shown explicitly in both the mass and energy equations and its pre-estimation is avoided. Secondly, in order to solve this new scheme more functionally, the development of the numerical scheme to the new scheme is described and a numerical algorithm appropriate to the numerical scheme is developed. In order to evaluate the new scheme of the frozen soil model and its relevant algorithm, a series of model evaluations are conducted by comparing numerical results from the new model scheme with three observational data sets. The comparisons show that the results from the model are in good agreement with these data sets in both the change trend of variables and their magnitude values, and the new scheme, together with the algorithm, is more efficient and saves more computer time.
基金support from the National High Technology Research and Development Program of China (863 Program) ( 2007AA06200)"Taishan Scholars" Construction Project (No. ts20070704)
文摘A three-dimensional,two-phase,five-component mathematical model has been developed to describe flow characteristics of clay particles and flocs in the profile control process,in which the clay particle suspension is injected into the formation to react with residual polymer.This model considers the reaction of clay particles with residual polymer,apparent viscosity of the mixture,retention of clay particles and flocs,as well as the decline in porosity and permeability caused by the retention of clay particles and flocs.A finite difference method is used to discretize the equation for each component in the model.The Runge-Kutta method is used to solve the polymer flow equation,and operator splitting algorithms are used to split the flow equation for clay particles into a hyperbolic equation for convection and a parabolic equation for diffusion,which effectively ensures excellent precision,high speed and good stability.The numerical simulation had been applied successfully in the 4-P1920 unit of the Lamadian Oilfield to forecast the blocking capacity of clay particle suspension and to optimize the injection parameters.
文摘A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.
