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.展开更多
Variational quantum algorithms are promising methods with the greatest potential to achieve quantum advantage,widely employed in the era of noisy intermediate-scale quantum computing.This study presents an advanced va...Variational quantum algorithms are promising methods with the greatest potential to achieve quantum advantage,widely employed in the era of noisy intermediate-scale quantum computing.This study presents an advanced variational hybrid algorithm(EVQLSE)that leverages both quantum and classical computing paradigms to address the solution of linear equation systems.Initially,an innovative loss function is proposed,drawing inspiration from the similarity measure between two quantum states.This function exhibits a substantial improvement in computational complexity when benchmarked against the variational quantum linear solver.Subsequently,a specialized parameterized quantum circuit structure is presented for small-scale linear systems,which exhibits powerful expressive capabilities.Through rigorous numerical analysis,the expressiveness of this circuit structure is quantitatively assessed using a variational quantum regression algorithm,and it obtained the best score compared to the others.Moreover,the expansion in system size is accompanied by an increase in the number of parameters,placing considerable strain on the training process for the algorithm.To address this challenge,an optimization strategy known as quantum parameter sharing is introduced,which proficiently minimizes parameter volume while adhering to exacting precision standards.Finally,EVQLSE is successfully implemented on a quantum computing platform provided by IBM for the resolution of large-scale problems characterized by a dimensionality of 220.展开更多
Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-ve...Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-vector product,matrix-matrix product,the sum of two matrices,and the calculation of determinant and inverse matrix.We encode the matrix entries into the probability amplitudes of the pure initial states of senders.After applying proper unitary transformation to the complete quantum system,the desired result can be found in certain blocks of the receiver’s density matrix.These quantum protocols can be used as subroutines in other quantum schemes.Furthermore,we present an alternative quantum algorithm for solving linear systems of equations.展开更多
Considering the coupled nonlinear Schr¨odinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite differ...Considering the coupled nonlinear Schr¨odinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference method, Fourier pseudospectral method and wavelet collocation method for spatial discretizations, a series of high accurate conservative algorithms are presented. The proposed algorithms can preserve the corresponding discrete charge and energy conservation laws exactly, which would guarantee their numerical stabilities during long time computations.Furthermore, several analogous multi-symplectic algorithms are constructed as comparison. Numerical experiments for the unstable plane waves will show the advantages of the proposed algorithms over long time and verify the theoretical analysis.展开更多
In this paper,we prove that Euclid's algorithm,Bezout's equation and Divi-sion algorithm are equivalent to each other.Our result shows that Euclid has preliminarily established the theory of divisibility and t...In this paper,we prove that Euclid's algorithm,Bezout's equation and Divi-sion algorithm are equivalent to each other.Our result shows that Euclid has preliminarily established the theory of divisibility and the greatest common divisor.We further provided several suggestions for teaching.展开更多
To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before...To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before.We present a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix A.In detail,we decompose the matrices A and A^(2)into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements.For the one-dimensional Poisson equation with different boundary conditions and the d-dimensional Poisson equation with Dirichlet boundary conditions,the number of decomposition terms is less than that reported in[Phys.Rev.A 2023108,032418].Based on the decomposition of the matrix,we design quantum circuits that efficiently evaluate the cost function.Additionally,numerical simulation verifies the feasibility of the proposed algorithm.Finally,the VQAs for linear systems of equations and matrix-vector multiplications with the K-banded Toeplitz matrix T_(n)^(K)are given,where T_(n)^(K)∈R^(n×n)and K∈O(ploylogn).展开更多
A complete mesh free adaptive algorithm （MFAA）, with solution adaptation and geometric adaptation, is developed to improve the resolution of flow features and to replace traditional global refinement techniques in s...A complete mesh free adaptive algorithm （MFAA）, with solution adaptation and geometric adaptation, is developed to improve the resolution of flow features and to replace traditional global refinement techniques in structured grids. Unnecessary redundant points and elements are avoided by using the mesh free local clouds refinement technology in shock influencing regions and regions near large curvature places on the boundary. Inviscid compressible flows over NACA0012 and RAE2822 airfoils are computed. Finally numerical results validate the accuracy of the above method.展开更多
This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.展开更多
A fast hybrid algorithm based on gridless method coupled with finite volume method （FVM） is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm i...A fast hybrid algorithm based on gridless method coupled with finite volume method （FVM） is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.展开更多
In this paper, we establish the existence of four families of simple wave solu- tion for two dimensional compressible full Euler system in the self-similar plane. For the 2 × 2 quasilinear non-reducible hyperboli...In this paper, we establish the existence of four families of simple wave solu- tion for two dimensional compressible full Euler system in the self-similar plane. For the 2 × 2 quasilinear non-reducible hyperbolic system, there not necessarily exists any simple wave solution. We prove the result that there are simple wave solutions for this 4× 4 non- reducible hyperbolic system, its simple wave flow is covered by four straight characteristics λ0 =λ1,λA2, λ3 and the solutions keep constants along these lines. We also investigate the existence of simple wave solution for the isentropic relativistic hydrodynamic system in the self-similar plane.展开更多
This paper introduces Soccer League Competition (SLC) algorithm as a new optimization technique for solving nonlinear systems of equations. Fundamental ideas of the method are inspired from soccer leagues and based on...This paper introduces Soccer League Competition (SLC) algorithm as a new optimization technique for solving nonlinear systems of equations. Fundamental ideas of the method are inspired from soccer leagues and based on the competitions among teams and players. Like other meta-heuristic methods, the proposed technique starts with an initial population. Population individuals called players are in two types: fixed players and substitutes that all together form some teams. The competition among teams to take the possession of the top ranked positions in the league table and the internal competitions between players in each team for personal improvements results in the convergence of population individuals to the global optimum. Results of applying the proposed algorithm in solving nonlinear systems of equations demonstrate that SLC converges to the answer more accurately and rapidly in comparison with other Meta-heuristic and Newton-type methods.展开更多
In this paper, we prove that if p, q are distinct primes, (p,q)≡(1,7) (mod 12) and Legendres symbol pq=1 , then the equation 1+p a=2 bq c+2 dp eq f has only solutions of the form (a,b,c,d,e,f)=...In this paper, we prove that if p, q are distinct primes, (p,q)≡(1,7) (mod 12) and Legendres symbol pq=1 , then the equation 1+p a=2 bq c+2 dp eq f has only solutions of the form (a,b,c,d,e,f)=(t,0,0,0,t,0), where t is a non negative integer. We also give all solutions of a kind of generalized Ramanujan Nagell equations by using the theories of imaginary quadratic field and Pells equation.展开更多
In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinea...In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.展开更多
In this article, the modified simple equation method has been extended to celebrate the exact solutions of nonlinear partial time-space differential equations of fractional order. Firstly, the fractional complex trans...In this article, the modified simple equation method has been extended to celebrate the exact solutions of nonlinear partial time-space differential equations of fractional order. Firstly, the fractional complex transformation has been implemented to convert nonlinear partial fractional differential equations into nonlinear ordinary differential equations. Afterwards, modified simple equation method has been implemented, to find the exact solutions of these equations, in the sense of modified Riemann-Liouville derivative. For applications, the exact solutions of time-space fractional derivative Burgers’ equation and time-space fractional derivative foam drainage equation have been discussed. Moreover, it can also be concluded that the proposed method is easy, direct and concise as compared to other existing methods.展开更多
In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared wi...In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.展开更多
The accuracy of numerical simulations and many other material design calculations, such as the rolling force, rollingtorque, etc., depends on the description of stress-strain relationship of the deformed materials. On...The accuracy of numerical simulations and many other material design calculations, such as the rolling force, rollingtorque, etc., depends on the description of stress-strain relationship of the deformed materials. One common methodof describing the stress-strain relationship is using constitutive equations, with the unknown parameters fitted byexperimental data obtained via plane strain compression (PSC). Due to the highly nonlinear behaviour of the constitutive equations and the noise included in the PSC data, determination of the model parameters is difficult. Inthis paper, genetic algorithms were exploited to optimise parameters for the constitutive equations based on thePSC data. The original PSC data were processed to generate the stress-strain data, and data pre-processing wascarried out to remove the noise contained in the original PSC data. Several genetic optimisation schemes have beeninvestigated, with different coding schemes and different genetic operators for selection, crossover and mutation.It was found that the real value coded genetic algorithms converged much faster and were more efficient for theparameter optimisation problem.展开更多
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.展开更多
Nonlinear equations systems(NESs)are widely used in real-world problems and they are difficult to solve due to their nonlinearity and multiple roots.Evolutionary algorithms(EAs)are one of the methods for solving NESs,...Nonlinear equations systems(NESs)are widely used in real-world problems and they are difficult to solve due to their nonlinearity and multiple roots.Evolutionary algorithms(EAs)are one of the methods for solving NESs,given their global search capabilities and ability to locate multiple roots of a NES simultaneously within one run.Currently,the majority of research on using EAs to solve NESs focuses on transformation techniques and improving the performance of the used EAs.By contrast,problem domain knowledge of NESs is investigated in this study,where we propose the incorporation of a variable reduction strategy(VRS)into EAs to solve NESs.The VRS makes full use of the systems of expressing a NES and uses some variables(i.e.,core variable)to represent other variables(i.e.,reduced variables)through variable relationships that exist in the equation systems.It enables the reduction of partial variables and equations and shrinks the decision space,thereby reducing the complexity of the problem and improving the search efficiency of the EAs.To test the effectiveness of VRS in dealing with NESs,this paper mainly integrates the VRS into two existing state-of-the-art EA methods(i.e.,MONES and DR-JADE)according to the integration framework of the VRS and EA,respectively.Experimental results show that,with the assistance of the VRS,the EA methods can produce better results than the original methods and other compared methods.Furthermore,extensive experiments regarding the influence of different reduction schemes and EAs substantiate that a better EA for solving a NES with more reduced variables tends to provide better performance.展开更多
In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of ti...In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of time-step and space-mesh. We prove that the, approximate solution has convergence independent of iteration times at each time-level. And the L^2 error estimates are given.展开更多
Jaumann rate, generalized Jaumann rate,Fu rate and Wu rate were incorporated into endochronic equations for finite plastic deformation to analyze simple shear finite deformation. The results show that an oscillatory s...Jaumann rate, generalized Jaumann rate,Fu rate and Wu rate were incorporated into endochronic equations for finite plastic deformation to analyze simple shear finite deformation. The results show that an oscillatory shear stress and normal stress response to a monotonically increasing shear strain occurs when Jaumann rate objective model is adopted for hypoelastic or endochronic materials. The oscillatory response is dependent on objective rate adopted,independent on elastoplastic models. Normal stress is unequal to zero during simple shear finite deformation.展开更多
基金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.
基金supported by the National Natural Science Foundation of China under Grant Nos.62172268 and 62302289the Shanghai Science and Technology Project under Grant Nos.21JC1402800 and 23YF1416200。
文摘Variational quantum algorithms are promising methods with the greatest potential to achieve quantum advantage,widely employed in the era of noisy intermediate-scale quantum computing.This study presents an advanced variational hybrid algorithm(EVQLSE)that leverages both quantum and classical computing paradigms to address the solution of linear equation systems.Initially,an innovative loss function is proposed,drawing inspiration from the similarity measure between two quantum states.This function exhibits a substantial improvement in computational complexity when benchmarked against the variational quantum linear solver.Subsequently,a specialized parameterized quantum circuit structure is presented for small-scale linear systems,which exhibits powerful expressive capabilities.Through rigorous numerical analysis,the expressiveness of this circuit structure is quantitatively assessed using a variational quantum regression algorithm,and it obtained the best score compared to the others.Moreover,the expansion in system size is accompanied by an increase in the number of parameters,placing considerable strain on the training process for the algorithm.To address this challenge,an optimization strategy known as quantum parameter sharing is introduced,which proficiently minimizes parameter volume while adhering to exacting precision standards.Finally,EVQLSE is successfully implemented on a quantum computing platform provided by IBM for the resolution of large-scale problems characterized by a dimensionality of 220.
基金supported by the National Natural Science Foundation of China(Grant No.12031004 and Grant No.12271474,61877054)the Fundamental Research Foundation for the Central Universities(Project No.K20210337)+1 种基金the Zhejiang University Global Partnership Fund,188170+194452119/003partially funded by a state task of Russian Fundamental Investigations(State Registration No.FFSG-2024-0002)。
文摘Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-vector product,matrix-matrix product,the sum of two matrices,and the calculation of determinant and inverse matrix.We encode the matrix entries into the probability amplitudes of the pure initial states of senders.After applying proper unitary transformation to the complete quantum system,the desired result can be found in certain blocks of the receiver’s density matrix.These quantum protocols can be used as subroutines in other quantum schemes.Furthermore,we present an alternative quantum algorithm for solving linear systems of equations.
基金Supported by the National Natural Science Foundation of China under Grant No.91130013Hunan Provincial Innovation Foundation under Grant No.CX2012B010+1 种基金the Innovation Fund of National University of Defense Technology under Grant No.B120205the Open Foundation of State Key Laboratory
文摘Considering the coupled nonlinear Schr¨odinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference method, Fourier pseudospectral method and wavelet collocation method for spatial discretizations, a series of high accurate conservative algorithms are presented. The proposed algorithms can preserve the corresponding discrete charge and energy conservation laws exactly, which would guarantee their numerical stabilities during long time computations.Furthermore, several analogous multi-symplectic algorithms are constructed as comparison. Numerical experiments for the unstable plane waves will show the advantages of the proposed algorithms over long time and verify the theoretical analysis.
基金Supported by the Natural Science Foundation of Chongqing(General Program,NO.CSTB2022NSCQ-MSX0884)Discipline Teaching Special Project of Yangtze Normal University(csxkjx14)。
文摘In this paper,we prove that Euclid's algorithm,Bezout's equation and Divi-sion algorithm are equivalent to each other.Our result shows that Euclid has preliminarily established the theory of divisibility and the greatest common divisor.We further provided several suggestions for teaching.
基金supported by the Shandong Provincial Natural Science Foundation for Quantum Science under Grant No.ZR2021LLZ002the Fundamental Research Funds for the Central Universities under Grant No.22CX03005A。
文摘To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before.We present a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix A.In detail,we decompose the matrices A and A^(2)into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements.For the one-dimensional Poisson equation with different boundary conditions and the d-dimensional Poisson equation with Dirichlet boundary conditions,the number of decomposition terms is less than that reported in[Phys.Rev.A 2023108,032418].Based on the decomposition of the matrix,we design quantum circuits that efficiently evaluate the cost function.Additionally,numerical simulation verifies the feasibility of the proposed algorithm.Finally,the VQAs for linear systems of equations and matrix-vector multiplications with the K-banded Toeplitz matrix T_(n)^(K)are given,where T_(n)^(K)∈R^(n×n)and K∈O(ploylogn).
文摘A complete mesh free adaptive algorithm （MFAA）, with solution adaptation and geometric adaptation, is developed to improve the resolution of flow features and to replace traditional global refinement techniques in structured grids. Unnecessary redundant points and elements are avoided by using the mesh free local clouds refinement technology in shock influencing regions and regions near large curvature places on the boundary. Inviscid compressible flows over NACA0012 and RAE2822 airfoils are computed. Finally numerical results validate the accuracy of the above method.
文摘This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.
基金Aeronautical Science Foundation of China (02A52002), National Natural Science Foundation of China(10372043)
文摘A fast hybrid algorithm based on gridless method coupled with finite volume method （FVM） is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.
文摘In this paper, we establish the existence of four families of simple wave solu- tion for two dimensional compressible full Euler system in the self-similar plane. For the 2 × 2 quasilinear non-reducible hyperbolic system, there not necessarily exists any simple wave solution. We prove the result that there are simple wave solutions for this 4× 4 non- reducible hyperbolic system, its simple wave flow is covered by four straight characteristics λ0 =λ1,λA2, λ3 and the solutions keep constants along these lines. We also investigate the existence of simple wave solution for the isentropic relativistic hydrodynamic system in the self-similar plane.
文摘This paper introduces Soccer League Competition (SLC) algorithm as a new optimization technique for solving nonlinear systems of equations. Fundamental ideas of the method are inspired from soccer leagues and based on the competitions among teams and players. Like other meta-heuristic methods, the proposed technique starts with an initial population. Population individuals called players are in two types: fixed players and substitutes that all together form some teams. The competition among teams to take the possession of the top ranked positions in the league table and the internal competitions between players in each team for personal improvements results in the convergence of population individuals to the global optimum. Results of applying the proposed algorithm in solving nonlinear systems of equations demonstrate that SLC converges to the answer more accurately and rapidly in comparison with other Meta-heuristic and Newton-type methods.
文摘In this paper, we prove that if p, q are distinct primes, (p,q)≡(1,7) (mod 12) and Legendres symbol pq=1 , then the equation 1+p a=2 bq c+2 dp eq f has only solutions of the form (a,b,c,d,e,f)=(t,0,0,0,t,0), where t is a non negative integer. We also give all solutions of a kind of generalized Ramanujan Nagell equations by using the theories of imaginary quadratic field and Pells equation.
基金supported by National Foundation of Natural Science under the Grant 11071216
文摘In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.
文摘In this article, the modified simple equation method has been extended to celebrate the exact solutions of nonlinear partial time-space differential equations of fractional order. Firstly, the fractional complex transformation has been implemented to convert nonlinear partial fractional differential equations into nonlinear ordinary differential equations. Afterwards, modified simple equation method has been implemented, to find the exact solutions of these equations, in the sense of modified Riemann-Liouville derivative. For applications, the exact solutions of time-space fractional derivative Burgers’ equation and time-space fractional derivative foam drainage equation have been discussed. Moreover, it can also be concluded that the proposed method is easy, direct and concise as compared to other existing methods.
基金Project supported by the National Natural Science Foundation of China （Grant No.10271074）, and the Special Funds for Major Specialities of Shanghai Education Commission （Grant No.J50101）
文摘In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.
文摘The accuracy of numerical simulations and many other material design calculations, such as the rolling force, rollingtorque, etc., depends on the description of stress-strain relationship of the deformed materials. One common methodof describing the stress-strain relationship is using constitutive equations, with the unknown parameters fitted byexperimental data obtained via plane strain compression (PSC). Due to the highly nonlinear behaviour of the constitutive equations and the noise included in the PSC data, determination of the model parameters is difficult. Inthis paper, genetic algorithms were exploited to optimise parameters for the constitutive equations based on thePSC data. The original PSC data were processed to generate the stress-strain data, and data pre-processing wascarried out to remove the noise contained in the original PSC data. Several genetic optimisation schemes have beeninvestigated, with different coding schemes and different genetic operators for selection, crossover and mutation.It was found that the real value coded genetic algorithms converged much faster and were more efficient for theparameter optimisation problem.
文摘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.
基金This work was supported by the National Natural Science Foundation of China(62073341)in part by the Natural Science Fund for Distinguished Young Scholars of Hunan Province(2019JJ20026).
文摘Nonlinear equations systems(NESs)are widely used in real-world problems and they are difficult to solve due to their nonlinearity and multiple roots.Evolutionary algorithms(EAs)are one of the methods for solving NESs,given their global search capabilities and ability to locate multiple roots of a NES simultaneously within one run.Currently,the majority of research on using EAs to solve NESs focuses on transformation techniques and improving the performance of the used EAs.By contrast,problem domain knowledge of NESs is investigated in this study,where we propose the incorporation of a variable reduction strategy(VRS)into EAs to solve NESs.The VRS makes full use of the systems of expressing a NES and uses some variables(i.e.,core variable)to represent other variables(i.e.,reduced variables)through variable relationships that exist in the equation systems.It enables the reduction of partial variables and equations and shrinks the decision space,thereby reducing the complexity of the problem and improving the search efficiency of the EAs.To test the effectiveness of VRS in dealing with NESs,this paper mainly integrates the VRS into two existing state-of-the-art EA methods(i.e.,MONES and DR-JADE)according to the integration framework of the VRS and EA,respectively.Experimental results show that,with the assistance of the VRS,the EA methods can produce better results than the original methods and other compared methods.Furthermore,extensive experiments regarding the influence of different reduction schemes and EAs substantiate that a better EA for solving a NES with more reduced variables tends to provide better performance.
基金This work was supported by Natural Science Foundation of China and Shandong Province.
文摘In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of time-step and space-mesh. We prove that the, approximate solution has convergence independent of iteration times at each time-level. And the L^2 error estimates are given.
文摘Jaumann rate, generalized Jaumann rate,Fu rate and Wu rate were incorporated into endochronic equations for finite plastic deformation to analyze simple shear finite deformation. The results show that an oscillatory shear stress and normal stress response to a monotonically increasing shear strain occurs when Jaumann rate objective model is adopted for hypoelastic or endochronic materials. The oscillatory response is dependent on objective rate adopted,independent on elastoplastic models. Normal stress is unequal to zero during simple shear finite deformation.
