Impulse noise removal is an important task in image restoration.In this paper,we introduce a general nonsmooth nonconvex model for recovering images degraded by blur and impulsive noise,which can easily include some p...Impulse noise removal is an important task in image restoration.In this paper,we introduce a general nonsmooth nonconvex model for recovering images degraded by blur and impulsive noise,which can easily include some prior information,such as box constraint or low rank,etc.To deal with the nonconvex problem,we employ the proximal linearized minimization algorithm.For the subproblem,we use the alternating direction method of multipliers to solve it.Furthermore,based on the assumption that the objective function satisfies the KurdykaLojasiewicz property,we prove the global convergence of the proposed algorithm.Numerical experiments demonstrate that our method outperforms both the l1TV and Nonconvex TV models in terms of subjective and objective quality measurements.展开更多
This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the...The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the band structures of TiO2 crystals doped with transition metal atoms become narrower. Interesting, an excursion towards high energy level with increasing atomic number in the same element period could be observed after doping with transition metal atoms.展开更多
The low-field nuclear magnetic resonance(NMR)technique has been used to probe the pore size distribution and the fluid composition in geophysical prospecting and related fields.However,the speed and accuracy of the ex...The low-field nuclear magnetic resonance(NMR)technique has been used to probe the pore size distribution and the fluid composition in geophysical prospecting and related fields.However,the speed and accuracy of the existing numerical inversion methods are still challenging due to the ill-posed nature of the first kind Fredholm integral equation and the contamination of the noises.This paper proposes a novel inversion algorithmto accelerate the convergence and enhance the precision using empirical truncated singular value decompositions(TSVD)and the linearized Bregman iteration.The L1 penalty term is applied to construct the objective function,and then the linearized Bregman iteration is utilized to obtain fast convergence.To reduce the complexity of the computation,empirical TSVD is proposed to compress the kernel matrix and determine the appropriate truncated position.This novel inversion method is validated using numerical simulations.The results indicate that the proposed novel method is significantly efficient and can achieve quick and effective data solutions with low signal-to-noise ratios.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
To numerically solve the initial-boundary value problem of the Korteweg-de Vries equation,an equivalent coupled system of nonlinear equations is obtained by the method of reduction of order.Then,a difference scheme is...To numerically solve the initial-boundary value problem of the Korteweg-de Vries equation,an equivalent coupled system of nonlinear equations is obtained by the method of reduction of order.Then,a difference scheme is constructed for the system.The new variable introduced can be separated from the difference scheme to obtain another difference scheme containing only the original variable.The energy method is applied to the theoretical analysis of the difference scheme.Results show that the difference scheme is uniquely solvable and satisfies the energy conservation law corresponding to the original problem.Moreover,the difference scheme converges when the step ratio satisfies a constraint condition,and the temporal and spatial convergence orders are both two.Numerical examples verify the convergence order and the invariant of the difference scheme.Furthermore,the step ratio constraint is unnecessary for the convergence of the difference scheme.Compared with a known two-level nonlinear difference scheme,the proposed difference scheme has more advantages in numerical calculation.展开更多
Steel catenary risers, (SCR) usually installed between seabed wellhead and floating platform are subjected to vortex shedding. These impose direct forces, hence cyclic stresses, and fatigue damage on the SCR. Riser fa...Steel catenary risers, (SCR) usually installed between seabed wellhead and floating platform are subjected to vortex shedding. These impose direct forces, hence cyclic stresses, and fatigue damage on the SCR. Riser failure has both economic and environmental consequences;hence the design life is usually greater than the field life, which is significantly reduced by vortex induced vibration (VIV). In this study, SCR and metOcean data from a field in Offshore Nigeria were substituted into linearized hydrodynamic models for simulations. The results showed that the hang off and touchdown regions were most susceptible to fatigue failure. Further analysis using Miner-Palm green models revealed that the fatigue life reduced from a design value of 20-years to 17.04-years, shortened by 2.96-years due to VIV. Furthermore, a maximum wave load of 5.154 kN was observed. The wave loads results corroborated with those obtained from finite element Orca Flex software, yielding a correlation coefficient of 0.975.展开更多
In order to design linear controller for nonlinear systems,a simple but efficient method of modeling a nonlinear system was proposed by means of multiple linearized models at different operating points in the entire r...In order to design linear controller for nonlinear systems,a simple but efficient method of modeling a nonlinear system was proposed by means of multiple linearized models at different operating points in the entire range of the expected changes of the operating points.The original nonlinear system was described by linear combination of these multiple linearized models,with the linear combination parameters being identified on line based on least squares method.Model Predictive Control,an optimization based technique,was used to design the linear controller.A sufficient condition for ensuring the existence of a linear controller for the original nonlinear system was also given.Good performance indicated by two simulated examples confirms the usefulness of the proposed method.展开更多
A double-input–multi-output linearized system is developed using the state-space method for dynamic analysis of methanation process of coke oven gas.The stability of reactor alone and reactor with feed-effluent heat ...A double-input–multi-output linearized system is developed using the state-space method for dynamic analysis of methanation process of coke oven gas.The stability of reactor alone and reactor with feed-effluent heat exchanger is compared through the dominant poles of the system transfer functions.With single or double disturbance of temperature and CO concentration at the reactor inlet,typical dynamic behavior in the reactor,including fast concentration response,slow temperature response and inverse response,is revealed for further understanding of the counteraction and synergy effects caused by simultaneous variation of concentration and temperature.Analysis results show that the stability of the reactor loop is more sensitive than that of reactor alone due to the positive heat feedback.Remarkably,with the decrease of heat exchange efficiency,the reactor system may display limit cycle behavior for a pair of complex conjugate poles across the imaginary axis.展开更多
In this paper a linearized and unified yield crierion of metals is presented, which is in a form of a set of linear functions with two pararneters. The parameters are ex- pressed in terms of tension yield stress and ...In this paper a linearized and unified yield crierion of metals is presented, which is in a form of a set of linear functions with two pararneters. The parameters are ex- pressed in terms of tension yield stress and so-called “shear-stretch ratio” and can bereadily determined from experimental data. It is shown that in stress space the set of yield functions is a set of polygons with twelve edges located between the Tresca’s hexagon and twin-shear-stress hexagon￣[1]. In this paper the present yield function isused to analyse the prestressiap loose running fit cylinders.展开更多
The explicit solution to Cauchy problem for linearized system of two-dimensional isentropic flow with axisymmetrical initial data in gas dynamics is given.
The Alternating Direction Multiplier Method (ADMM) is widely used in various fields, and different variables are customized in the literature for different application scenarios [1] [2] [3] [4]. Among them, the linear...The Alternating Direction Multiplier Method (ADMM) is widely used in various fields, and different variables are customized in the literature for different application scenarios [1] [2] [3] [4]. Among them, the linearized alternating direction multiplier method (LADMM) has received extensive attention because of its effectiveness and ease of implementation. This paper mainly discusses the application of ADMM in dictionary learning (non-convex problem). Many numerical experiments show that to achieve higher convergence accuracy, the convergence speed of ADMM is slower, especially near the optimal solution. Therefore, we introduce the linearized alternating direction multiplier method (LADMM) to accelerate the convergence speed of ADMM. Specifically, the problem is solved by linearizing the quadratic term of the subproblem, and the convergence of the algorithm is proved. Finally, there is a brief summary of the full text.展开更多
The paper is concerned with the problem of reduction of the general relativity theory to the Newton gravitation theory for a gravitation field with relatively low intensity. This problem is traditionally solved on the...The paper is concerned with the problem of reduction of the general relativity theory to the Newton gravitation theory for a gravitation field with relatively low intensity. This problem is traditionally solved on the basis of linearized equations of general relativity which, being matched to the Newton theory equations, allow us to link the classical gravitation constant with the constant entering the general relativity equations. Analysis of the linearized general relativity equations shows that it can be done only for empty space in which the energy tensor is zero. In solids, the set of linearized general relativity equations is not consistent and is not reduced to the Newton theory equations. Specific features of the problem are demonstrated with the spherically symmetric static problem of general relativity which has the closed-form solution.展开更多
Laboratory tests were carried out to study the breakage kinetics of diasporic bauxite and determine its breakage distribution function. Non-first order breakage with different deceleration rates for different size int...Laboratory tests were carried out to study the breakage kinetics of diasporic bauxite and determine its breakage distribution function. Non-first order breakage with different deceleration rates for different size intervals is found, which is most probably caused by the heterogeneity of the ore. Piecewise linearization method is proposed to describe the non-first order breakage according to its characteristics. In the method, grinding time is divided into several intervals and breakage is assumed to be first order in each interval. So, the breakage rates are calculated by taking the product of the last interval as feed and then established as a function of particle size and grinding time. Based on the predetermined breakage rate function, the breakage distribution of the ore is back-calculated from the experimental data using the population balance model (PBM). Finally, the obtained breakage parameters are validated and the simulated data are in good agreement with the experimental data. The obtained breakage distribution and the method for breakage rate description are both significant for modeling the full scale ball milling process of bauxite.展开更多
A linearized rock physics inversion method is proposed to deal with two important issues, rock physical model and inversion algorithm, which restrict the accuracy of rock physics inversion. In this method, first, the ...A linearized rock physics inversion method is proposed to deal with two important issues, rock physical model and inversion algorithm, which restrict the accuracy of rock physics inversion. In this method, first, the complex rock physics model is expanded into Taylor series to get the first-order approximate expression of the inverse problem of rock physics;then the damped least square method is used to solve the linearized rock physics inverse problem directly to get the analytical solution of the rock physics inverse problem. This method does not need global optimization or random sampling, but directly calculates the inverse operation, with high computational efficiency. The theoretical model analysis shows that the linearized rock physical model can be used to approximate the complex rock physics model. The application of actual logging data and seismic data shows that the linearized rock physics inversion method can obtain accurate physical parameters. This method is suitable for linear or slightly non-linear rock physics model, but may not be suitable for highly non-linear rock physics model.展开更多
This paper focuses on the optimal error analysis of a linearized Crank-Nicolson finite element scheme for the time-dependent penetrative convection problem,where the mini element and piecewise linear finite element ar...This paper focuses on the optimal error analysis of a linearized Crank-Nicolson finite element scheme for the time-dependent penetrative convection problem,where the mini element and piecewise linear finite element are used to approximate the velocity field,the pressure,and the temperature,respectively.We proved that the proposed finite element scheme is unconditionally stable and the optimal error estimates in\(L^2\)-norm are derived.Finally,numerical results are presented to confirm the theoretical analysis.展开更多
Integrated high-linearity modulators are crucial for high dynamic-range microwave photonic(MWP)systems.Conventional linearization schemes usually involve the fine tuning of radio-frequency(RF)power distribution,which ...Integrated high-linearity modulators are crucial for high dynamic-range microwave photonic(MWP)systems.Conventional linearization schemes usually involve the fine tuning of radio-frequency(RF)power distribution,which is rather inconvenient for practical applications and can hardly be implemented on the integrated photonics chip.In this paper,we propose an elegant scheme to linearize a silicon-based modulator in which the active tuning of RF power is eliminated.The device consists of two carrier-depletion-based Mach-Zehnder modulators(MZMs),which are connected in series by a 1×2 thermal optical switch(OS).The OS is used to adjust the ratio between the modulation depths of the two sub-MZMs.Under a proper ratio,the complementary third-order intermodulation distortion(IMD3)of the two sub-MZMs can effectively cancel each other out.The measured spurious-free dynamic ranges for IMD3 are 131,127,118,110,and 109 d B·Hz^(6∕7)at frequencies of 1,10,20,30,and 40 GHz,respectively,which represent the highest linearities ever reached by the integrated modulator chips on all available material platforms.展开更多
We consider a wide range of non-convex regularized minimization problems, where the non-convex regularization term is composite with a linear function engaged in sparse learning. Recent theoretical investigations have...We consider a wide range of non-convex regularized minimization problems, where the non-convex regularization term is composite with a linear function engaged in sparse learning. Recent theoretical investigations have demonstrated their superiority over their convex counterparts. The computational challenge lies in the fact that the proximal mapping associated with non-convex regularization is not easily obtained due to the imposed linear composition. Fortunately, the problem structure allows one to introduce an auxiliary variable and reformulate it as an optimization problem with linear constraints, which can be solved using the Linearized Alternating Direction Method of Multipliers （LADMM）. Despite the success of LADMM in practice, it remains unknown whether LADMM is convergent in solving such non-convex compositely regularized optimizations. In this research, we first present a detailed convergence analysis of the LADMM algorithm for solving a non-convex compositely regularized optimization problem with a large class of non-convex penalties. Furthermore, we propose an Adaptive LADMM （AdaLADMM） algorithm with a line-search criterion. Experimental results on different genres of datasets validate the efficacy of the proposed algorithm.展开更多
The proximal alternating linearized minimization(PALM)method suits well for solving blockstructured optimization problems,which are ubiquitous in real applications.In the cases where subproblems do not have closed-for...The proximal alternating linearized minimization(PALM)method suits well for solving blockstructured optimization problems,which are ubiquitous in real applications.In the cases where subproblems do not have closed-form solutions,e.g.,due to complex constraints,infeasible subsolvers are indispensable,giving rise to an infeasible inexact PALM(PALM-I).Numerous efforts have been devoted to analyzing the feasible PALM,while little attention has been paid to the PALM-I.The usage of the PALM-I thus lacks a theoretical guarantee.The essential difficulty of analysis consists in the objective value nonmonotonicity induced by the infeasibility.We study in the present work the convergence properties of the PALM-I.In particular,we construct a surrogate sequence to surmount the nonmonotonicity issue and devise an implementable inexact criterion.Based upon these,we manage to establish the stationarity of any accumulation point,and moreover,show the iterate convergence and the asymptotic convergence rates under the assumption of the Lojasiewicz property.The prominent advantages of the PALM-I on CPU time are illustrated via numerical experiments on problems arising from quantum physics and 3-dimensional anisotropic frictional contact.展开更多
基金Supported by the National Natural Science Foundations of China(Grant No.12061045,12031003)the Guangzhou Education Scientific Research Project 2024(Grant No.202315829)the Natural Science Foundation of Jiangxi Province(Grant No.20224ACB211004)。
文摘Impulse noise removal is an important task in image restoration.In this paper,we introduce a general nonsmooth nonconvex model for recovering images degraded by blur and impulsive noise,which can easily include some prior information,such as box constraint or low rank,etc.To deal with the nonconvex problem,we employ the proximal linearized minimization algorithm.For the subproblem,we use the alternating direction method of multipliers to solve it.Furthermore,based on the assumption that the objective function satisfies the KurdykaLojasiewicz property,we prove the global convergence of the proposed algorithm.Numerical experiments demonstrate that our method outperforms both the l1TV and Nonconvex TV models in terms of subjective and objective quality measurements.
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
文摘The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the band structures of TiO2 crystals doped with transition metal atoms become narrower. Interesting, an excursion towards high energy level with increasing atomic number in the same element period could be observed after doping with transition metal atoms.
基金support by the National Nature Science Foundation of China(42174142)CNPC Innovation Found(2021DQ02-0402)National Key Foundation for Exploring Scientific Instrument of China(2013YQ170463).
文摘The low-field nuclear magnetic resonance(NMR)technique has been used to probe the pore size distribution and the fluid composition in geophysical prospecting and related fields.However,the speed and accuracy of the existing numerical inversion methods are still challenging due to the ill-posed nature of the first kind Fredholm integral equation and the contamination of the noises.This paper proposes a novel inversion algorithmto accelerate the convergence and enhance the precision using empirical truncated singular value decompositions(TSVD)and the linearized Bregman iteration.The L1 penalty term is applied to construct the objective function,and then the linearized Bregman iteration is utilized to obtain fast convergence.To reduce the complexity of the computation,empirical TSVD is proposed to compress the kernel matrix and determine the appropriate truncated position.This novel inversion method is validated using numerical simulations.The results indicate that the proposed novel method is significantly efficient and can achieve quick and effective data solutions with low signal-to-noise ratios.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金The National Natural Science Foundation of China(No.11671081).
文摘To numerically solve the initial-boundary value problem of the Korteweg-de Vries equation,an equivalent coupled system of nonlinear equations is obtained by the method of reduction of order.Then,a difference scheme is constructed for the system.The new variable introduced can be separated from the difference scheme to obtain another difference scheme containing only the original variable.The energy method is applied to the theoretical analysis of the difference scheme.Results show that the difference scheme is uniquely solvable and satisfies the energy conservation law corresponding to the original problem.Moreover,the difference scheme converges when the step ratio satisfies a constraint condition,and the temporal and spatial convergence orders are both two.Numerical examples verify the convergence order and the invariant of the difference scheme.Furthermore,the step ratio constraint is unnecessary for the convergence of the difference scheme.Compared with a known two-level nonlinear difference scheme,the proposed difference scheme has more advantages in numerical calculation.
文摘Steel catenary risers, (SCR) usually installed between seabed wellhead and floating platform are subjected to vortex shedding. These impose direct forces, hence cyclic stresses, and fatigue damage on the SCR. Riser failure has both economic and environmental consequences;hence the design life is usually greater than the field life, which is significantly reduced by vortex induced vibration (VIV). In this study, SCR and metOcean data from a field in Offshore Nigeria were substituted into linearized hydrodynamic models for simulations. The results showed that the hang off and touchdown regions were most susceptible to fatigue failure. Further analysis using Miner-Palm green models revealed that the fatigue life reduced from a design value of 20-years to 17.04-years, shortened by 2.96-years due to VIV. Furthermore, a maximum wave load of 5.154 kN was observed. The wave loads results corroborated with those obtained from finite element Orca Flex software, yielding a correlation coefficient of 0.975.
文摘In order to design linear controller for nonlinear systems,a simple but efficient method of modeling a nonlinear system was proposed by means of multiple linearized models at different operating points in the entire range of the expected changes of the operating points.The original nonlinear system was described by linear combination of these multiple linearized models,with the linear combination parameters being identified on line based on least squares method.Model Predictive Control,an optimization based technique,was used to design the linear controller.A sufficient condition for ensuring the existence of a linear controller for the original nonlinear system was also given.Good performance indicated by two simulated examples confirms the usefulness of the proposed method.
基金Supported by the Major Research plan of the National Natural Science Foundation of China(91334101)the National Basic Research Program of China(2009CB219906)the National Natural Science Foundation of China(21276203)
文摘A double-input–multi-output linearized system is developed using the state-space method for dynamic analysis of methanation process of coke oven gas.The stability of reactor alone and reactor with feed-effluent heat exchanger is compared through the dominant poles of the system transfer functions.With single or double disturbance of temperature and CO concentration at the reactor inlet,typical dynamic behavior in the reactor,including fast concentration response,slow temperature response and inverse response,is revealed for further understanding of the counteraction and synergy effects caused by simultaneous variation of concentration and temperature.Analysis results show that the stability of the reactor loop is more sensitive than that of reactor alone due to the positive heat feedback.Remarkably,with the decrease of heat exchange efficiency,the reactor system may display limit cycle behavior for a pair of complex conjugate poles across the imaginary axis.
文摘In this paper a linearized and unified yield crierion of metals is presented, which is in a form of a set of linear functions with two pararneters. The parameters are ex- pressed in terms of tension yield stress and so-called “shear-stretch ratio” and can bereadily determined from experimental data. It is shown that in stress space the set of yield functions is a set of polygons with twelve edges located between the Tresca’s hexagon and twin-shear-stress hexagon￣[1]. In this paper the present yield function isused to analyse the prestressiap loose running fit cylinders.
文摘The explicit solution to Cauchy problem for linearized system of two-dimensional isentropic flow with axisymmetrical initial data in gas dynamics is given.
文摘The Alternating Direction Multiplier Method (ADMM) is widely used in various fields, and different variables are customized in the literature for different application scenarios [1] [2] [3] [4]. Among them, the linearized alternating direction multiplier method (LADMM) has received extensive attention because of its effectiveness and ease of implementation. This paper mainly discusses the application of ADMM in dictionary learning (non-convex problem). Many numerical experiments show that to achieve higher convergence accuracy, the convergence speed of ADMM is slower, especially near the optimal solution. Therefore, we introduce the linearized alternating direction multiplier method (LADMM) to accelerate the convergence speed of ADMM. Specifically, the problem is solved by linearizing the quadratic term of the subproblem, and the convergence of the algorithm is proved. Finally, there is a brief summary of the full text.
文摘The paper is concerned with the problem of reduction of the general relativity theory to the Newton gravitation theory for a gravitation field with relatively low intensity. This problem is traditionally solved on the basis of linearized equations of general relativity which, being matched to the Newton theory equations, allow us to link the classical gravitation constant with the constant entering the general relativity equations. Analysis of the linearized general relativity equations shows that it can be done only for empty space in which the energy tensor is zero. In solids, the set of linearized general relativity equations is not consistent and is not reduced to the Newton theory equations. Specific features of the problem are demonstrated with the spherically symmetric static problem of general relativity which has the closed-form solution.
基金Supported by the Fundamental Research Funds for the Central Universities (2012QNZT069)the Postdoctoral Science Foundation of China (2012M521413)+1 种基金the National Science Fund for Distinguished Young Scholars of China (61025015)the National Natural Science Foundation of China (61273187, 61273159)
文摘Laboratory tests were carried out to study the breakage kinetics of diasporic bauxite and determine its breakage distribution function. Non-first order breakage with different deceleration rates for different size intervals is found, which is most probably caused by the heterogeneity of the ore. Piecewise linearization method is proposed to describe the non-first order breakage according to its characteristics. In the method, grinding time is divided into several intervals and breakage is assumed to be first order in each interval. So, the breakage rates are calculated by taking the product of the last interval as feed and then established as a function of particle size and grinding time. Based on the predetermined breakage rate function, the breakage distribution of the ore is back-calculated from the experimental data using the population balance model (PBM). Finally, the obtained breakage parameters are validated and the simulated data are in good agreement with the experimental data. The obtained breakage distribution and the method for breakage rate description are both significant for modeling the full scale ball milling process of bauxite.
基金Supported by the China National Science and Technology Major Project(2017ZX05049-002,2016ZX05027004-001)the National Natural Science Foundation of China(41874146,41674130)+2 种基金Fundamental Research Funds for the Central University(18CX02061A)Innovative Fund Project of China National Petroleum Corporation(2016D-5007-0301)Scientific Research&Technology Development Project of China National Petroleum Corporation(2017D-3504).
文摘A linearized rock physics inversion method is proposed to deal with two important issues, rock physical model and inversion algorithm, which restrict the accuracy of rock physics inversion. In this method, first, the complex rock physics model is expanded into Taylor series to get the first-order approximate expression of the inverse problem of rock physics;then the damped least square method is used to solve the linearized rock physics inverse problem directly to get the analytical solution of the rock physics inverse problem. This method does not need global optimization or random sampling, but directly calculates the inverse operation, with high computational efficiency. The theoretical model analysis shows that the linearized rock physical model can be used to approximate the complex rock physics model. The application of actual logging data and seismic data shows that the linearized rock physics inversion method can obtain accurate physical parameters. This method is suitable for linear or slightly non-linear rock physics model, but may not be suitable for highly non-linear rock physics model.
基金supported by the National Natural Science Foundation of China(No.11771337)the Natural Science Foundation of Zhejiang Province of China(No.LY23A010002).
文摘This paper focuses on the optimal error analysis of a linearized Crank-Nicolson finite element scheme for the time-dependent penetrative convection problem,where the mini element and piecewise linear finite element are used to approximate the velocity field,the pressure,and the temperature,respectively.We proved that the proposed finite element scheme is unconditionally stable and the optimal error estimates in\(L^2\)-norm are derived.Finally,numerical results are presented to confirm the theoretical analysis.
基金National Natural Science Foundation of China(62305303,62205164)Science and Technology Plan Project of Zhejiang(2022C01108)。
文摘Integrated high-linearity modulators are crucial for high dynamic-range microwave photonic(MWP)systems.Conventional linearization schemes usually involve the fine tuning of radio-frequency(RF)power distribution,which is rather inconvenient for practical applications and can hardly be implemented on the integrated photonics chip.In this paper,we propose an elegant scheme to linearize a silicon-based modulator in which the active tuning of RF power is eliminated.The device consists of two carrier-depletion-based Mach-Zehnder modulators(MZMs),which are connected in series by a 1×2 thermal optical switch(OS).The OS is used to adjust the ratio between the modulation depths of the two sub-MZMs.Under a proper ratio,the complementary third-order intermodulation distortion(IMD3)of the two sub-MZMs can effectively cancel each other out.The measured spurious-free dynamic ranges for IMD3 are 131,127,118,110,and 109 d B·Hz^(6∕7)at frequencies of 1,10,20,30,and 40 GHz,respectively,which represent the highest linearities ever reached by the integrated modulator chips on all available material platforms.
基金supported by the National Natural Science Foundation of China(Nos.61303264,61202482,and 61202488)Guangxi Cooperative Innovation Center of Cloud Computing and Big Data(No.YD16505)Distinguished Young Scientist Promotion of National University of Defense Technology
文摘We consider a wide range of non-convex regularized minimization problems, where the non-convex regularization term is composite with a linear function engaged in sparse learning. Recent theoretical investigations have demonstrated their superiority over their convex counterparts. The computational challenge lies in the fact that the proximal mapping associated with non-convex regularization is not easily obtained due to the imposed linear composition. Fortunately, the problem structure allows one to introduce an auxiliary variable and reformulate it as an optimization problem with linear constraints, which can be solved using the Linearized Alternating Direction Method of Multipliers （LADMM）. Despite the success of LADMM in practice, it remains unknown whether LADMM is convergent in solving such non-convex compositely regularized optimizations. In this research, we first present a detailed convergence analysis of the LADMM algorithm for solving a non-convex compositely regularized optimization problem with a large class of non-convex penalties. Furthermore, we propose an Adaptive LADMM （AdaLADMM） algorithm with a line-search criterion. Experimental results on different genres of datasets validate the efficacy of the proposed algorithm.
基金supported by National Natural Science Foundation of China(Grant Nos.12125108,11971466,11991021,11991020,12021001 and 12288201)Key Research Program of Frontier Sciences,Chinese Academy of Sciences(Grant No.ZDBS-LY-7022)CAS(the Chinese Academy of Sciences)AMSS(Academy of Mathematics and Systems Science)-PolyU(The Hong Kong Polytechnic University)Joint Laboratory of Applied Mathematics.
文摘The proximal alternating linearized minimization(PALM)method suits well for solving blockstructured optimization problems,which are ubiquitous in real applications.In the cases where subproblems do not have closed-form solutions,e.g.,due to complex constraints,infeasible subsolvers are indispensable,giving rise to an infeasible inexact PALM(PALM-I).Numerous efforts have been devoted to analyzing the feasible PALM,while little attention has been paid to the PALM-I.The usage of the PALM-I thus lacks a theoretical guarantee.The essential difficulty of analysis consists in the objective value nonmonotonicity induced by the infeasibility.We study in the present work the convergence properties of the PALM-I.In particular,we construct a surrogate sequence to surmount the nonmonotonicity issue and devise an implementable inexact criterion.Based upon these,we manage to establish the stationarity of any accumulation point,and moreover,show the iterate convergence and the asymptotic convergence rates under the assumption of the Lojasiewicz property.The prominent advantages of the PALM-I on CPU time are illustrated via numerical experiments on problems arising from quantum physics and 3-dimensional anisotropic frictional contact.
