For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows....For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.展开更多
Compressed sensing(CS) has achieved great success in single noise removal. However, it cannot restore the images contaminated with mixed noise efficiently. This paper introduces nonlocal similarity and cosparsity insp...Compressed sensing(CS) has achieved great success in single noise removal. However, it cannot restore the images contaminated with mixed noise efficiently. This paper introduces nonlocal similarity and cosparsity inspired by compressed sensing to overcome the difficulties in mixed noise removal, in which nonlocal similarity explores the signal sparsity from similar patches, and cosparsity assumes that the signal is sparse after a possibly redundant transform. Meanwhile, an adaptive scheme is designed to keep the balance between mixed noise removal and detail preservation based on local variance. Finally, IRLSM and RACoSaMP are adopted to solve the objective function. Experimental results demonstrate that the proposed method is superior to conventional CS methods, like K-SVD and state-of-art method nonlocally centralized sparse representation(NCSR), in terms of both visual results and quantitative measures.展开更多
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.展开更多
Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in prac...Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
A new method to reconstruct 3D scene points from nonparallel stereo is proposed. From a pair of conjugate images in an arbitrarily configured stereo system that has been calibrated, coordinates of 3D scene points can ...A new method to reconstruct 3D scene points from nonparallel stereo is proposed. From a pair of conjugate images in an arbitrarily configured stereo system that has been calibrated, coordinates of 3D scene points can be computed directly using the method, bypassing the process of rectifying images or iterative solution involved in existing methods. Experiment results from both simulated data and real images validate the method. Practical application to surgical navigator shows that the method has advantages to improve efficiency and accuracy of 3D reconstruction from nonparallel stereo system in comparison with the conventional method that employs algorithm for standard parallel axes stereo geometry.展开更多
AIM To evaluate the results of the Ilizarov method in management of cases with bone loss, soft tissue loss and infection.METHODS Twenty eight patients with severe leg trauma complicated by bone loss, soft tissue loss ...AIM To evaluate the results of the Ilizarov method in management of cases with bone loss, soft tissue loss and infection.METHODS Twenty eight patients with severe leg trauma complicated by bone loss, soft tissue loss and infection were managed by distraction osteogenesis in our institution. After radical debridement of all the infected and dead tissues the Ilizarov frame was applied, corticotomy was done and bone transport started. The wounds were left open to drain. Partial limb shortening was done in seven cases to reduce the size of both the skeletal and soft tissue defects. The average follow up period was 39 mo(range 27-56 mo). RESULTS The infection was eradicated in all cases. All the soft tissue defects healed during bone transport and plastic surgery was only required in 2 cases. Skeletal defects were treated in all cases. All patients required another surgery at the docking site to fashion the soft tissue and to cover the bone ends. The external fixation time ranged from 9 to 17 mo with an average of 13 mo. The complications included pin tract infection in 16 cases, wire breakage in 2 cases, unstable scar in 4 cases and chronic edema in 3 cases. According to the association for study and application of methods of Ilizarov score the bone results were excellent in 10, good in 16 and fair in 2 cases while the functional results were excellent in 8, good in 17 and fair in 3 cases. CONCLUSION Distraction osteogenesis is a good method that can treat the three problems of this triad simultaneously.展开更多
The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectra...The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.展开更多
Online reactivity monitoring plays an important role in operation and safety analyses of fission reactor systems. The inverse kinetics method, which is based on a point kinetics model, is the most widely used method f...Online reactivity monitoring plays an important role in operation and safety analyses of fission reactor systems. The inverse kinetics method, which is based on a point kinetics model, is the most widely used method for reactivity reconstruction of critical water reactors. However, this method is seldom applied to the reactivity reconstruction of subcritical reactors. In this study, an inverse kinetics method was employed for the reactivity reconstruction of a lead-based reactor under different initial reactivity states(ρ_0= 0,-2786,-5486,-8367, and-12,371 pcm). The results showed that the deviation in the reactivity of the lead-based subcritical reactor was greater when ρ_0 became smaller. The reactivity reconstructed using the inverse kinetics method was globally underestimated. At a given reactivity perturbation, the relative and absolute errors increased with the decrease in the initial reactivity. At a given initial reactivity, with the increase in the reactivity perturbation, the absolute error increased, whereas the relative error remained the same.This deviation is due to the variation in the external neutron source, spatial-spectral effects, and sub-diffusive effects, which require further study.展开更多
A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluati...A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation,their characters,advantages and shortcomings are compared,analyzed and commended in particular.展开更多
This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior p...This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest （ROI） accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection （FBP）.In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.展开更多
基金supported by the National Natural Science Foundation of China(Nos.11702329,12102247)the Shanghai Key Lab of Vehicle Aerodynamics and Vehicle Thermal Management Systems,China(No.VATLAB-2021-01)。
文摘For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.
基金supported by the National Natural Science Foundation of China(Nos.61403146 and 61603105)the Fundamental Research Funds for the Central Universities(No.2015ZM128)the Science and Technology Program of Guangzhou in China(Nos.201707010054 and 201704030072)
文摘Compressed sensing(CS) has achieved great success in single noise removal. However, it cannot restore the images contaminated with mixed noise efficiently. This paper introduces nonlocal similarity and cosparsity inspired by compressed sensing to overcome the difficulties in mixed noise removal, in which nonlocal similarity explores the signal sparsity from similar patches, and cosparsity assumes that the signal is sparse after a possibly redundant transform. Meanwhile, an adaptive scheme is designed to keep the balance between mixed noise removal and detail preservation based on local variance. Finally, IRLSM and RACoSaMP are adopted to solve the objective function. Experimental results demonstrate that the proposed method is superior to conventional CS methods, like K-SVD and state-of-art method nonlocally centralized sparse representation(NCSR), in terms of both visual results and quantitative measures.
基金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.
基金the National High Technology Research and Development Program of China(Grant No.2012AA011603)
文摘Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
基金The National Natural Science Foundation of China(No60675017)
文摘A new method to reconstruct 3D scene points from nonparallel stereo is proposed. From a pair of conjugate images in an arbitrarily configured stereo system that has been calibrated, coordinates of 3D scene points can be computed directly using the method, bypassing the process of rectifying images or iterative solution involved in existing methods. Experiment results from both simulated data and real images validate the method. Practical application to surgical navigator shows that the method has advantages to improve efficiency and accuracy of 3D reconstruction from nonparallel stereo system in comparison with the conventional method that employs algorithm for standard parallel axes stereo geometry.
文摘AIM To evaluate the results of the Ilizarov method in management of cases with bone loss, soft tissue loss and infection.METHODS Twenty eight patients with severe leg trauma complicated by bone loss, soft tissue loss and infection were managed by distraction osteogenesis in our institution. After radical debridement of all the infected and dead tissues the Ilizarov frame was applied, corticotomy was done and bone transport started. The wounds were left open to drain. Partial limb shortening was done in seven cases to reduce the size of both the skeletal and soft tissue defects. The average follow up period was 39 mo(range 27-56 mo). RESULTS The infection was eradicated in all cases. All the soft tissue defects healed during bone transport and plastic surgery was only required in 2 cases. Skeletal defects were treated in all cases. All patients required another surgery at the docking site to fashion the soft tissue and to cover the bone ends. The external fixation time ranged from 9 to 17 mo with an average of 13 mo. The complications included pin tract infection in 16 cases, wire breakage in 2 cases, unstable scar in 4 cases and chronic edema in 3 cases. According to the association for study and application of methods of Ilizarov score the bone results were excellent in 10, good in 16 and fair in 2 cases while the functional results were excellent in 8, good in 17 and fair in 3 cases. CONCLUSION Distraction osteogenesis is a good method that can treat the three problems of this triad simultaneously.
基金Research was supported in part by NSF grant DMS-0800612Research was supported by Applied Mathematics program of the US DOE Office of Advanced Scientific Computing ResearchThe Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830
文摘The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.
基金supported by the Strategic Priority Science and Technology Program of the Chinese Academy of Sciences(No.XDA03040000)the National Natural Science Foundation of China(NSFC)(Nos.11305205,11305203,and 11405204)+3 种基金the Special Program for Informatization of the Chinese Academy of Sciences(No.XXH12504-1-09)the Anhui Provincial Special project for High Technology Industrythe Special Project of Youth Innovation Promotion Association of Chinese Academy of Sciencesthe Industrialization Fund
文摘Online reactivity monitoring plays an important role in operation and safety analyses of fission reactor systems. The inverse kinetics method, which is based on a point kinetics model, is the most widely used method for reactivity reconstruction of critical water reactors. However, this method is seldom applied to the reactivity reconstruction of subcritical reactors. In this study, an inverse kinetics method was employed for the reactivity reconstruction of a lead-based reactor under different initial reactivity states(ρ_0= 0,-2786,-5486,-8367, and-12,371 pcm). The results showed that the deviation in the reactivity of the lead-based subcritical reactor was greater when ρ_0 became smaller. The reactivity reconstructed using the inverse kinetics method was globally underestimated. At a given reactivity perturbation, the relative and absolute errors increased with the decrease in the initial reactivity. At a given initial reactivity, with the increase in the reactivity perturbation, the absolute error increased, whereas the relative error remained the same.This deviation is due to the variation in the external neutron source, spatial-spectral effects, and sub-diffusive effects, which require further study.
文摘A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation,their characters,advantages and shortcomings are compared,analyzed and commended in particular.
基金supported by the National Natural Science Foundation of China (Grant No.60872116)
文摘This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest （ROI） accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection （FBP）.In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.
