A generalization of the direct method of Clarkson and Kruskal for finding similarity reductions of partial differential equations with arbitrary functions is found and discussed for the generalized Burgers equation. T...A generalization of the direct method of Clarkson and Kruskal for finding similarity reductions of partial differential equations with arbitrary functions is found and discussed for the generalized Burgers equation. The corresponding reductions and the exact solutions due to the methods of the ordinary differential equations are then given by the methods. The results given here answer partially an open problem proposed by Clarkson, that is how to develop the direct method to seek symmetry reductions of nonlinear PDEs with arbitrary functions.展开更多
This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis techniq...This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.展开更多
Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks....Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.展开更多
Magnesium and its alloys are promising candidates for a new generation of biodegradable metals in orthopaedic applications due to their excellent biocompatibility,biodegradability,and mechanical properties that are si...Magnesium and its alloys are promising candidates for a new generation of biodegradable metals in orthopaedic applications due to their excellent biocompatibility,biodegradability,and mechanical properties that are similar to natural bone.However,direct in vitro assessment of these materials in the presence of cells is complicated by degradation products from the alloy that lead to a false positive for the most commonly used cell adhesion and cell proliferation assays.In this paper,a cyanine dye was used to quantitatively evaluate the in vitro biocompatibility of a Mg AZ31 alloy by both direct and indirect methods.The cytotoxicity of the corrosion products was evaluated via an indirect method;a 25%decrease in cell viability compared to control samples was observed.Moreover,direct assessment of cell adhesion and proliferation showed a statistically significant increase in cell number at the surface after 72 h.In addition,the degradation rate and surface characteristics of the Mg AZ31 alloy were evaluated for both direct and indirect tests.The degradation rate was unaffected by the presence of cells while evidence of an increase in calcium phosphate deposition on the magnesium alloy surface in the presence of cells was observed.This study demonstrates that a cyanine dye based assay provides a more accurate assessment of the overall in vitro biocompatibility of biodegradable metals than the more commonly used assays reported in the literature to date.展开更多
This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite differenc...This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.展开更多
A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysi...A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.展开更多
In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear conve...In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear convection-diffusion problems.In the spatial discretization,both the original DDG methods and the refined DDG methods with interface corrections are considered.In the time discretization,the convection term is treated explicitly and the diffusion term implicitly.By the energy method,we show that the corresponding fully discrete schemes are unconditionally stable,in the sense that the time-step is only required to be upper bounded by a constant which is independent of the mesh size h.Opti-mal error estimate is also obtained by the aid of a special global projection.Numerical experiments are given to verify the stability and accuracy of the proposed schemes.展开更多
<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the...<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the optimal error estimates. The numerical experiments imply that <em>L</em><sup>2 </sup>norms error estimates can reach to order <em>k</em> + 1 by using time discretization methods. </div>展开更多
The fixed-time synchronization and preassigned-time synchronization are investigated for a class of quaternion-valued neural networks with time-varying delays and discontinuous activation functions. Unlike previous ef...The fixed-time synchronization and preassigned-time synchronization are investigated for a class of quaternion-valued neural networks with time-varying delays and discontinuous activation functions. Unlike previous efforts that employed separation analysis and the real-valued control design, based on the quaternion-valued signum function and several related properties, a direct analytical method is proposed here and the quaternion-valued controllers are designed in order to discuss the fixed-time synchronization for the relevant quaternion-valued neural networks. In addition, the preassigned-time synchronization is investigated based on a quaternion-valued control design, where the synchronization time is preassigned and the control gains are finite. Compared with existing results, the direct method without separation developed in this article is beneficial in terms of simplifying theoretical analysis, and the proposed quaternion-valued control schemes are simpler and more effective than the traditional design, which adds four real-valued controllers. Finally, two numerical examples are given in order to support the theoretical results.展开更多
FL-Online(http://fanlab.ac.cn) is an out-of-box modern web service featuring a user-friendly interface and simplified parameters, providing academic users with access to a series of online programs for biomolecular cr...FL-Online(http://fanlab.ac.cn) is an out-of-box modern web service featuring a user-friendly interface and simplified parameters, providing academic users with access to a series of online programs for biomolecular crystallography, including SAPI-online, OASIS-online, C-IPCAS-online and a series of upcoming software releases. Meanwhile, it is a highly scalable and maintainable web application framework that provides a powerful and flexible solution for academic web development needs. All the codes are open-source under MIT licenses in GitHub.展开更多
To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direc...To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.展开更多
The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And...The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.展开更多
This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic So...This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.展开更多
A novel algorithm, i.e. the fast alternating direction method of multipliers （ADMM）, is applied to solve the classical total-variation （ TV ）-based model for image reconstruction. First, the TV-based model is refo...A novel algorithm, i.e. the fast alternating direction method of multipliers （ADMM）, is applied to solve the classical total-variation （ TV ）-based model for image reconstruction. First, the TV-based model is reformulated as a linear equality constrained problem where the objective function is separable. Then, by introducing the augmented Lagrangian function, the two variables are alternatively minimized by the Gauss-Seidel idea. Finally, the dual variable is updated. Because the approach makes full use of the special structure of the problem and decomposes the original problem into several low-dimensional sub-problems, the per iteration computational complexity of the approach is dominated by two fast Fourier transforms. Elementary experimental results indicate that the proposed approach is more stable and efficient compared with some state-of-the-art algorithms.展开更多
In practice,simultaneous impact localization and time history reconstruction can hardly be achieved,due to the illposed and under-determined problems induced by the constrained and harsh measuring conditions.Although ...In practice,simultaneous impact localization and time history reconstruction can hardly be achieved,due to the illposed and under-determined problems induced by the constrained and harsh measuring conditions.Although l_(1) regularization can be used to obtain sparse solutions,it tends to underestimate solution amplitudes as a biased estimator.To address this issue,a novel impact force identification method with l_(p) regularization is proposed in this paper,using the alternating direction method of multipliers(ADMM).By decomposing the complex primal problem into sub-problems solvable in parallel via proximal operators,ADMM can address the challenge effectively.To mitigate the sensitivity to regularization parameters,an adaptive regularization parameter is derived based on the K-sparsity strategy.Then,an ADMM-based sparse regularization method is developed,which is capable of handling l_(p) regularization with arbitrary p values using adaptively-updated parameters.The effectiveness and performance of the proposed method are validated on an aircraft skin-like composite structure.Additionally,an investigation into the optimal p value for achieving high-accuracy solutions via l_(p) regularization is conducted.It turns out that l_(0.6)regularization consistently yields sparser and more accurate solutions for impact force identification compared to the classic l_(1) regularization method.The impact force identification method proposed in this paper can simultaneously reconstruct impact time history with high accuracy and accurately localize the impact using an under-determined sensor configuration.展开更多
Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in te...Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.展开更多
The paper generalizes the direct method of moving planes to the Logarithmic Laplacian system.Firstly,some key ingredients of the method are discussed,for example,Narrow region principle and Decay at infinity.Then,the ...The paper generalizes the direct method of moving planes to the Logarithmic Laplacian system.Firstly,some key ingredients of the method are discussed,for example,Narrow region principle and Decay at infinity.Then,the radial symmetry of the solution of the Logarithmic Laplacian system is obtained.展开更多
In this paper,we design an efficient,multi-stage image segmentation framework that incorporates a weighted difference of anisotropic and isotropic total variation(AITV).The segmentation framework generally consists of...In this paper,we design an efficient,multi-stage image segmentation framework that incorporates a weighted difference of anisotropic and isotropic total variation(AITV).The segmentation framework generally consists of two stages:smoothing and thresholding,thus referred to as smoothing-and-thresholding(SaT).In the first stage,a smoothed image is obtained by an AITV-regularized Mumford-Shah(MS)model,which can be solved efficiently by the alternating direction method of multipliers(ADMMs)with a closed-form solution of a proximal operator of the l_(1)-αl_(2) regularizer.The convergence of the ADMM algorithm is analyzed.In the second stage,we threshold the smoothed image by K-means clustering to obtain the final segmentation result.Numerical experiments demonstrate that the proposed segmentation framework is versatile for both grayscale and color images,effcient in producing high-quality segmentation results within a few seconds,and robust to input images that are corrupted with noise,blur,or both.We compare the AITV method with its original convex TV and nonconvex TVP(O<p<1)counterparts,showcasing the qualitative and quantitative advantages of our proposed method.展开更多
Mechanical engineering structures and structural components are often subjected to cyclic thermomechanical loading which stresses their material beyond its elastic limits well inside the inelastic regime.Depending on ...Mechanical engineering structures and structural components are often subjected to cyclic thermomechanical loading which stresses their material beyond its elastic limits well inside the inelastic regime.Depending on the level of loading inelastic strains may lead either to failure,due to low cycle fatigue or ratcheting,or to safety,through elastic shakedown.Thus,it is important to estimate the asymptotic stress state of such structures.This state may be determined by cumbersome incremental time-stepping calculations.Direct methods,alternatively,have big computational advantages as they focus on the characteristics of these states and try to establish them,in a direct way,right from the beginning of the calculations.Among the very few such general-purpose direct methods,a powerful direct method which has been called RSDM has appeared in the literature.The method may directly predict any asymptotic state when the exact time history of the loading is known.The advantage of the method is due to the fact that it addresses the physics of the asymptotic cycle and exploits the cyclic nature of its expected residual stress distribution.Based on RSDM a method for the shakedown analysis of structures,called RSDM-S has also been developed.Despite most direct methods for shakedown,RSDM-S does not need an optimization algorithm for its implementation.Both RSDM and RSDM-S may be implemented in any Finite Element Code.A thorough review of both these methods,together with examples of implementation are presented herein.展开更多
One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respec...One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respectively at φh″＋｜△φh｜ and φh″-｜△φh｜. The probability for △φh being positive （P＋） can be derived based on the Cochran distribution in direct methods. Hence the SAD phase ambiguity can be resolved by multiplying the Gaussian function peaked at φh″＋｜△φh｜ with P＋ and multiplying the Gaussian function peaked at φh″-｜△φh｜ with P_ （=1- P＋）. The direct-method SAD h phasing has been proved powerful in breaking SAD phase ambiguities, in particular when anomalous-scattering signals are weak. However, the approximation of bimodal phase distributions by the sum of two Gaussian functions introduces considerable errors. In this paper we show that a much better approximation can be achieved by replacing the two Gaussian functions with two von Mises distributions. Test results showed that this leads to significant improvement on the efficiency of direct-method SAD-phasing.展开更多
基金the National Natural Science Foundation of China(1 990 1 0 2 7)
文摘A generalization of the direct method of Clarkson and Kruskal for finding similarity reductions of partial differential equations with arbitrary functions is found and discussed for the generalized Burgers equation. The corresponding reductions and the exact solutions due to the methods of the ordinary differential equations are then given by the methods. The results given here answer partially an open problem proposed by Clarkson, that is how to develop the direct method to seek symmetry reductions of nonlinear PDEs with arbitrary functions.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12071214)the Natural Science Foundation for Colleges and Universities of Jiangsu Province of China(Grant No.20KJB110011)+1 种基金supported by the National Science Foundation(Grant No.DMS-1620335)and the Simons Foundation(Grant No.637716)supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12272347).
文摘This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.
文摘Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.
基金We gratefully acknowledge the Ministry of Education in the Kingdom of Saudi Arabia for financial support.
文摘Magnesium and its alloys are promising candidates for a new generation of biodegradable metals in orthopaedic applications due to their excellent biocompatibility,biodegradability,and mechanical properties that are similar to natural bone.However,direct in vitro assessment of these materials in the presence of cells is complicated by degradation products from the alloy that lead to a false positive for the most commonly used cell adhesion and cell proliferation assays.In this paper,a cyanine dye was used to quantitatively evaluate the in vitro biocompatibility of a Mg AZ31 alloy by both direct and indirect methods.The cytotoxicity of the corrosion products was evaluated via an indirect method;a 25%decrease in cell viability compared to control samples was observed.Moreover,direct assessment of cell adhesion and proliferation showed a statistically significant increase in cell number at the surface after 72 h.In addition,the degradation rate and surface characteristics of the Mg AZ31 alloy were evaluated for both direct and indirect tests.The degradation rate was unaffected by the presence of cells while evidence of an increase in calcium phosphate deposition on the magnesium alloy surface in the presence of cells was observed.This study demonstrates that a cyanine dye based assay provides a more accurate assessment of the overall in vitro biocompatibility of biodegradable metals than the more commonly used assays reported in the literature to date.
文摘This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.
文摘A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.
基金the NSFC grant 11871428the Nature Science Research Program for Colleges and Universities of Jiangsu Province grant 20KJB110011Qiang Zhang:Research supported by the NSFC grant 11671199。
文摘In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear convection-diffusion problems.In the spatial discretization,both the original DDG methods and the refined DDG methods with interface corrections are considered.In the time discretization,the convection term is treated explicitly and the diffusion term implicitly.By the energy method,we show that the corresponding fully discrete schemes are unconditionally stable,in the sense that the time-step is only required to be upper bounded by a constant which is independent of the mesh size h.Opti-mal error estimate is also obtained by the aid of a special global projection.Numerical experiments are given to verify the stability and accuracy of the proposed schemes.
文摘<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the optimal error estimates. The numerical experiments imply that <em>L</em><sup>2 </sup>norms error estimates can reach to order <em>k</em> + 1 by using time discretization methods. </div>
基金supported by the National Natural Science Foundation of China (61963033, 61866036, 62163035)the Key Project of Natural Science Foundation of Xinjiang (2021D01D10)+1 种基金the Xinjiang Key Laboratory of Applied Mathematics (XJDX1401)the Special Project for Local Science and Technology Development Guided by the Central Government (ZYYD2022A05)。
文摘The fixed-time synchronization and preassigned-time synchronization are investigated for a class of quaternion-valued neural networks with time-varying delays and discontinuous activation functions. Unlike previous efforts that employed separation analysis and the real-valued control design, based on the quaternion-valued signum function and several related properties, a direct analytical method is proposed here and the quaternion-valued controllers are designed in order to discuss the fixed-time synchronization for the relevant quaternion-valued neural networks. In addition, the preassigned-time synchronization is investigated based on a quaternion-valued control design, where the synchronization time is preassigned and the control gains are finite. Compared with existing results, the direct method without separation developed in this article is beneficial in terms of simplifying theoretical analysis, and the proposed quaternion-valued control schemes are simpler and more effective than the traditional design, which adds four real-valued controllers. Finally, two numerical examples are given in order to support the theoretical results.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.32371280 and T2350011)。
文摘FL-Online(http://fanlab.ac.cn) is an out-of-box modern web service featuring a user-friendly interface and simplified parameters, providing academic users with access to a series of online programs for biomolecular crystallography, including SAPI-online, OASIS-online, C-IPCAS-online and a series of upcoming software releases. Meanwhile, it is a highly scalable and maintainable web application framework that provides a powerful and flexible solution for academic web development needs. All the codes are open-source under MIT licenses in GitHub.
基金financially supported by the National Natural Science Foundation of China(No.41574127 and 41174104)the National Key Technology R&D Program for the 13th five-year plan(No.2016ZX05018006-006)
文摘To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.
文摘The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.
基金Project supported by the Natural Science Foundation of China(10371009)Research Fund for the Doctoral Program Higher Education
文摘This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.
基金The Scientific Research Foundation of Nanjing University of Posts and Telecommunications(No.NY210049)
文摘A novel algorithm, i.e. the fast alternating direction method of multipliers （ADMM）, is applied to solve the classical total-variation （ TV ）-based model for image reconstruction. First, the TV-based model is reformulated as a linear equality constrained problem where the objective function is separable. Then, by introducing the augmented Lagrangian function, the two variables are alternatively minimized by the Gauss-Seidel idea. Finally, the dual variable is updated. Because the approach makes full use of the special structure of the problem and decomposes the original problem into several low-dimensional sub-problems, the per iteration computational complexity of the approach is dominated by two fast Fourier transforms. Elementary experimental results indicate that the proposed approach is more stable and efficient compared with some state-of-the-art algorithms.
基金Supported by National Natural Science Foundation of China (Grant Nos.52305127,52075414)China Postdoctoral Science Foundation (Grant No.2021M702595)。
文摘In practice,simultaneous impact localization and time history reconstruction can hardly be achieved,due to the illposed and under-determined problems induced by the constrained and harsh measuring conditions.Although l_(1) regularization can be used to obtain sparse solutions,it tends to underestimate solution amplitudes as a biased estimator.To address this issue,a novel impact force identification method with l_(p) regularization is proposed in this paper,using the alternating direction method of multipliers(ADMM).By decomposing the complex primal problem into sub-problems solvable in parallel via proximal operators,ADMM can address the challenge effectively.To mitigate the sensitivity to regularization parameters,an adaptive regularization parameter is derived based on the K-sparsity strategy.Then,an ADMM-based sparse regularization method is developed,which is capable of handling l_(p) regularization with arbitrary p values using adaptively-updated parameters.The effectiveness and performance of the proposed method are validated on an aircraft skin-like composite structure.Additionally,an investigation into the optimal p value for achieving high-accuracy solutions via l_(p) regularization is conducted.It turns out that l_(0.6)regularization consistently yields sparser and more accurate solutions for impact force identification compared to the classic l_(1) regularization method.The impact force identification method proposed in this paper can simultaneously reconstruct impact time history with high accuracy and accurately localize the impact using an under-determined sensor configuration.
基金supported by the NSFC Major Research Plan--Interpretable and Generalpurpose Next-generation Artificial Intelligence(No.92370205).
文摘Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.
基金Supported by the National Natural Science Foundation of China(11501342,12001344)。
文摘The paper generalizes the direct method of moving planes to the Logarithmic Laplacian system.Firstly,some key ingredients of the method are discussed,for example,Narrow region principle and Decay at infinity.Then,the radial symmetry of the solution of the Logarithmic Laplacian system is obtained.
基金partially supported by the NSF grants DMS-1854434,DMS-1952644,DMS-2151235,DMS-2219904,and CAREER 1846690。
文摘In this paper,we design an efficient,multi-stage image segmentation framework that incorporates a weighted difference of anisotropic and isotropic total variation(AITV).The segmentation framework generally consists of two stages:smoothing and thresholding,thus referred to as smoothing-and-thresholding(SaT).In the first stage,a smoothed image is obtained by an AITV-regularized Mumford-Shah(MS)model,which can be solved efficiently by the alternating direction method of multipliers(ADMMs)with a closed-form solution of a proximal operator of the l_(1)-αl_(2) regularizer.The convergence of the ADMM algorithm is analyzed.In the second stage,we threshold the smoothed image by K-means clustering to obtain the final segmentation result.Numerical experiments demonstrate that the proposed segmentation framework is versatile for both grayscale and color images,effcient in producing high-quality segmentation results within a few seconds,and robust to input images that are corrupted with noise,blur,or both.We compare the AITV method with its original convex TV and nonconvex TVP(O<p<1)counterparts,showcasing the qualitative and quantitative advantages of our proposed method.
文摘Mechanical engineering structures and structural components are often subjected to cyclic thermomechanical loading which stresses their material beyond its elastic limits well inside the inelastic regime.Depending on the level of loading inelastic strains may lead either to failure,due to low cycle fatigue or ratcheting,or to safety,through elastic shakedown.Thus,it is important to estimate the asymptotic stress state of such structures.This state may be determined by cumbersome incremental time-stepping calculations.Direct methods,alternatively,have big computational advantages as they focus on the characteristics of these states and try to establish them,in a direct way,right from the beginning of the calculations.Among the very few such general-purpose direct methods,a powerful direct method which has been called RSDM has appeared in the literature.The method may directly predict any asymptotic state when the exact time history of the loading is known.The advantage of the method is due to the fact that it addresses the physics of the asymptotic cycle and exploits the cyclic nature of its expected residual stress distribution.Based on RSDM a method for the shakedown analysis of structures,called RSDM-S has also been developed.Despite most direct methods for shakedown,RSDM-S does not need an optimization algorithm for its implementation.Both RSDM and RSDM-S may be implemented in any Finite Element Code.A thorough review of both these methods,together with examples of implementation are presented herein.
基金Project supported by the Innovation Foundation of the Chinese Academy of Sciences and by the National Basic Research Program of China(Grant No.2002CB713801)
文摘One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respectively at φh″＋｜△φh｜ and φh″-｜△φh｜. The probability for △φh being positive （P＋） can be derived based on the Cochran distribution in direct methods. Hence the SAD phase ambiguity can be resolved by multiplying the Gaussian function peaked at φh″＋｜△φh｜ with P＋ and multiplying the Gaussian function peaked at φh″-｜△φh｜ with P_ （=1- P＋）. The direct-method SAD h phasing has been proved powerful in breaking SAD phase ambiguities, in particular when anomalous-scattering signals are weak. However, the approximation of bimodal phase distributions by the sum of two Gaussian functions introduces considerable errors. In this paper we show that a much better approximation can be achieved by replacing the two Gaussian functions with two von Mises distributions. Test results showed that this leads to significant improvement on the efficiency of direct-method SAD-phasing.
