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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
The variational method is applied to calculate the dispersion characteristics of disc-loaded waveguide slow-wave structures. The parameters describing the waveguide discontinuities in disc-loaded waveguide are calcula...The variational method is applied to calculate the dispersion characteristics of disc-loaded waveguide slow-wave structures. The parameters describing the waveguide discontinuities in disc-loaded waveguide are calculated by the variational method. Then the dispersion characteristics of slow-wave structures are obtained using lossless microwave quadrupole theory. Good agreement was observed between results of the Variational method and those of field matching method and high frequency structure simulator. In the case of broad band, results of the variational method are better than those of field matching method.展开更多
In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary condi...In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.展开更多
In this work, we propose an efficient method of reducing the computational effort of variational calculation with a Hylleraas-like trial wavefunction. The method consists of introducing integral transforms for the ter...In this work, we propose an efficient method of reducing the computational effort of variational calculation with a Hylleraas-like trial wavefunction. The method consists of introducing integral transforms for the terms as r12^κ exp (-λτ12) which provide the calculation of the expectation value of energy and the relevant matrix elements to be done analytically over single-electron coordinates instead of Hylleraas coordinates. We have used this method to calculate the ground state energy of a two-electron system in a spherical dot and a disk-like quantum dot separately. Under parabolic confinement potential and within effective mass approximation size and shape effects of quantum dots on the ground state energy of two electrons have been investigated. The calculation shows that our results even with a small number of basis states are in good agreement with previous theoretical results.展开更多
Image segmentation is a hot topic in image science. In this paper we present a new variational segmentation model based on the theory of Mumford-Shah model. The aim of our model is to divide noised image, according to...Image segmentation is a hot topic in image science. In this paper we present a new variational segmentation model based on the theory of Mumford-Shah model. The aim of our model is to divide noised image, according to a certain criterion, into homogeneous and smooth regions that should correspond to structural units in the scene or objects of interest. The proposed region-based model uses total variation as a regularization term, and different fidelity term can be used for image segmentation in the cases of physical noise, such as Gaussian, Poisson and multiplicative speckle noise. Our model consists of five weighted terms, two of them are responsible for image denoising based on fidelity term and total variation term, the others assure that the three conditions of adherence to the data, smoothing, and discontinuity detection are met at once. We also develop a primal-dual hybrid gradient algorithm for our model. Numerical results on various synthetic and real images are provided to compare our method with others, these results show that our proposed model and algorithms are effective.展开更多
This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification o...This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.展开更多
文摘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.
基金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.
基金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.
文摘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 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 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.
文摘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.
文摘The variational method is applied to calculate the dispersion characteristics of disc-loaded waveguide slow-wave structures. The parameters describing the waveguide discontinuities in disc-loaded waveguide are calculated by the variational method. Then the dispersion characteristics of slow-wave structures are obtained using lossless microwave quadrupole theory. Good agreement was observed between results of the Variational method and those of field matching method and high frequency structure simulator. In the case of broad band, results of the variational method are better than those of field matching method.
基金Project supported by the Cumhuriyet University National MOVPE Crystal Growth and Characterization Laboratory,DPT-K120,TUBITAK (Grant Nos TBAG 105T492,TBAG 107T012,and TBAG-108T015)
文摘In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.
基金supported by the Cumhuriyet University,National MOVPE Crystal Growth and Characterization Laboratory,Turkey,DPT-K 120,TUBITAK (Grant Nos TBAG 105T492,TBAG 107T012,and TBAG-108T015)
文摘In this work, we propose an efficient method of reducing the computational effort of variational calculation with a Hylleraas-like trial wavefunction. The method consists of introducing integral transforms for the terms as r12^κ exp (-λτ12) which provide the calculation of the expectation value of energy and the relevant matrix elements to be done analytically over single-electron coordinates instead of Hylleraas coordinates. We have used this method to calculate the ground state energy of a two-electron system in a spherical dot and a disk-like quantum dot separately. Under parabolic confinement potential and within effective mass approximation size and shape effects of quantum dots on the ground state energy of two electrons have been investigated. The calculation shows that our results even with a small number of basis states are in good agreement with previous theoretical results.
基金Supported in part by the NNSF of China(11301129,11271323,91330105,11326033)the Zhejiang Provincial Natural Science Foundation of China(LQ13A010025,LZ13A010002)
文摘Image segmentation is a hot topic in image science. In this paper we present a new variational segmentation model based on the theory of Mumford-Shah model. The aim of our model is to divide noised image, according to a certain criterion, into homogeneous and smooth regions that should correspond to structural units in the scene or objects of interest. The proposed region-based model uses total variation as a regularization term, and different fidelity term can be used for image segmentation in the cases of physical noise, such as Gaussian, Poisson and multiplicative speckle noise. Our model consists of five weighted terms, two of them are responsible for image denoising based on fidelity term and total variation term, the others assure that the three conditions of adherence to the data, smoothing, and discontinuity detection are met at once. We also develop a primal-dual hybrid gradient algorithm for our model. Numerical results on various synthetic and real images are provided to compare our method with others, these results show that our proposed model and algorithms are effective.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10771019 and 10826107)
文摘This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.