The accuracy of gradient reconstruction methods on unstructured meshes is analyzed both mathematically and numerically.Mathematical derivations reveal that,for gradient reconstruction based on the Green-Gauss theorem(...The accuracy of gradient reconstruction methods on unstructured meshes is analyzed both mathematically and numerically.Mathematical derivations reveal that,for gradient reconstruction based on the Green-Gauss theorem(the GG methods),if the summation of first-and-lower-order terms does not counterbalance in the discretized integral process,which rarely occurs,second-order accurate approximation of face midpoint value is necessary to produce at least first-order accurate gradient.However,gradient reconstruction based on the least-squares approach(the LSQ methods)is at least first-order on arbitrary unstructured grids.Verifications are performed on typical isotropic grid stencils by analyzing the relationship between the discretization error of gradient reconstruction and the discretization error of the face midpoint value approximation of a given analytic function.Meanwhile,the numerical accuracy of gradient reconstruction methods is examined with grid convergence study on typical isotropic grids.Results verify the phenomenon of accuracy degradation for the GG methods when the face midpoint value condition is not satisfied.The LSQ methods are proved to be at least first-order on all tested isotropic grids.To study gradient accuracy effects on inviscid flow simulation,solution errors are quantified using the Method of Manufactured Solutions(MMS)which was validated before adoption by comparing with an exact solution case,i.e.,the 2-dimensional(2D)inviscid isentropic vortex.Numerical results demonstrate that the order of accuracy(OOA)of gradient reconstruction is crucial in determining the OOA of numerical solutions.Solution accuracy deteriorates seriously if gradient reconstruction does not reach first-order.展开更多
The additional sparse prior of images has been the subject of much research in problems of sparse-view computed tomography(CT) reconstruction. A method employing the image gradient sparsity is often used to reduce t...The additional sparse prior of images has been the subject of much research in problems of sparse-view computed tomography(CT) reconstruction. A method employing the image gradient sparsity is often used to reduce the sampling rate and is shown to remove the unwanted artifacts while preserve sharp edges, but may cause blocky or patchy artifacts.To eliminate this drawback, we propose a novel sparsity exploitation-based model for CT image reconstruction. In the presented model, the sparse representation and sparsity exploitation of both gradient and nonlocal gradient are investigated.The new model is shown to offer the potential for better results by introducing a similarity prior information of the image structure. Then, an effective alternating direction minimization algorithm is developed to optimize the objective function with a robust convergence result. Qualitative and quantitative evaluations have been carried out both on the simulation and real data in terms of accuracy and resolution properties. The results indicate that the proposed method can be applied for achieving better image-quality potential with the theoretically expected detailed feature preservation.展开更多
This paper is focused on studying an important concept of the system analysis, which is the regional enlarged observability or constrained observability of the gradient for distributed parabolic systems evolving in th...This paper is focused on studying an important concept of the system analysis, which is the regional enlarged observability or constrained observability of the gradient for distributed parabolic systems evolving in the spatial domain Ω We will explore an approach based on the Hilbert Uniqueness Method (HUM), which can reconstruct the initial gradient state between two prescribed functions f1 and f2 only in a critical subregion ω of Ω without the knowledge of the state. Finally, the obtained results are illustrated by numerical simulations.展开更多
The aim of this work is to study the notion of the gradient observability on a subregion?ω of the evolution domain?Ω for a class of semilinear hyperbolic systems. We show, under some hypothesis, that the gradient re...The aim of this work is to study the notion of the gradient observability on a subregion?ω of the evolution domain?Ω for a class of semilinear hyperbolic systems. We show, under some hypothesis, that the gradient reconstruction is achieved following sectorial approach combined with fixed point techniques. The obtained results lead to an algorithm which can be implemented numerically.展开更多
This paper presents a method to deal with an extension of regional gradient observability developed for parabolic system [1,2] to hyperbolic one. This concerns the reconstruction of the state gradient only on a subreg...This paper presents a method to deal with an extension of regional gradient observability developed for parabolic system [1,2] to hyperbolic one. This concerns the reconstruction of the state gradient only on a subregion of the system domain. Then necessary conditions for sensors structure are established in order to obtain regional gradient observability. An approach is developed which allows the reconstruction of the system state gradient on a given subregion. The obtained results are illustrated by numerical examples and simulations.展开更多
A fast converging sparse reconstruction algorithm in ghost imaging is presented. It utilizes total variation regularization and its formulation is based on the Karush-Kuhn-Tucker (KKT) theorem in the theory of convex ...A fast converging sparse reconstruction algorithm in ghost imaging is presented. It utilizes total variation regularization and its formulation is based on the Karush-Kuhn-Tucker (KKT) theorem in the theory of convex optimization. Tests using experimental data show that, compared with the algorithm of Gradient Projection for Sparse Reconstruction (GPSR), the proposed algorithm yields better results with less computation work.展开更多
The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode disc...The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode discharge with different auxiliary heatings are reconstructed by using magnetic diagnostics,kinetic diagnostics and TRANSP code.It is found that the fast-ion pressure might be up to one-third of the plasma pressure and the contribution is mainly in the core plasma due to the neutral beam injection power is primarily deposited in the core region.The fast-ion current contributes mainly in the core region while contributes little to the pedestal current.A steep pressure gradient in the pedestal is observed which gives rise to a strong edge current.It is proved that the fast ion effects cannot be ignored and should be considered in the future study of EAST.展开更多
The aim of this paper is to study the notion of the gradient observability on a subregion w of the evolution domain W and also we consider the case where the subregion of interest is a boundary part of the system evol...The aim of this paper is to study the notion of the gradient observability on a subregion w of the evolution domain W and also we consider the case where the subregion of interest is a boundary part of the system evolution domain for the class of semilinear hyperbolic systems. We show, under some hypotheses, that the flux reconstruction is guaranteed by means of the sectorial approach combined with fixed point techniques. This leads to several interesting results which are performed through numerical examples and simulations.展开更多
基金National Natural Science Foundation of China[grant numbers 11532016,91530325].
文摘The accuracy of gradient reconstruction methods on unstructured meshes is analyzed both mathematically and numerically.Mathematical derivations reveal that,for gradient reconstruction based on the Green-Gauss theorem(the GG methods),if the summation of first-and-lower-order terms does not counterbalance in the discretized integral process,which rarely occurs,second-order accurate approximation of face midpoint value is necessary to produce at least first-order accurate gradient.However,gradient reconstruction based on the least-squares approach(the LSQ methods)is at least first-order on arbitrary unstructured grids.Verifications are performed on typical isotropic grid stencils by analyzing the relationship between the discretization error of gradient reconstruction and the discretization error of the face midpoint value approximation of a given analytic function.Meanwhile,the numerical accuracy of gradient reconstruction methods is examined with grid convergence study on typical isotropic grids.Results verify the phenomenon of accuracy degradation for the GG methods when the face midpoint value condition is not satisfied.The LSQ methods are proved to be at least first-order on all tested isotropic grids.To study gradient accuracy effects on inviscid flow simulation,solution errors are quantified using the Method of Manufactured Solutions(MMS)which was validated before adoption by comparing with an exact solution case,i.e.,the 2-dimensional(2D)inviscid isentropic vortex.Numerical results demonstrate that the order of accuracy(OOA)of gradient reconstruction is crucial in determining the OOA of numerical solutions.Solution accuracy deteriorates seriously if gradient reconstruction does not reach first-order.
基金Manuscript received February 13, 2016 accepted December 7, 2016. This work was supported by the National Natural Science Foundation of China (61362001, 61661031), Jiangxi Province Innovation Projects for Postgraduate Funds (YC2016-S006), the International Postdoctoral Exchange Fellowship Program, and Jiangxi Advanced Project for Post-Doctoral Research Fund (2014KY02).
基金Project supported by the National Natural Science Foundation of China(Grant No.61372172)
文摘The additional sparse prior of images has been the subject of much research in problems of sparse-view computed tomography(CT) reconstruction. A method employing the image gradient sparsity is often used to reduce the sampling rate and is shown to remove the unwanted artifacts while preserve sharp edges, but may cause blocky or patchy artifacts.To eliminate this drawback, we propose a novel sparsity exploitation-based model for CT image reconstruction. In the presented model, the sparse representation and sparsity exploitation of both gradient and nonlocal gradient are investigated.The new model is shown to offer the potential for better results by introducing a similarity prior information of the image structure. Then, an effective alternating direction minimization algorithm is developed to optimize the objective function with a robust convergence result. Qualitative and quantitative evaluations have been carried out both on the simulation and real data in terms of accuracy and resolution properties. The results indicate that the proposed method can be applied for achieving better image-quality potential with the theoretically expected detailed feature preservation.
文摘This paper is focused on studying an important concept of the system analysis, which is the regional enlarged observability or constrained observability of the gradient for distributed parabolic systems evolving in the spatial domain Ω We will explore an approach based on the Hilbert Uniqueness Method (HUM), which can reconstruct the initial gradient state between two prescribed functions f1 and f2 only in a critical subregion ω of Ω without the knowledge of the state. Finally, the obtained results are illustrated by numerical simulations.
文摘The aim of this work is to study the notion of the gradient observability on a subregion?ω of the evolution domain?Ω for a class of semilinear hyperbolic systems. We show, under some hypothesis, that the gradient reconstruction is achieved following sectorial approach combined with fixed point techniques. The obtained results lead to an algorithm which can be implemented numerically.
文摘This paper presents a method to deal with an extension of regional gradient observability developed for parabolic system [1,2] to hyperbolic one. This concerns the reconstruction of the state gradient only on a subregion of the system domain. Then necessary conditions for sensors structure are established in order to obtain regional gradient observability. An approach is developed which allows the reconstruction of the system state gradient on a given subregion. The obtained results are illustrated by numerical examples and simulations.
基金Supported by the Hi-Tech Research and Development Program of China (No. 2011AA120102)
文摘A fast converging sparse reconstruction algorithm in ghost imaging is presented. It utilizes total variation regularization and its formulation is based on the Karush-Kuhn-Tucker (KKT) theorem in the theory of convex optimization. Tests using experimental data show that, compared with the algorithm of Gradient Projection for Sparse Reconstruction (GPSR), the proposed algorithm yields better results with less computation work.
基金supported by National Key R&D Program of China under Grant No.2017YFE0300400National Natural Science Foundation of China under Grant Nos.11475220,11405218,11575248+1 种基金the National Magnetic Confinement Fusion Science Program of China under Contracts No.2014GB106001sponsored in part by Youth Innovation Promotion Association Chinese Academy of Sciences (Grant No.2016384)
文摘The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode discharge with different auxiliary heatings are reconstructed by using magnetic diagnostics,kinetic diagnostics and TRANSP code.It is found that the fast-ion pressure might be up to one-third of the plasma pressure and the contribution is mainly in the core plasma due to the neutral beam injection power is primarily deposited in the core region.The fast-ion current contributes mainly in the core region while contributes little to the pedestal current.A steep pressure gradient in the pedestal is observed which gives rise to a strong edge current.It is proved that the fast ion effects cannot be ignored and should be considered in the future study of EAST.
文摘The aim of this paper is to study the notion of the gradient observability on a subregion w of the evolution domain W and also we consider the case where the subregion of interest is a boundary part of the system evolution domain for the class of semilinear hyperbolic systems. We show, under some hypotheses, that the flux reconstruction is guaranteed by means of the sectorial approach combined with fixed point techniques. This leads to several interesting results which are performed through numerical examples and simulations.