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.展开更多
Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface tem...Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface temperature(SST),sea surface height anomaly(SSHA),and sea surface salinity(SSS).This study employs a variational method to reconstruct the three-dimensional thermohaline structure of the Arctic Ocean.Compared to the Regional Arctic Reanalysis(RARE),the reconstruction well captures both the horizontal and vertical temperature and salinity structures in the Arctic.It demonstrates superior skill over RARE,when compared with Argo profiles and Ice-Tethered Profiler(ITP)observations.The reconstruction is particularly effective in ice-covered regions,where it more accurately captures the transition from Pacific water to Atlantic water compared to RARE.These findings underscore the potential of applying Arctic satellite data to reconstruct vertical thermohaline structures in the Arctic,particularly in areas due to lack of the subsurface observation reanalysis data exhibit significant biases.As Arctic satellite observations continue to advance,the applications of this method are becoming increasingly promising,which is useful for monitoring the ice-covered region environment and can be applied to oceanographic research.展开更多
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.展开更多
In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseu...In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].展开更多
Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational pr...Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational principles with multi-variables without arty variational crisis phenomenon. The method is to construct an energy trial-functional with an unknown function F, which can be readily identified by making the trial-functional stationary and using known constraint equations. As a result generalized variational principles with two kinds of independent variables (such as well-known Hellinger-Reissner variational principle and Hu-Washizu principle) and generalized variational principles with three kinds of independent variables (such as Chien's generalized variational principles) in elasticity have been deduced without using Lagrange multiplier method. By semi-inverse method, the author has also proved that Hu-Washizu principle is actually a variational principle with only two kinds of independent variables, stress-strain relations are still its constraints.展开更多
In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight env...In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight envelop.First,a coupled nonlinear six-degree-of-freedom equation of motion for a guided spinning rocket is developed,and the lateral acceleration motion is modeled as a control plant with time-varying matched uncertainties and unmodeled dynamics.Then,a robust adaptive control method is proposed by combining Bregman divergence and variational method to achieve fast adaption and maintain bounded tracking.The stability of the resulting closed-loop system is proved,and the ultimate bound and convergence rate are also analyzed.Finally,numerical simulations are performed for a single operating point and the whole flight trajectory to show the robustness and adaptability of the proposed method with respect to timevarying uncertainties and unmodeled dynamics.展开更多
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.展开更多
Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving m...Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense.As the state-of-the-art moving mesh method,the variational mesh adaptation approach has been introduced to CFD calculation.However,quickly estimating the flow field on the updated meshes during the iterative algorithm is challenging.A mesh optimization method,which embeds a machine learning regression model into the variational mesh adaptation,is proposed.The regression model captures the mapping between the initial mesh nodes and the flow field,so that the variational method could move mesh nodes iteratively by solving the mesh functional which is built from the estimated flow field on the updated mesh via the regression model.After the optimization,the density of the nodes in the high gradient area increases while the density in the low gradient area decreases.Benchmark examples are first used to verify the feasibility and effectiveness of the proposed method.And then we use the steady subsonic and transonic flows over cylinder and NACA0012 airfoil on unstructured triangular meshes to test our method.Results show that the proposed method significantly improves the accuracy of the local flow features on the adaptive meshes.Our work indicates that the proposed mesh optimization approach is promising for improving the accuracy and efficiency of CFD computation.展开更多
The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisi...The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.展开更多
Simulations of contact problems involving at least one plastic solid may be costly due to their strong nonlinearity and requirements of stability.In this work,we develop an explicit asynchronous variational integrator...Simulations of contact problems involving at least one plastic solid may be costly due to their strong nonlinearity and requirements of stability.In this work,we develop an explicit asynchronous variational integrator(AVI)for inelastic non-frictional contact problems involving a plastic solid.The AVI assigns each element in the mesh an independent time step and updates the solution at the elements and nodes asynchronously.This asynchrony makes the AVI highly efficient in solving such bi-material problems.Taking advantage of the AVI,the constitutive update is locally performed in one element at a time,and contact constraints are also enforced on only one element.The time step of the contact element is subdivided into multiple segments,and the fields are updated accordingly.During a contact event,only one element involving a few degrees of freedom is considered,leading to high efficiency.The proposed formulation is first verified with a pure elastodynamics benchmark and further applied to a contact problem involving an elastoplastic solid with non-associative volumetric hardening.The numerical results indicate that the AVI exhibits excellent energy behaviors and has high computational efficiency.展开更多
The main purpose of this paper is to introduce and deal with a self adaptive inertial subgradient extragradient iterative algorithm with a new and interesting stepsize rule in real Hilbert spaces.Under some proper con...The main purpose of this paper is to introduce and deal with a self adaptive inertial subgradient extragradient iterative algorithm with a new and interesting stepsize rule in real Hilbert spaces.Under some proper control conditions imposed on the coefficients and operators,we prove a new strong convergence result for solving variational inequalities with regard to pseudomonotone and Lipschitzian operators.Moreover,some numerical simulation results are given to show the rationality and validity of our algorithm.展开更多
Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new metho...Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new method possesses a high efficiency. insolving such applied engineering problems as obstacle problems and .free boundary.problems arising in fluid lubrications.展开更多
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.展开更多
The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods fo...The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.展开更多
In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient ext...In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient extragradient method and the projection contraction method.Some very recent papers have considered different inertial algorithms which allowed the inertial factor is chosen in[0;1].The purpose of this work is to continue working in this direction,we propose another inertial subgradient extragradient method that the inertial factor can be chosen in a special case to be 1.Under suitable mild conditions,we establish the weak convergence of the proposed algorithm.Moreover,linear convergence is obtained under strong pseudomonotonicity and Lipschitz continuity assumptions.Finally,some numerical illustrations are given to confirm the theoretical analysis.展开更多
The generalized variational data assimilation for non-differential dynamical systems is studied.There is no tangent linear model for non-differential systems and thus the general adjoint model can not be derived in th...The generalized variational data assimilation for non-differential dynamical systems is studied.There is no tangent linear model for non-differential systems and thus the general adjoint model can not be derived in the traditional way.The weak form of the original system was introduced, and then the generalized adjoint model was derived. The generalized variational data assimilation methods were developed for non-differential low dimensional system and non-differential high dimensional system with global and local observations. Furthermore, ideas in inverse problems are introduced to 4DVAR (Four-dimensional variational) of non-differential partial differential system with local observations.展开更多
This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODE...This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.展开更多
The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational in...The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational inequality problem based on a new smoothing approximation function.The proposed algorithm is proved to be well defined and convergent globally under weaker conditions.展开更多
基金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.
基金The National Key R&D Program of China under contract No.2022YFE0106400the China Scholarship Council under contract No.202206710071+2 种基金the Postgraduate Research&Practice Innovation Program of Jiangsu Province under contract No.KYCX23_0657the Special Founds for Creative Research under contract No.2022C61540the Opening Project of the Key Laboratory of Marine Environmental Information Technology under contract No.521037412.
文摘Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface temperature(SST),sea surface height anomaly(SSHA),and sea surface salinity(SSS).This study employs a variational method to reconstruct the three-dimensional thermohaline structure of the Arctic Ocean.Compared to the Regional Arctic Reanalysis(RARE),the reconstruction well captures both the horizontal and vertical temperature and salinity structures in the Arctic.It demonstrates superior skill over RARE,when compared with Argo profiles and Ice-Tethered Profiler(ITP)observations.The reconstruction is particularly effective in ice-covered regions,where it more accurately captures the transition from Pacific water to Atlantic water compared to RARE.These findings underscore the potential of applying Arctic satellite data to reconstruct vertical thermohaline structures in the Arctic,particularly in areas due to lack of the subsurface observation reanalysis data exhibit significant biases.As Arctic satellite observations continue to advance,the applications of this method are becoming increasingly promising,which is useful for monitoring the ice-covered region environment and can be applied to oceanographic research.
基金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.
基金supported by Scientific Research Fund of Sichuan Provincial Education Department (09ZB102)Scientific Research Fund of Science and Technology Deportment of Sichuan Provincial (2011JYZ011)
文摘In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].
文摘Semi-inverse method, which is an integration and an extension of Hu's try-and-error method, Chien's veighted residual method and Liu's systematic method, is proposed to establish generalized variational principles with multi-variables without arty variational crisis phenomenon. The method is to construct an energy trial-functional with an unknown function F, which can be readily identified by making the trial-functional stationary and using known constraint equations. As a result generalized variational principles with two kinds of independent variables (such as well-known Hellinger-Reissner variational principle and Hu-Washizu principle) and generalized variational principles with three kinds of independent variables (such as Chien's generalized variational principles) in elasticity have been deduced without using Lagrange multiplier method. By semi-inverse method, the author has also proved that Hu-Washizu principle is actually a variational principle with only two kinds of independent variables, stress-strain relations are still its constraints.
基金supported by the National Natural Science Foundation of China (No. 11532002)。
文摘In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight envelop.First,a coupled nonlinear six-degree-of-freedom equation of motion for a guided spinning rocket is developed,and the lateral acceleration motion is modeled as a control plant with time-varying matched uncertainties and unmodeled dynamics.Then,a robust adaptive control method is proposed by combining Bregman divergence and variational method to achieve fast adaption and maintain bounded tracking.The stability of the resulting closed-loop system is proved,and the ultimate bound and convergence rate are also analyzed.Finally,numerical simulations are performed for a single operating point and the whole flight trajectory to show the robustness and adaptability of the proposed method with respect to timevarying uncertainties and unmodeled dynamics.
文摘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.
基金co-supported by the Key Laboratory of Aerodynamic Noise Control,China(No.ANCL20190103)the State Key Laboratory of Aerodynamics,China(No.SKLA20180102)the Aeronautical Science Foundation of China(Nos.2018ZA52002 and 2019ZA052011)。
文摘Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense.As the state-of-the-art moving mesh method,the variational mesh adaptation approach has been introduced to CFD calculation.However,quickly estimating the flow field on the updated meshes during the iterative algorithm is challenging.A mesh optimization method,which embeds a machine learning regression model into the variational mesh adaptation,is proposed.The regression model captures the mapping between the initial mesh nodes and the flow field,so that the variational method could move mesh nodes iteratively by solving the mesh functional which is built from the estimated flow field on the updated mesh via the regression model.After the optimization,the density of the nodes in the high gradient area increases while the density in the low gradient area decreases.Benchmark examples are first used to verify the feasibility and effectiveness of the proposed method.And then we use the steady subsonic and transonic flows over cylinder and NACA0012 airfoil on unstructured triangular meshes to test our method.Results show that the proposed method significantly improves the accuracy of the local flow features on the adaptive meshes.Our work indicates that the proposed mesh optimization approach is promising for improving the accuracy and efficiency of CFD computation.
文摘The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.
基金support from the Hui-Chun Chin and Tsung-Dao Lee Chinese Undergraduate Research Endowment(CURE).
文摘Simulations of contact problems involving at least one plastic solid may be costly due to their strong nonlinearity and requirements of stability.In this work,we develop an explicit asynchronous variational integrator(AVI)for inelastic non-frictional contact problems involving a plastic solid.The AVI assigns each element in the mesh an independent time step and updates the solution at the elements and nodes asynchronously.This asynchrony makes the AVI highly efficient in solving such bi-material problems.Taking advantage of the AVI,the constitutive update is locally performed in one element at a time,and contact constraints are also enforced on only one element.The time step of the contact element is subdivided into multiple segments,and the fields are updated accordingly.During a contact event,only one element involving a few degrees of freedom is considered,leading to high efficiency.The proposed formulation is first verified with a pure elastodynamics benchmark and further applied to a contact problem involving an elastoplastic solid with non-associative volumetric hardening.The numerical results indicate that the AVI exhibits excellent energy behaviors and has high computational efficiency.
文摘The main purpose of this paper is to introduce and deal with a self adaptive inertial subgradient extragradient iterative algorithm with a new and interesting stepsize rule in real Hilbert spaces.Under some proper control conditions imposed on the coefficients and operators,we prove a new strong convergence result for solving variational inequalities with regard to pseudomonotone and Lipschitzian operators.Moreover,some numerical simulation results are given to show the rationality and validity of our algorithm.
文摘Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new method possesses a high efficiency. insolving such applied engineering problems as obstacle problems and .free boundary.problems arising in fluid lubrications.
基金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.
基金supported by the Major Research Plan on West-Pacific Earth System Multispheric Interactions (Nos.91858215,91958206)the National Natural Science Foundation of China (NSFC)Shiptime Sharing Project (No.41949581)the Key Research and Development Program of Shandong Province (No.2019GHY112019)。
文摘The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.
基金funded by the University of Science,Vietnam National University,Hanoi under project number TN.21.01。
文摘In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient extragradient method and the projection contraction method.Some very recent papers have considered different inertial algorithms which allowed the inertial factor is chosen in[0;1].The purpose of this work is to continue working in this direction,we propose another inertial subgradient extragradient method that the inertial factor can be chosen in a special case to be 1.Under suitable mild conditions,we establish the weak convergence of the proposed algorithm.Moreover,linear convergence is obtained under strong pseudomonotonicity and Lipschitz continuity assumptions.Finally,some numerical illustrations are given to confirm the theoretical analysis.
文摘The generalized variational data assimilation for non-differential dynamical systems is studied.There is no tangent linear model for non-differential systems and thus the general adjoint model can not be derived in the traditional way.The weak form of the original system was introduced, and then the generalized adjoint model was derived. The generalized variational data assimilation methods were developed for non-differential low dimensional system and non-differential high dimensional system with global and local observations. Furthermore, ideas in inverse problems are introduced to 4DVAR (Four-dimensional variational) of non-differential partial differential system with local observations.
文摘This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.
基金Supported by the NNSF of China（11071041）Supported by the Fujian Natural Science Foundation（2009J01002）Supported by the Fujian Department of Education Foundation（JA11270）
文摘The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational inequality problem based on a new smoothing approximation function.The proposed algorithm is proved to be well defined and convergent globally under weaker conditions.
