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.展开更多
Using the (2+1)-dimensional Broer-Kaup equation as an simple example, a new direct method is developed to find symmetry groups and symmetry algebras and then exact solutions of nonlinear mathematical physical equations.
A novel control method for a general class of nonlinear systems using fuzzy logic systems (FLSs) is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a ...A novel control method for a general class of nonlinear systems using fuzzy logic systems (FLSs) is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a high-gain observer is used to estimate the derivatives of the system output. The closed-loop system is proven to be semiglobally uniformly ultimately bounded. In addition, it is shown that if the approximation accuracy of the fuzzy logic system is high enough and the observer gain is chosen sufficiently large, an arbitrarily small tracking error can be achieved. Simulation results verify the effectiveness of the newly designed scheme and the theoretical discussion.展开更多
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 direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions an...This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.展开更多
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.展开更多
The present study aims at the design and making of measuring instrument of whole direct method for bed shear stress under two-dimensional water-flow co-action. The instrument combines the traditional strain gauge with...The present study aims at the design and making of measuring instrument of whole direct method for bed shear stress under two-dimensional water-flow co-action. The instrument combines the traditional strain gauge with a precise pressure gauge, and adopts the method directly measuring the difference between the lateral hydrodynamic pressure and different head pressures on both sides of the force plate. As a result, such an instrument solves a technical puzzle of the past strain gauge, i.e. the difficulty to set apart shear stress and lateral force. Static force test and sink test both prove that the instrument is precise, stable and applicable to the measurement of rough beds with different shear stresses.展开更多
In this article, by introducing characteristic singular integral operator and associate singular integral equations (SIEs), the authors discuss the direct method of solution for a class of singular integral equation...In this article, by introducing characteristic singular integral operator and associate singular integral equations (SIEs), the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.展开更多
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.展开更多
In this paper,based on the symbolic computing system Maple,the direct method for Lie symmetry groupspresented by Sen-Yue Lou [J.Phys.A:Math.Gen.38 (2005) L129] is extended from the continuous differential equationsto ...In this paper,based on the symbolic computing system Maple,the direct method for Lie symmetry groupspresented by Sen-Yue Lou [J.Phys.A:Math.Gen.38 (2005) L129] is extended from the continuous differential equationsto the differential-difference equations.With the extended method,we study the well-known differential-difference KPequation,KZ equation and (2+1)-dimensional ANNV system,and both the Lie point symmetry groups and the non-Liesymmetry groups are obtained.展开更多
A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, th...A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.展开更多
In this paper, first, we employ classic Lie symmetry groups approach to obtain the Lie symmetry groupsof the well-known (2+1)-dimensional Generalized Sasa-Satsuma (GSS) equation. Second, based on a modified directmeth...In this paper, first, we employ classic Lie symmetry groups approach to obtain the Lie symmetry groupsof the well-known (2+1)-dimensional Generalized Sasa-Satsuma (GSS) equation. Second, based on a modified directmethod proposed by Lou [J. Phys. A: Math. Gen. 38 (2005) L129], more general symmetry groups are obtained andthe relationship between the new solution and known solution is set up. At the same time, the Lie symmetry groupsobtained are only special cases of the more general symmetry groups. At last, some exact solutions of GSS equationsare constructed by the relationship obtained in the paper between the new solution and known solution.展开更多
The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolat...The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.展开更多
The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor ...The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor both methods, series reduction solutions are consequently derived.Higher order similarity reduction equations arelinear variable coefficients ordinary differential equations.By comparison, it is find that the results generated from theapproximate direct method are more general than the results generated from the approximate symmetry perturbationmethod.展开更多
This paper attempts to further extend the so-called Trefftz direct method(TDM), which has been developed and applied to a wide variety of boundary value problems in recent years. In this approach, the complete system ...This paper attempts to further extend the so-called Trefftz direct method(TDM), which has been developed and applied to a wide variety of boundary value problems in recent years. In this approach, the complete system of solutions of the partial differential equations is used as weighting functions and a non-singular boundary integral equation is used as the starting formulation. Several relative problems are discussed here. They are the problems of necessary and sufficient conditions for Trefftz-type boundary integral equations: the relationship between the Trefftz method and the variational principle and the infrequently encountered but possible singular problems in the Trefftz method and the treatment in practice. Numerical results are presented to illustrate the computational efficiency of the approach and to demonstrate its advantages.展开更多
Forest canopy water storage (S), direct throughfall fraction (p) and mean evaporation rate to mean rainfall intensity ratio (E/R) vary between storms and seasonally. Typically, researchers only quantify the mean growi...Forest canopy water storage (S), direct throughfall fraction (p) and mean evaporation rate to mean rainfall intensity ratio (E/R) vary between storms and seasonally. Typically, researchers only quantify the mean growing and dormant season values of S, p and E/R for deciduous forests, thereby ignoring seasonal changes S, p and E/R .Past researchers adapted the mean method, which is usually used to estimate S, p and E/R on an annual or seasonal basis, to estimate the same canopy variables on a per storm basis (individual storm (IS) method). The disadvantage of the IS method is that it requires more expensive equipment and the calculation of the canopy variables is more labor intensive relative to the mean method. The goal of this study was to explore the use of the IS method for northern hardwood forests and to determine whether estimates of S, p and E/R derived by the IS method produce more accurate estimates of rainfall interception loss (In) using the Gash model relative to estimates derived by the mean method. The IS method estimated that S increased from approximately 0.11 mm in the early spring to 1.2 mm in the summer and then declined to 0.24 mm after fall senescence. Direct throughfall decreased from 0.4 in the early spring to 0.11 in the summer, and then increased to 0.4 after leaf senescence. When measurement period estimates of p, S and E/R derived by the IS and mean methods were applied to the Gash model, the modeled estimates of In differed from the measured values by 14.0 mm and 1.3 mm, respectively. Therefore, because the mean method provided more accurate estimates of In, the extra effort and expense required by the IS method is not advantageous for studies in northern hardwood forests that only need to model annual or seasonal estimates of In.展开更多
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.展开更多
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.展开更多
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.展开更多
文摘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.
文摘Using the (2+1)-dimensional Broer-Kaup equation as an simple example, a new direct method is developed to find symmetry groups and symmetry algebras and then exact solutions of nonlinear mathematical physical equations.
基金This project was supported by the National Natural Science Foundation of China (90405011).
文摘A novel control method for a general class of nonlinear systems using fuzzy logic systems (FLSs) is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a high-gain observer is used to estimate the derivatives of the system output. The closed-loop system is proven to be semiglobally uniformly ultimately bounded. In addition, it is shown that if the approximation accuracy of the fuzzy logic system is high enough and the observer gain is chosen sufficiently large, an arbitrarily small tracking error can be achieved. Simulation results verify the effectiveness of the newly designed scheme and the theoretical discussion.
基金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.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
基金Project supported by the National Natural Science Foundations of China(Grant Nos.10735030,10475055,10675065 and 90503006)the National Basic Research Program of China(Grant No.2007CB814800)
文摘This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.
文摘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.
基金supported by the National Major Special Project for Scientific Instruments and Equipment Development“Intelligent Measurement and Control System Development for Large-Scale River Model Tests in China”(Grant No.2011YQ070055)
文摘The present study aims at the design and making of measuring instrument of whole direct method for bed shear stress under two-dimensional water-flow co-action. The instrument combines the traditional strain gauge with a precise pressure gauge, and adopts the method directly measuring the difference between the lateral hydrodynamic pressure and different head pressures on both sides of the force plate. As a result, such an instrument solves a technical puzzle of the past strain gauge, i.e. the difficulty to set apart shear stress and lateral force. Static force test and sink test both prove that the instrument is precise, stable and applicable to the measurement of rough beds with different shear stresses.
基金Project was supported by RFDP of Higher Education and NNSF of China, SF of Wuhan University
文摘In this article, by introducing characteristic singular integral operator and associate singular integral equations (SIEs), the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.
基金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.
文摘In this paper,based on the symbolic computing system Maple,the direct method for Lie symmetry groupspresented by Sen-Yue Lou [J.Phys.A:Math.Gen.38 (2005) L129] is extended from the continuous differential equationsto the differential-difference equations.With the extended method,we study the well-known differential-difference KPequation,KZ equation and (2+1)-dimensional ANNV system,and both the Lie point symmetry groups and the non-Liesymmetry groups are obtained.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10672044)
文摘A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.
基金Supported by the National Natural Science Foundation of China under Grant No. 10735030Shanghai Leading Academic Discipline Project under Grant No. B412+2 种基金National Natural Science Foundation of China under Grant No. 90718041Program for Changjiang Scholars and Innovative Research Team in University under Grant No. IRT0734K.C. Wong Magna Fund in Ningbo University
文摘In this paper, first, we employ classic Lie symmetry groups approach to obtain the Lie symmetry groupsof the well-known (2+1)-dimensional Generalized Sasa-Satsuma (GSS) equation. Second, based on a modified directmethod proposed by Lou [J. Phys. A: Math. Gen. 38 (2005) L129], more general symmetry groups are obtained andthe relationship between the new solution and known solution is set up. At the same time, the Lie symmetry groupsobtained are only special cases of the more general symmetry groups. At last, some exact solutions of GSS equationsare constructed by the relationship obtained in the paper between the new solution and known solution.
文摘The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.
基金Supported by the National Natural Science Foundation of China under Grant Nos.10735030,10475055,10675065,and 90503006National Basic Research Program of China (973 Program 2007CB814800)
文摘The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor both methods, series reduction solutions are consequently derived.Higher order similarity reduction equations arelinear variable coefficients ordinary differential equations.By comparison, it is find that the results generated from theapproximate direct method are more general than the results generated from the approximate symmetry perturbationmethod.
基金Supported by National Natural Science Foundation of China(No.19872019).
文摘This paper attempts to further extend the so-called Trefftz direct method(TDM), which has been developed and applied to a wide variety of boundary value problems in recent years. In this approach, the complete system of solutions of the partial differential equations is used as weighting functions and a non-singular boundary integral equation is used as the starting formulation. Several relative problems are discussed here. They are the problems of necessary and sufficient conditions for Trefftz-type boundary integral equations: the relationship between the Trefftz method and the variational principle and the infrequently encountered but possible singular problems in the Trefftz method and the treatment in practice. Numerical results are presented to illustrate the computational efficiency of the approach and to demonstrate its advantages.
文摘Forest canopy water storage (S), direct throughfall fraction (p) and mean evaporation rate to mean rainfall intensity ratio (E/R) vary between storms and seasonally. Typically, researchers only quantify the mean growing and dormant season values of S, p and E/R for deciduous forests, thereby ignoring seasonal changes S, p and E/R .Past researchers adapted the mean method, which is usually used to estimate S, p and E/R on an annual or seasonal basis, to estimate the same canopy variables on a per storm basis (individual storm (IS) method). The disadvantage of the IS method is that it requires more expensive equipment and the calculation of the canopy variables is more labor intensive relative to the mean method. The goal of this study was to explore the use of the IS method for northern hardwood forests and to determine whether estimates of S, p and E/R derived by the IS method produce more accurate estimates of rainfall interception loss (In) using the Gash model relative to estimates derived by the mean method. The IS method estimated that S increased from approximately 0.11 mm in the early spring to 1.2 mm in the summer and then declined to 0.24 mm after fall senescence. Direct throughfall decreased from 0.4 in the early spring to 0.11 in the summer, and then increased to 0.4 after leaf senescence. When measurement period estimates of p, S and E/R derived by the IS and mean methods were applied to the Gash model, the modeled estimates of In differed from the measured values by 14.0 mm and 1.3 mm, respectively. Therefore, because the mean method provided more accurate estimates of In, the extra effort and expense required by the IS method is not advantageous for studies in northern hardwood forests that only need to model annual or seasonal estimates of In.
基金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.
基金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 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.