In this paper we investigate a class of impulsive differential equations with Dirichlet boundary conditions. Firstly, we define new inner product of <img src="Edit_890fce38-e82b-4f36-be40-9d05e8119b88.png"...In this paper we investigate a class of impulsive differential equations with Dirichlet boundary conditions. Firstly, we define new inner product of <img src="Edit_890fce38-e82b-4f36-be40-9d05e8119b88.png" width="40" height="17" alt="" /> and prove that the norm which is deduced by the inner product is equivalent to the usual norm. Secondly, we construct the lower and upper solutions of (1.1). Thirdly, we obtain the existence of a positive solution, a negative solution and a sign-changing solution by using critical point theory and variational methods. Finally, an example is presented to illustrate the application of our main result.展开更多
Multiplicative noise removal problems have attracted much attention in recent years.Unlike additive noise,multiplicative noise destroys almost all information of the original image,especially for texture images.Motiva...Multiplicative noise removal problems have attracted much attention in recent years.Unlike additive noise,multiplicative noise destroys almost all information of the original image,especially for texture images.Motivated by the TV-Stokes model,we propose a new two-step variational model to denoise the texture images corrupted by multiplicative noise with a good geometry explanation in this paper.In the first step,we convert the multiplicative denoising problem into an additive one by the logarithm transform and propagate the isophote directions in the tangential field smoothing.Once the isophote directions are constructed,an image is restored to fit the constructed directions in the second step.The existence and uniqueness of the solution to the variational problems are proved.In these two steps,we use the gradient descent method and construct finite difference schemes to solve the problems.Especially,the augmented Lagrangian method and the fast Fourier transform are adopted to accelerate the calculation.Experimental results show that the proposed model can remove the multiplicative noise efficiently and protect the texture well.展开更多
The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-...The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.展开更多
The variational principle of minimum free energy(MFEVP) has been widely used in research of soft matter statics.The MFEVP can be used not only to derive equilibrium equations(including both bulk equations and boundary...The variational principle of minimum free energy(MFEVP) has been widely used in research of soft matter statics.The MFEVP can be used not only to derive equilibrium equations(including both bulk equations and boundary conditions),but also to develop direct variational methods(such as Ritz method) to find approximate solutions to these equilibrium equations. We apply these variational methods to study long-range force transmission in nonlinear elastic biopolymer gels.It is shown that the slow decay of cell-induced displacements measured experimentally for fibroblast spheroids in threedimensional fibrin gels can be well explained by variational approximations based on the three-chain model of biopolymer gels.展开更多
The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this ...The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical,Biological,Radiological and Nuclear(CBRN)applications.The code relies on the Finite Element Method(FEM)for bo...A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical,Biological,Radiological and Nuclear(CBRN)applications.The code relies on the Finite Element Method(FEM)for both the fluid and the dispersed solid phases.Starting from the Navier-Stokes equations and a general description of the FEM strategy,the Streamline Upwind Petrov-Galerkin(SUPG)method is formulated putting some emphasis on the related assembly matrix and stabilization coefficients.Then,the Variational Multiscale Method(VMS)is presented together with a detailed illustration of its algorithm and hierarchy of computational steps.It is demonstrated that the VMS can be considered as a more general version of the SUPG method.The final part of the work is used to assess the reliability of the implemented predictor/multicorrector solution strategy.展开更多
In this paper,we are concerned with solutions to the fractional Schrodinger-Poisson system■ with prescribed mass ∫_(R^(3))|u|^(2)dx=a^(2),where a> 0 is a prescribed number,μ> 0 is a paremeter,s ∈(0,1),2 <...In this paper,we are concerned with solutions to the fractional Schrodinger-Poisson system■ with prescribed mass ∫_(R^(3))|u|^(2)dx=a^(2),where a> 0 is a prescribed number,μ> 0 is a paremeter,s ∈(0,1),2 <q <2_(s)^(*),and 2_(s)^(*)=6/(3-2s) is the fractional critical Sobolev exponent.In the L2-subcritical case,we show the existence of multiple normalized solutions by using the genus theory and the truncation technique;in the L^(2)-supercritical case,we obtain a couple of normalized solutions by developing a fiber map.Under both cases,to recover the loss of compactness of the energy functional caused by the doubly critical growth,we need to adopt the concentration-compactness principle.Our results complement and improve upon some existing studies on the fractional Schrodinger-Poisson system with a nonlocal critical term.展开更多
For any s∈(0,1),let the nonlocal Sobolev space X^(s)(ℝ^(N))be the linear space of Lebesgue measure functions fromℝN toℝsuch that any function u in X^(s)(ℝ^(N))belongs to L2(ℝN)and the function(x,y)\longmapsto\big(u(x...For any s∈(0,1),let the nonlocal Sobolev space X^(s)(ℝ^(N))be the linear space of Lebesgue measure functions fromℝN toℝsuch that any function u in X^(s)(ℝ^(N))belongs to L2(ℝN)and the function(x,y)\longmapsto\big(u(x)-u(y)\big)\sqrt{K(x-y)}is in L2(ℝN,ℝN).First,we show,for a coercive function V(x),the subspace E:=\bigg\{u\in X^s(\mathbb{R}^N):\int_{\mathbb{R}^N}V(x)u^2{\rm d}x<+\infty\bigg\}of X^(s)(ℝ^(N))is embedded compactly into L^(p)(ℝ^(N))for p\in[2,2_s^*),where 2_s^*is the fractional Sobolev critical exponent.In terms of applications,the existence of a least energy sign-changing solution and infinitely many sign-changing solutions of the nonlocal Schrödinger equation-{\cal{L}_K}u+V(x)u=f(x,u),\x\in\\mathbb{R}^N are obtained,where-{\cal{L}_K}is an integro-differential operator and V is coercive at infinity.展开更多
We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved e...We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved exactly,yielding a set of Bessel vortices.Additionally,based on linear solution and using variational approximation,the solutions for the full nonlinear system and their ground state phase diagrams are derived,including the vortex states with quantum numbers m=0,1,as well as mixed states.In this work,mixed states in spin-1 spin-orbit coupling(SOC)BEC are interpreted for the first time as weighted superpositions of three vortex states.Furthermore,the results also indicate that under strong Zeeman splitting,the system cannot form localized states.The variational solutions align well with numerical simulations,showing stable evolution and meeting the criteria for long-term observation in experiments.展开更多
In this paper, we present a new integration algorithm based on the discrete Pfaff-Birkhoff principle for Birkhoffian systems. It is proved that the new algorithm can preserve the general symplectic geometric structure...In this paper, we present a new integration algorithm based on the discrete Pfaff-Birkhoff principle for Birkhoffian systems. It is proved that the new algorithm can preserve the general symplectic geometric structures of Birkhoffian systems. A numerical experiment for a damping oscillator system is conducted. The result shows that the new algorithm can better simulate the energy dissipation than the R-K method, which illustrates that we can numerically solve the dynamical equations by the discrete variational method in a Birkhoffian framework for the systems with a general symplectic structure. Furthermore, it is demonstrated that the results of the numerical experiments are determined not by the constructing methods of Birkhoffian functions but by whether the numerical method can preserve the inherent nature of the dynamical system.展开更多
A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO mode...A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.展开更多
A class of E1 Niйo atmospheric physics oscillation model is considered. The E1 Niйo atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conce...A class of E1 Niйo atmospheric physics oscillation model is considered. The E1 Niйo atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conceptual oscillator model should consider the variations of both the eastern and western Pacific anomaly patterns. An E1 Niйo atmospheric physics model is proposed using a method for the variational iteration theory. Using the variational iteration method, the approximate expansions of the solution of corresponding problem are constructed. That is, firstly, introducing a set of functional and accounting their variationals, the Lagrange multiplicators are counted, and then the variational iteration is defined, finally, the approximate solution is obtained. From approximate expansions of the solution, the zonal sea surface temperature anomaly in the equatorial eastern Pacific and the thermocline depth anomaly of the sea-air oscillation for E1 Niйo atmospheric physics model can be analyzed. E1 Niйo is a very complicated natural phenomenon. Hence basic models need to be reduced for the sea-air oscillator and are solved. The variational iteration is a simple and valid approximate method.展开更多
The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inv...The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.展开更多
One of the advantages of the variational iteration method is the free choice of initial guess. In this paper we use the basic idea of the Jacobian-function method to construct a generalized trial function with some un...One of the advantages of the variational iteration method is the free choice of initial guess. In this paper we use the basic idea of the Jacobian-function method to construct a generalized trial function with some unknown parameters. The Jaulent-Miodek equations are used to illustrate effectiveness and convenience of this method, some new explicit exact travelling wave solutions have been obtained, which include bell-type soliton solution, kink-type soliton solutions, solitary wave solutions, and doubly periodic wave solutions.展开更多
From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given....From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and the initial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.展开更多
In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through min...In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.展开更多
In this work, the effects of quantum confinement on the ground state energy of a correlated electron-hole pair in a spherical and in a disc-like quantum dot have been investigated as a function of quantum dot size. Un...In this work, the effects of quantum confinement on the ground state energy of a correlated electron-hole pair in a spherical and in a disc-like quantum dot have been investigated as a function of quantum dot size. Under parabolic confinement potential and within effective mass approximation Ritz's variational method is applied to Hylleraas-like trial wavefunction. An efficient method for reducing the main effort of the calculation of terms like τekh exp (-λτeh) is introduced. The main contribution of the present work is the introduction of integral transforms which provide the calculation of expectation value of energy and the related matrix elements to be done analytically over single-particle coordinates instead of Hvlleraas coordinates.展开更多
In this paper, the shallow water problem is discussed. By treating the incompressible condition as the constraint, a constrained Hamilton variational principle is presented for the shallow water problem. Based on the ...In this paper, the shallow water problem is discussed. By treating the incompressible condition as the constraint, a constrained Hamilton variational principle is presented for the shallow water problem. Based on the constrained Hamilton variational principle, a shallow water equation based on displacement and pressure (SWE-DP) is developed. A hybrid numerical method combining the finite element method for spa- tial discretization and the Zu-class method for time integration is created for the SWE- DP. The correctness of the proposed SWE-DP is verified by numerical comparisons with two existing shallow water equations (SWEs). The effectiveness of the hybrid numerical method proposed for the SWE-DP is also verified by numerical experiments. Moreover, the numerical experiments demonstrate that the Zu-class method shows excellent perfor- mance with respect to simulating the long time evolution of the shallow water.展开更多
Aimed at the problem that the state estimation in the measurement update of the simultaneous localization and mapping(SLAM)method is incorrect or even not convergent because of the non-Gaussian measurement noise,outli...Aimed at the problem that the state estimation in the measurement update of the simultaneous localization and mapping(SLAM)method is incorrect or even not convergent because of the non-Gaussian measurement noise,outliers,or unknown and time-varying noise statistical characteristics,a robust SLAM method based on the improved variational Bayesian adaptive Kalman filtering(IVBAKF)is proposed.First,the measurement noise covariance is estimated using the variable Bayesian adaptive filtering algorithm.Then,the estimated covariance matrix is robustly processed through the weight function constructed in the form of a reweighted average.Finally,the system updates are iterated multiple times to further gradually correct the state estimation error.Furthermore,to observe features at different depths,a feature measurement model containing depth parameters is constructed.Experimental results show that when the measurement noise does not obey the Gaussian distribution and there are outliers in the measurement information,compared with the variational Bayesian adaptive SLAM method,the positioning accuracy of the proposed method is improved by 17.23%,20.46%,and 17.76%,which has better applicability and robustness to environmental disturbance.展开更多
文摘In this paper we investigate a class of impulsive differential equations with Dirichlet boundary conditions. Firstly, we define new inner product of <img src="Edit_890fce38-e82b-4f36-be40-9d05e8119b88.png" width="40" height="17" alt="" /> and prove that the norm which is deduced by the inner product is equivalent to the usual norm. Secondly, we construct the lower and upper solutions of (1.1). Thirdly, we obtain the existence of a positive solution, a negative solution and a sign-changing solution by using critical point theory and variational methods. Finally, an example is presented to illustrate the application of our main result.
文摘Multiplicative noise removal problems have attracted much attention in recent years.Unlike additive noise,multiplicative noise destroys almost all information of the original image,especially for texture images.Motivated by the TV-Stokes model,we propose a new two-step variational model to denoise the texture images corrupted by multiplicative noise with a good geometry explanation in this paper.In the first step,we convert the multiplicative denoising problem into an additive one by the logarithm transform and propagate the isophote directions in the tangential field smoothing.Once the isophote directions are constructed,an image is restored to fit the constructed directions in the second step.The existence and uniqueness of the solution to the variational problems are proved.In these two steps,we use the gradient descent method and construct finite difference schemes to solve the problems.Especially,the augmented Lagrangian method and the fast Fourier transform are adopted to accelerate the calculation.Experimental results show that the proposed model can remove the multiplicative noise efficiently and protect the texture well.
基金supported by the National Key Research and Development Program of China(No.2020YFB1901900)the National Natural Science Foundation of China(Nos.U20B2011,12175138)the Shanghai Rising-Star Program。
文摘The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.
基金supported by the National Science Foundation for Young Scientists of China (Grant No. 12004082)Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme (2019), 2020 Li Ka Shing Foundation Cross-Disciplinary Research (Grant No. 2020LKSFG08A)+3 种基金Provincial Science Foundation of Guangdong (Grant No. 2019A1515110809)Guangdong Basic and Applied Basic Research Foundation (Grant No. 2020B1515310005)Featured Innovative Projects (Grant No. 2018KTSCX282)Youth Talent Innovative Platforms (Grant No. 2018KQNCX318) in Universities in Guangdong Province。
文摘The variational principle of minimum free energy(MFEVP) has been widely used in research of soft matter statics.The MFEVP can be used not only to derive equilibrium equations(including both bulk equations and boundary conditions),but also to develop direct variational methods(such as Ritz method) to find approximate solutions to these equilibrium equations. We apply these variational methods to study long-range force transmission in nonlinear elastic biopolymer gels.It is shown that the slow decay of cell-induced displacements measured experimentally for fibroblast spheroids in threedimensional fibrin gels can be well explained by variational approximations based on the three-chain model of biopolymer gels.
文摘The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金The authors received the funding of the Royal Higher Institute for Defence(MSP16-06).
文摘A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical,Biological,Radiological and Nuclear(CBRN)applications.The code relies on the Finite Element Method(FEM)for both the fluid and the dispersed solid phases.Starting from the Navier-Stokes equations and a general description of the FEM strategy,the Streamline Upwind Petrov-Galerkin(SUPG)method is formulated putting some emphasis on the related assembly matrix and stabilization coefficients.Then,the Variational Multiscale Method(VMS)is presented together with a detailed illustration of its algorithm and hierarchy of computational steps.It is demonstrated that the VMS can be considered as a more general version of the SUPG method.The final part of the work is used to assess the reliability of the implemented predictor/multicorrector solution strategy.
基金supported by the BIT Research and Innovation Promoting Project(2023YCXY046)the NSFC(11771468,11971027,11971061,12171497 and 12271028)+1 种基金the BNSF(1222017)the Fundamental Research Funds for the Central Universities。
文摘In this paper,we are concerned with solutions to the fractional Schrodinger-Poisson system■ with prescribed mass ∫_(R^(3))|u|^(2)dx=a^(2),where a> 0 is a prescribed number,μ> 0 is a paremeter,s ∈(0,1),2 <q <2_(s)^(*),and 2_(s)^(*)=6/(3-2s) is the fractional critical Sobolev exponent.In the L2-subcritical case,we show the existence of multiple normalized solutions by using the genus theory and the truncation technique;in the L^(2)-supercritical case,we obtain a couple of normalized solutions by developing a fiber map.Under both cases,to recover the loss of compactness of the energy functional caused by the doubly critical growth,we need to adopt the concentration-compactness principle.Our results complement and improve upon some existing studies on the fractional Schrodinger-Poisson system with a nonlocal critical term.
基金supported by the NSFC(12261107)Yunnan Key Laboratory of Modern Analytical Mathematics and Applications(202302AN360007).
文摘For any s∈(0,1),let the nonlocal Sobolev space X^(s)(ℝ^(N))be the linear space of Lebesgue measure functions fromℝN toℝsuch that any function u in X^(s)(ℝ^(N))belongs to L2(ℝN)and the function(x,y)\longmapsto\big(u(x)-u(y)\big)\sqrt{K(x-y)}is in L2(ℝN,ℝN).First,we show,for a coercive function V(x),the subspace E:=\bigg\{u\in X^s(\mathbb{R}^N):\int_{\mathbb{R}^N}V(x)u^2{\rm d}x<+\infty\bigg\}of X^(s)(ℝ^(N))is embedded compactly into L^(p)(ℝ^(N))for p\in[2,2_s^*),where 2_s^*is the fractional Sobolev critical exponent.In terms of applications,the existence of a least energy sign-changing solution and infinitely many sign-changing solutions of the nonlocal Schrödinger equation-{\cal{L}_K}u+V(x)u=f(x,u),\x\in\\mathbb{R}^N are obtained,where-{\cal{L}_K}is an integro-differential operator and V is coercive at infinity.
基金supported by the Guangdong Basic and Applied Basic Research Foundation(Grant No.2023A1515110198)the Natural Science Foundation of Guangdong Province,China(Grant Nos.2024A1515030131 and 2021A1515010214)+2 种基金the National Natural Science Foundation of China(Grant Nos.12274077,11905032,and 12475014)the Research Fund of the Guangdong-Hong Kong-Macao Joint Laboratory for Intelligent Micro-Nano Optoelectronic Technology(Grant No.2020B1212030010)the Israel Science Foundation(Grant No.1695/22).
文摘We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved exactly,yielding a set of Bessel vortices.Additionally,based on linear solution and using variational approximation,the solutions for the full nonlinear system and their ground state phase diagrams are derived,including the vortex states with quantum numbers m=0,1,as well as mixed states.In this work,mixed states in spin-1 spin-orbit coupling(SOC)BEC are interpreted for the first time as weighted superpositions of three vortex states.Furthermore,the results also indicate that under strong Zeeman splitting,the system cannot form localized states.The variational solutions align well with numerical simulations,showing stable evolution and meeting the criteria for long-term observation in experiments.
基金supported by the National Natural Science Foundation of China(Grant Nos.11301350,11172120,and 11202090)the Liaoning University Prereporting Fund Natural Projects(Grant No.2013LDGY02)
文摘In this paper, we present a new integration algorithm based on the discrete Pfaff-Birkhoff principle for Birkhoffian systems. It is proved that the new algorithm can preserve the general symplectic geometric structures of Birkhoffian systems. A numerical experiment for a damping oscillator system is conducted. The result shows that the new algorithm can better simulate the energy dissipation than the R-K method, which illustrates that we can numerically solve the dynamical equations by the discrete variational method in a Birkhoffian framework for the systems with a general symplectic structure. Furthermore, it is demonstrated that the results of the numerical experiments are determined not by the constructing methods of Birkhoffian functions but by whether the numerical method can preserve the inherent nature of the dynamical system.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90111011 and 10471039), the National Key Basic Research Special Foundation of China (Grant Nos 2003CB415101-03 and 2004CB418304), the Key Basic Research Foundation of the Chinese Academy of Sciences (Grant No KZCX3-SW-221) and in part by E-Institutes of Shanghai Municipal Education Commission (Grant No N.E03004).
文摘A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.
文摘A class of E1 Niйo atmospheric physics oscillation model is considered. The E1 Niйo atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conceptual oscillator model should consider the variations of both the eastern and western Pacific anomaly patterns. An E1 Niйo atmospheric physics model is proposed using a method for the variational iteration theory. Using the variational iteration method, the approximate expansions of the solution of corresponding problem are constructed. That is, firstly, introducing a set of functional and accounting their variationals, the Lagrange multiplicators are counted, and then the variational iteration is defined, finally, the approximate solution is obtained. From approximate expansions of the solution, the zonal sea surface temperature anomaly in the equatorial eastern Pacific and the thermocline depth anomaly of the sea-air oscillation for E1 Niйo atmospheric physics model can be analyzed. E1 Niйo is a very complicated natural phenomenon. Hence basic models need to be reduced for the sea-air oscillator and are solved. The variational iteration is a simple and valid approximate method.
基金Project supported by the National Natural Science Foundation of China(Grant No.41175025)
文摘The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.
基金National Natural Science Foundation of China under Grant No.10172056
文摘One of the advantages of the variational iteration method is the free choice of initial guess. In this paper we use the basic idea of the Jacobian-function method to construct a generalized trial function with some unknown parameters. The Jaulent-Miodek equations are used to illustrate effectiveness and convenience of this method, some new explicit exact travelling wave solutions have been obtained, which include bell-type soliton solution, kink-type soliton solutions, solitary wave solutions, and doubly periodic wave solutions.
文摘From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and the initial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFC1510400)the National Natural Science Foundation of China(Grant Nos.41975054 and 41930967)the Special Fund for Forecasters of China Meteorological Administration(Grant No.CMAYBY2018-040)。
文摘In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.
文摘In this work, the effects of quantum confinement on the ground state energy of a correlated electron-hole pair in a spherical and in a disc-like quantum dot have been investigated as a function of quantum dot size. Under parabolic confinement potential and within effective mass approximation Ritz's variational method is applied to Hylleraas-like trial wavefunction. An efficient method for reducing the main effort of the calculation of terms like τekh exp (-λτeh) is introduced. The main contribution of the present work is the introduction of integral transforms which provide the calculation of expectation value of energy and the related matrix elements to be done analytically over single-particle coordinates instead of Hvlleraas coordinates.
基金Project supported by the National Natural Science Foundation of China(No.11472067)
文摘In this paper, the shallow water problem is discussed. By treating the incompressible condition as the constraint, a constrained Hamilton variational principle is presented for the shallow water problem. Based on the constrained Hamilton variational principle, a shallow water equation based on displacement and pressure (SWE-DP) is developed. A hybrid numerical method combining the finite element method for spa- tial discretization and the Zu-class method for time integration is created for the SWE- DP. The correctness of the proposed SWE-DP is verified by numerical comparisons with two existing shallow water equations (SWEs). The effectiveness of the hybrid numerical method proposed for the SWE-DP is also verified by numerical experiments. Moreover, the numerical experiments demonstrate that the Zu-class method shows excellent perfor- mance with respect to simulating the long time evolution of the shallow water.
基金Primary Research and Development Plan of Jiangsu Province(No.BE2022389)Jiangsu Province Agricultural Science and Technology Independent Innovation Fund Project(No.CX(22)3091)the National Natural Science Foundation of China(No.61773113)。
文摘Aimed at the problem that the state estimation in the measurement update of the simultaneous localization and mapping(SLAM)method is incorrect or even not convergent because of the non-Gaussian measurement noise,outliers,or unknown and time-varying noise statistical characteristics,a robust SLAM method based on the improved variational Bayesian adaptive Kalman filtering(IVBAKF)is proposed.First,the measurement noise covariance is estimated using the variable Bayesian adaptive filtering algorithm.Then,the estimated covariance matrix is robustly processed through the weight function constructed in the form of a reweighted average.Finally,the system updates are iterated multiple times to further gradually correct the state estimation error.Furthermore,to observe features at different depths,a feature measurement model containing depth parameters is constructed.Experimental results show that when the measurement noise does not obey the Gaussian distribution and there are outliers in the measurement information,compared with the variational Bayesian adaptive SLAM method,the positioning accuracy of the proposed method is improved by 17.23%,20.46%,and 17.76%,which has better applicability and robustness to environmental disturbance.