This paper is devoted to the investigation of the Landau–Ginzburg–Higgs equation(LGHe),which serves as a mathematical model to understand phenomena such as superconductivity and cyclotron waves.The LGHe finds applic...This paper is devoted to the investigation of the Landau–Ginzburg–Higgs equation(LGHe),which serves as a mathematical model to understand phenomena such as superconductivity and cyclotron waves.The LGHe finds applications in various scientific fields,including fluid dynamics,plasma physics,biological systems,and electricity-electronics.The study adopts Lie symmetry analysis as the primary framework for exploration.This analysis involves the identification of Lie point symmetries that are admitted by the differential equation.By leveraging these Lie point symmetries,symmetry reductions are performed,leading to the discovery of group invariant solutions.To obtain explicit solutions,several mathematical methods are applied,including Kudryashov's method,the extended Jacobi elliptic function expansion method,the power series method,and the simplest equation method.These methods yield solutions characterized by exponential,hyperbolic,and elliptic functions.The obtained solutions are visually represented through 3D,2D,and density plots,which effectively illustrate the nature of the solutions.These plots depict various patterns,such as kink-shaped,singular kink-shaped,bell-shaped,and periodic solutions.Finally,the paper employs the multiplier method and the conservation theorem introduced by Ibragimov to derive conserved vectors.These conserved vectors play a crucial role in the study of physical quantities,such as the conservation of energy and momentum,and contribute to the understanding of the underlying physics of the system.展开更多
The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction di...The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.展开更多
The generalized product bi-conjugate gradient(GPBiCG(m,l))method has been recently proposed as a hybrid variant of the GPBi CG and the Bi CGSTAB methods to solve the linear system Ax=b with non-symmetric coefficient m...The generalized product bi-conjugate gradient(GPBiCG(m,l))method has been recently proposed as a hybrid variant of the GPBi CG and the Bi CGSTAB methods to solve the linear system Ax=b with non-symmetric coefficient matrix,and its attractive convergence behavior has been authenticated in many numerical experiments.By means of the Kronecker product and the vectorization operator,this paper aims to develop the GPBi CG(m,l)method to solve the general matrix equation■ and the general discrete-time periodic matrix equations■ which include the well-known Lyapunov,Stein,and Sylvester matrix equations that arise in a wide variety of applications in engineering,communications and scientific computations.The accuracy and efficiency of the extended GPBi CG(m,l)method assessed against some existing iterative methods are illustrated by several numerical experiments.展开更多
For electromagnetic governing equations formulated by magnetic vector potential and electric scalar potential,its detailed numerical implementation is achieved by using meshless method and Galerkin approach.And essent...For electromagnetic governing equations formulated by magnetic vector potential and electric scalar potential,its detailed numerical implementation is achieved by using meshless method and Galerkin approach.And essential boundary and interface condition of electromagnetic field are imposed by means of Lagrange multiplier method.Furthermore,the influences of interpolation point number at essential boundary and interface on computational results are also discussed.Examples are given to validate the effects of meshless method and Lagrange multiplier approach for electromagnetic field.展开更多
We solve the Duffin-Kemmer-Petiau （DKP） equation with a non-minimal vector Yukawa potential in （1＋1）- dimensional spa^e-time for spin-1 particles. The Nikiforov Uvarov method is used in the calculations, and the ...We solve the Duffin-Kemmer-Petiau （DKP） equation with a non-minimal vector Yukawa potential in （1＋1）- dimensional spa^e-time for spin-1 particles. The Nikiforov Uvarov method is used in the calculations, and the eigen- functions as well as the energy eigenvalues are obtained in a proper Pekeris-type approximation.展开更多
We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is app...We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is applied for time discretization.Compared with the average vector field method on the uniform mesh,the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation.This is due to the moving mesh method,which can concentrate the grid points more densely where the solution changes drastically.Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.展开更多
The general approach for solving the nonlinear equations is linearizing the equations and forming various iterative procedures, then executing the numerical simulation. For the strongly nonlinear problems, the solutio...The general approach for solving the nonlinear equations is linearizing the equations and forming various iterative procedures, then executing the numerical simulation. For the strongly nonlinear problems, the solution obtained in the iterative process is always difficult, even divergent due to the numerical instability. It can not fulfill the engineering requirements. Newton's method and its variants can not settle this problem. As a result, the application of numerical simulation for the strongly nonlinear problems is limited. An auto-adjustable damping method has been presented in this paper. This is a further improvement of Newton's method with damping factor. A set of vector of damping factor is introduced. This set of vector can be adjusted continuously during the iterative process in accordance with the judgement and adjustment. An effective convergence coefficient and quichening coefficient are employed to relax the restricted requirements for the initial values and to shorten the iterative process. Then, the numerical stability will be ensured for the solution of complicated strongly nonlinear equations. Using this method, some complicated strongly nonlinear heat transfer problems in airplanes and aeroengines have been numerically simulated successfully. It can be used for the numerical simulation of strongly nonlinear problems in engineering such as nonlinear hydrodynamics and aerodynamics, heat transfer and structural dynamic response etc.展开更多
By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at t...By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.展开更多
In this paper, iterative or successive approximation methods for the Hamilton-Jacobi-Bellman-lsaacs equations （HJBIEs） arising in both deterministic and stochastic optimal control for affine nonlinear systems are de...In this paper, iterative or successive approximation methods for the Hamilton-Jacobi-Bellman-lsaacs equations （HJBIEs） arising in both deterministic and stochastic optimal control for affine nonlinear systems are developed. Convergence of the methods are established under fairly mild assumptions, and examples are solved to demonstrate the effectiveness of the methods. However, the results presented in the paper are preliminary, and do not yet imply in anyway that the solutions computed will be stabilizing. More improvements and experimentation will be required before a satisfactory algorithm is developed.展开更多
In this paper,some existence results for the fourth order nonlinear subelliptic equations on the Heisenberg group are given by means of variational methods.
New two-component vector breather solution of the modified Benjamin-Bona-Mahony(MBBM)equation is considered.Using the generalized perturbation reduction method,the MBBM equation is reduced to the coupled nonlinear Sch...New two-component vector breather solution of the modified Benjamin-Bona-Mahony(MBBM)equation is considered.Using the generalized perturbation reduction method,the MBBM equation is reduced to the coupled nonlinear Schr¨odinger equations for auxiliary functions.Explicit analytical expressions for the profile and parameters of the vector breather oscillating with the sum and difference of the frequencies and wavenumbers are presented.The two-component vector breather and single-component scalar breather of the MBBM equation is compared.展开更多
This paper is concerned with the following second-order vector boundary value problem ：x^R=f（t,Sx,x,x＇）,0〈t〈1,x（0）=A,g（x（1）,x＇（1））=B,where x,f,g,A and B are n-vectors. Under appropriate assumptions,exis...This paper is concerned with the following second-order vector boundary value problem ：x^R=f（t,Sx,x,x＇）,0〈t〈1,x（0）=A,g（x（1）,x＇（1））=B,where x,f,g,A and B are n-vectors. Under appropriate assumptions,existence and uniqueness of solutions are obtained by using upper and lower solutions method.展开更多
The multi-cavity soft actuator is assembled from single-cavity soft actuator through a reasonable geometric distribution.It has the characteristic that the pneumatic soft actuator is driven by its own deformation and ...The multi-cavity soft actuator is assembled from single-cavity soft actuator through a reasonable geometric distribution.It has the characteristic that the pneumatic soft actuator is driven by its own deformation and has more degrees of freedom.Pneumatic soft actuator is widely used as an emerging discipline and its strong compliance has been greatly developed and applied.However,as the most application potential type of soft actuators,there is still a lack of simple and effective deformation prediction methods for studying the spatial deformation of multi-cavity soft actuators.To solve this problem,a vector equation method is proposed based on the analysis of the principle of the space deformation of the two-cavity,three-cavity and four-cavity soft actuators.Furthermore,a nonlinear mathematical model of the air pressure,space position and deformation trajectory of the soft actuator end is established by combining the vector equation method.Finally,the three-channel soft actuator is verified through experiments.The results show that the mathematical model can better predict the space deformation trajectory of the soft actuator,which provides a new research method for studying the space deformation of the multi-channel soft actuator.展开更多
The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain ...The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain can be expressed through its values on the boundary. Boundary integral equations of the second kind for solving boundary-valued problems of the first and second kinds are developed. The result has been also generalised to the case of solenoidal vector fields with potential vorticity. It is shown that the resulting integral equations are Fredholm integral equations of the second kind and allow effective numerical solving of corresponding boundary-valued problems. Examples of numerical solutions for a sphere and an ellipsoid are given for demonstration of efficiency of the offered method.展开更多
In the presence of spin and pseudospin （p-spin） symmetries, the approximate analytical bound states of the Dirac equation for scalar-vector-tensor Hulth6n potentials are obtained with any arbitrary spin-orbit coupli...In the presence of spin and pseudospin （p-spin） symmetries, the approximate analytical bound states of the Dirac equation for scalar-vector-tensor Hulth6n potentials are obtained with any arbitrary spin-orbit coupling number using the Pekeris approximation. The Hulth6n tensor interaction is studied instead of the commonly used Coulomb or linear terms. The generalized parametric Nikiforov-Uvarov （NU） method is used to obtain energy eigenvalues and corresponding wave functions in their closed forms. It is shown that tensor interaction removes degeneracy between spin and p-spin doublets. Some numerical results are also given.展开更多
Alinhac solved a long-standing open problem in 2001 and established that quasilinear wave equations in two space dimensions with quadratic null nonlinearities admit global-in-time solutions,provided that the initial d...Alinhac solved a long-standing open problem in 2001 and established that quasilinear wave equations in two space dimensions with quadratic null nonlinearities admit global-in-time solutions,provided that the initial data are compactly supported and sufficiently small in Sobolev norm.In this work,Alinhac obtained an upper bound with polynomial growth in time for the top-order energy of the solutions.A natural question then arises whether the time-growth is a true phenomenon,despite the possible conservation of basic energy.In the present paper,we establish that the top-order energy of the solutions in Alinhac theorem remains globally bounded in time.展开更多
基金the South African National Space Agency (SANSA) for funding this work
文摘This paper is devoted to the investigation of the Landau–Ginzburg–Higgs equation(LGHe),which serves as a mathematical model to understand phenomena such as superconductivity and cyclotron waves.The LGHe finds applications in various scientific fields,including fluid dynamics,plasma physics,biological systems,and electricity-electronics.The study adopts Lie symmetry analysis as the primary framework for exploration.This analysis involves the identification of Lie point symmetries that are admitted by the differential equation.By leveraging these Lie point symmetries,symmetry reductions are performed,leading to the discovery of group invariant solutions.To obtain explicit solutions,several mathematical methods are applied,including Kudryashov's method,the extended Jacobi elliptic function expansion method,the power series method,and the simplest equation method.These methods yield solutions characterized by exponential,hyperbolic,and elliptic functions.The obtained solutions are visually represented through 3D,2D,and density plots,which effectively illustrate the nature of the solutions.These plots depict various patterns,such as kink-shaped,singular kink-shaped,bell-shaped,and periodic solutions.Finally,the paper employs the multiplier method and the conservation theorem introduced by Ibragimov to derive conserved vectors.These conserved vectors play a crucial role in the study of physical quantities,such as the conservation of energy and momentum,and contribute to the understanding of the underlying physics of the system.
基金supported by the National Natural Science Foundation of China(Grant No.91130013)the Open Foundation of State Key Laboratory of HighPerformance Computing of China
文摘The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.
基金Supported by the National Natural Sciences Foundation of China(Grant Nos.11501079 11571061)Part by the Higher Education Commission of Egypt
文摘The generalized product bi-conjugate gradient(GPBiCG(m,l))method has been recently proposed as a hybrid variant of the GPBi CG and the Bi CGSTAB methods to solve the linear system Ax=b with non-symmetric coefficient matrix,and its attractive convergence behavior has been authenticated in many numerical experiments.By means of the Kronecker product and the vectorization operator,this paper aims to develop the GPBi CG(m,l)method to solve the general matrix equation■ and the general discrete-time periodic matrix equations■ which include the well-known Lyapunov,Stein,and Sylvester matrix equations that arise in a wide variety of applications in engineering,communications and scientific computations.The accuracy and efficiency of the extended GPBi CG(m,l)method assessed against some existing iterative methods are illustrated by several numerical experiments.
基金the National Natural Science Foundation of China(No.50875169)
文摘For electromagnetic governing equations formulated by magnetic vector potential and electric scalar potential,its detailed numerical implementation is achieved by using meshless method and Galerkin approach.And essential boundary and interface condition of electromagnetic field are imposed by means of Lagrange multiplier method.Furthermore,the influences of interpolation point number at essential boundary and interface on computational results are also discussed.Examples are given to validate the effects of meshless method and Lagrange multiplier approach for electromagnetic field.
文摘We solve the Duffin-Kemmer-Petiau （DKP） equation with a non-minimal vector Yukawa potential in （1＋1）- dimensional spa^e-time for spin-1 particles. The Nikiforov Uvarov method is used in the calculations, and the eigen- functions as well as the energy eigenvalues are obtained in a proper Pekeris-type approximation.
基金the National Key R&D Program of China(Grant No.2020YFA0709800)the National Natural Science Foundation of China(Grant Nos.11901577,11971481,12071481,and 12001539)+3 种基金the Natural Science Foundation of Hunan,China(Grant Nos.S2017JJQNJJ0764 and 2020JJ5652)the fund from Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering(Grant No.2018MMAEZD004)the Basic Research Foundation of National Numerical Wind Tunnel Project,China(Grant No.NNW2018-ZT4A08)the Research Fund of National University of Defense Technology(Grant No.ZK19-37)。
文摘We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is applied for time discretization.Compared with the average vector field method on the uniform mesh,the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation.This is due to the moving mesh method,which can concentrate the grid points more densely where the solution changes drastically.Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.
文摘The general approach for solving the nonlinear equations is linearizing the equations and forming various iterative procedures, then executing the numerical simulation. For the strongly nonlinear problems, the solution obtained in the iterative process is always difficult, even divergent due to the numerical instability. It can not fulfill the engineering requirements. Newton's method and its variants can not settle this problem. As a result, the application of numerical simulation for the strongly nonlinear problems is limited. An auto-adjustable damping method has been presented in this paper. This is a further improvement of Newton's method with damping factor. A set of vector of damping factor is introduced. This set of vector can be adjusted continuously during the iterative process in accordance with the judgement and adjustment. An effective convergence coefficient and quichening coefficient are employed to relax the restricted requirements for the initial values and to shorten the iterative process. Then, the numerical stability will be ensured for the solution of complicated strongly nonlinear equations. Using this method, some complicated strongly nonlinear heat transfer problems in airplanes and aeroengines have been numerically simulated successfully. It can be used for the numerical simulation of strongly nonlinear problems in engineering such as nonlinear hydrodynamics and aerodynamics, heat transfer and structural dynamic response etc.
基金Project supported by the National Natural Science Foundation of China (No. 10172038).
文摘By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.
文摘In this paper, iterative or successive approximation methods for the Hamilton-Jacobi-Bellman-lsaacs equations （HJBIEs） arising in both deterministic and stochastic optimal control for affine nonlinear systems are developed. Convergence of the methods are established under fairly mild assumptions, and examples are solved to demonstrate the effectiveness of the methods. However, the results presented in the paper are preliminary, and do not yet imply in anyway that the solutions computed will be stabilizing. More improvements and experimentation will be required before a satisfactory algorithm is developed.
文摘In this paper,some existence results for the fourth order nonlinear subelliptic equations on the Heisenberg group are given by means of variational methods.
文摘New two-component vector breather solution of the modified Benjamin-Bona-Mahony(MBBM)equation is considered.Using the generalized perturbation reduction method,the MBBM equation is reduced to the coupled nonlinear Schr¨odinger equations for auxiliary functions.Explicit analytical expressions for the profile and parameters of the vector breather oscillating with the sum and difference of the frequencies and wavenumbers are presented.The two-component vector breather and single-component scalar breather of the MBBM equation is compared.
文摘This paper is concerned with the following second-order vector boundary value problem ：x^R=f（t,Sx,x,x＇）,0〈t〈1,x（0）=A,g（x（1）,x＇（1））=B,where x,f,g,A and B are n-vectors. Under appropriate assumptions,existence and uniqueness of solutions are obtained by using upper and lower solutions method.
基金the National Natural Science Foundation of China(No.11604205)。
文摘The multi-cavity soft actuator is assembled from single-cavity soft actuator through a reasonable geometric distribution.It has the characteristic that the pneumatic soft actuator is driven by its own deformation and has more degrees of freedom.Pneumatic soft actuator is widely used as an emerging discipline and its strong compliance has been greatly developed and applied.However,as the most application potential type of soft actuators,there is still a lack of simple and effective deformation prediction methods for studying the spatial deformation of multi-cavity soft actuators.To solve this problem,a vector equation method is proposed based on the analysis of the principle of the space deformation of the two-cavity,three-cavity and four-cavity soft actuators.Furthermore,a nonlinear mathematical model of the air pressure,space position and deformation trajectory of the soft actuator end is established by combining the vector equation method.Finally,the three-channel soft actuator is verified through experiments.The results show that the mathematical model can better predict the space deformation trajectory of the soft actuator,which provides a new research method for studying the space deformation of the multi-channel soft actuator.
文摘The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain can be expressed through its values on the boundary. Boundary integral equations of the second kind for solving boundary-valued problems of the first and second kinds are developed. The result has been also generalised to the case of solenoidal vector fields with potential vorticity. It is shown that the resulting integral equations are Fredholm integral equations of the second kind and allow effective numerical solving of corresponding boundary-valued problems. Examples of numerical solutions for a sphere and an ellipsoid are given for demonstration of efficiency of the offered method.
文摘In the presence of spin and pseudospin （p-spin） symmetries, the approximate analytical bound states of the Dirac equation for scalar-vector-tensor Hulth6n potentials are obtained with any arbitrary spin-orbit coupling number using the Pekeris approximation. The Hulth6n tensor interaction is studied instead of the commonly used Coulomb or linear terms. The generalized parametric Nikiforov-Uvarov （NU） method is used to obtain energy eigenvalues and corresponding wave functions in their closed forms. It is shown that tensor interaction removes degeneracy between spin and p-spin doublets. Some numerical results are also given.
基金supported by the China Postdoctoral Science Foundation(2021M690702)The author Z.L.was in part supported by NSFC(11725102)+2 种基金Sino-German Center(M-0548)the National Key R&D Program of China(2018AAA0100303)National Support Program for Young Top-Notch TalentsShanghai Science and Technology Program[21JC1400600 and No.19JC1420101].
文摘Alinhac solved a long-standing open problem in 2001 and established that quasilinear wave equations in two space dimensions with quadratic null nonlinearities admit global-in-time solutions,provided that the initial data are compactly supported and sufficiently small in Sobolev norm.In this work,Alinhac obtained an upper bound with polynomial growth in time for the top-order energy of the solutions.A natural question then arises whether the time-growth is a true phenomenon,despite the possible conservation of basic energy.In the present paper,we establish that the top-order energy of the solutions in Alinhac theorem remains globally bounded in time.
