In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patte...Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.展开更多
The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the p...The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.展开更多
Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly effi...Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly efficient calculations for O-CNOPs are still challenging in the field of ensemble forecasting.In this study,we combine a gradient-based iterative idea with the Gram‒Schmidt orthogonalization,and propose an iterative optimization method to compute O-CNOPs.This method is different from the original sequential optimization method,and allows parallel computations of O-CNOPs,thus saving a large amount of computational time.We evaluate this method by using the Lorenz-96 model on the basis of the ensemble forecasting ability achieved and on the time consumed for computing O-CNOPs.The results demonstrate that the parallel iterative method causes O-CNOPs to yield reliable ensemble members and to achieve ensemble forecasting skills similar to or even slightly higher than those produced by the sequential method.Moreover,the parallel method significantly reduces the computational time for O-CNOPs.Therefore,the parallel iterative method provides a highly effective and efficient approach for calculating O-CNOPs for ensemble forecasts.Expectedly,it can play an important role in the application of the O-CNOPs to realistic ensemble forecasts for high-impact weather and climate events.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
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.展开更多
Lithium-ion batteries(LIBs)have evolved into the mainstream power source of ene rgy sto rage equipment by reason of their advantages such as high energy density,high power,long cycle life and less pollution.With the e...Lithium-ion batteries(LIBs)have evolved into the mainstream power source of ene rgy sto rage equipment by reason of their advantages such as high energy density,high power,long cycle life and less pollution.With the expansion of their applications in deep-sea exploration,aerospace and military equipment,special working conditions have placed higher demands on the low-temperature performance of LIBs.However,at low temperatures,the severe polarization and inferior electrochemical activity of electrode materials cause the acute capacity fading upon cycling,which greatly hindered the further development of LIBs.In this review,we summarize the recent important progress of LIBs in low-temperature operations and introduce the key methods and the related action mechanisms for enhancing the capacity of the various cathode and anode materials.It aims to promote the development of high-performance electrode materials and broaden the application range of LIBs.展开更多
This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are glob...This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.展开更多
Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary a...Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms （EAs） and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.展开更多
In order to thoroughly investigate the diversity of glacier microorganisms, four DNA extraction methods with differem lysis pat- terns were tested and two screened methods （the Bosshard-Bano method and the Zhou metho...In order to thoroughly investigate the diversity of glacier microorganisms, four DNA extraction methods with differem lysis pat- terns were tested and two screened methods （the Bosshard-Bano method and the Zhou method） were optimized for the most ef- fective form of the filter membrane （cut vs. uncut）, the DNA extraction method, and the precipitation method. The two optimized methods were then compared with the commercial Mo-Bio DNA extraction kit, and the results showed that the kit was generally suitable for extraction of microorganism DNA fi＇om glacier surface snow. Procedurally, it was found that a modified Boss- hard-Bano method （i.e., cutting the filter membrane into pieces, using a specific lysis pattern [lysozyme （5 mg/mL）-protease K （ 1 mg/mL）-CTAB （ 1%）-SDS （ 1%）], performing the extraction only once by chloroform-isoamyl alcohol （24： 1）, and conducting DNA precipitation by pure ethanol） was also an effective and less expensive method for extraction of microorganism DNA from glacier surface snow.展开更多
In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model f...In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.展开更多
The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine...The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.展开更多
This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order...This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order to reduce experimental cost,the D-optimal criteria with constraints are constructed.The turntable positions and the number of test points are chosen to build D-optimal experimental designs.The D-optimal experimental designs are tested by multi-position calibration experiment for tactical-grade DTGs.Test results show that,with the same cost,the fit uncertainty is reduced by about 50% by using the D-optimal 8-position experimental procedure,compared to using a defacto standard experimental procedure in ANSI/IEEE Std 813-1988.Furthermore,the new experimental procedure almost achieves optimal accuracy with only 12-position which is half the cost of the widely adopted 24-position experimental procedure for achieving optimal accuracy.展开更多
The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimiz...The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.展开更多
The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practi...The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practice.An improved optimal elemental method is presented that defines a new objective function,and as a byproduct,circumvents the need for mass normalized modal shapes,which are also not readily available in practice.To solve the group of nonlinear equations created by the improved optimal method,the Lagrange multiplier method and Matlab function fmincon are employed.To deal with actual complex structures, the float-encoding genetic algorithm(FGA)is introduced to enhance the capability of the improved method.Two examples,a 7- degree of freedom(DOF)mass-spring system and a 53-DOF planar frame,respectively,are updated using the improved method. The example results demonstrate the advantages of the improved method over existing optimal methods,and show that the genetic algorithm is an effective way to update the models used for actual complex structures.展开更多
In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is...In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.展开更多
Based on the widely used DRASTIC method, a fuzzy pattern recognition and optimization method was proposed and applied to the fissured-karstic aquifer of Zhangji area for assessing groundwater vulnerability to pollutio...Based on the widely used DRASTIC method, a fuzzy pattern recognition and optimization method was proposed and applied to the fissured-karstic aquifer of Zhangji area for assessing groundwater vulnerability to pollution. The result is compared with DRASTIC method. It is shown that by taking the fuzziness into consideration, the fuzzy pattern recognition and optimization method reflects more efficiently the fuzzy nature of the groundwater vulnerability to pollution and is more applicable in reality.展开更多
Many effective optimization algorithms require partial derivatives of objective functions,while some optimization problems'objective functions have no derivatives.According to former research studies,some search d...Many effective optimization algorithms require partial derivatives of objective functions,while some optimization problems'objective functions have no derivatives.According to former research studies,some search directions are obtained using the quadratic hypothesis of objective functions.Based on derivatives,quadratic function assumptions,and directional derivatives,the computational formulas of numerical first-order partial derivatives,second-order partial derivatives,and numerical second-order mixed partial derivatives were constructed.Based on the coordinate transformation relation,a set of orthogonal vectors in the fixed coordinate system was established according to the optimization direction.A numerical algorithm was proposed,taking the second order approximation direction as an example.A large stepsize numerical algorithm based on coordinate transformation was proposed.Several algorithms were validated by an unconstrained optimization of the two-dimensional Rosenbrock objective function.The numerical second order approximation direction with the numerical mixed partial derivatives showed good results.Its calculated amount is 0.2843%of that of without second-order mixed partial derivative.In the process of rotating the local coordinate system 360°,because the objective function is more complex than the quadratic function,if the numerical direction derivative is used instead of the analytic partial derivative,the optimization direction varies with a range of 103.05°.Because theoretical error is in the numerical negative gradient direction,the calculation with the coordinate transformation is 94.71%less than the calculation without coordinate transformation.If there is no theoretical error in the numerical negative gradient direction or in the large-stepsize numerical optimization algorithm based on the coordinate transformation,the sawtooth phenomenon occurs.When each numerical mixed partial derivative takes more than one point,the optimization results cannot be improved.The numerical direction based on the quadratic hypothesis only requires the objective function to be obtained,but does not require derivability and does not take into account truncation error and rounding error.Thus,the application scopes of many optimization methods are extended.展开更多
ICM （Independent Continuous Mapping） method can solve topological optimization problems with the minimized weight as the objective and subjected to displacement constraints. To get a clearer topological configuratio...ICM （Independent Continuous Mapping） method can solve topological optimization problems with the minimized weight as the objective and subjected to displacement constraints. To get a clearer topological configuration, by introducing the discrete condition of topological variables and integrating with the original objective, an optimal model with multi-objectives is formulated to make the topological variables approach 0 or 1 as near as possible, and the model reduces the effect of deleting rate on the result. The image-filtering method is employed to eliminate the checkerboard patterns and mesh dependence that occurred in the topology optimization of a continuum structure. The computational efficiency is enhanced through selecting quasi-active displacement constraints and a design region. Numerical examples indicate that this algorithm is robust and practicable, though the number of iterations is slightly increased with respect to the original algorithm.展开更多
Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verifi...Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verified the performance of the optimized airfoil.Using CRA09-A as the baseline rotor airfoil,the CRA09-B optimized rotor airfoil was designed successfully.Combined with the foundation of high-precision rotor airfoil stationary test technology,the CRA09-B and CRA09-A rotor airfoils were tested in the S3 MA high-speed wind tunnel of ONERA.In order to correct the aerodynamic data,a single parameter linear wall pressure method is used to consider the tunnel effects.The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable,and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.展开更多
基金Supported by the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
基金the project PID2022-139202OB-I00Neural Networks and Optimization Techniques for the Design and Safe Maintenance of Transportation Infrastructures:Volcanic Rock Geotechnics and Slope Stability(IA-Pyroslope),funded by the Spanish State Research Agency of the Ministry of Science,Innovation and Universities of Spain and the European Regional Development Fund,MCIN/AEI/10.13039/501100011033/FEDER,EU。
文摘Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.
基金Supported by National Natural Science Foundation of China(Grant No.51975025)National Key Research and Development Program of China(Grant No.2019YFB2004500)。
文摘The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.
基金sponsored by the National Natural Science Foundation of China(Grant Nos.41930971,42330111,and 42405061)the National Key Scientific and Technological Infrastructure project“Earth System Numerical Simulation Facility”(Earth Lab).
文摘Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly efficient calculations for O-CNOPs are still challenging in the field of ensemble forecasting.In this study,we combine a gradient-based iterative idea with the Gram‒Schmidt orthogonalization,and propose an iterative optimization method to compute O-CNOPs.This method is different from the original sequential optimization method,and allows parallel computations of O-CNOPs,thus saving a large amount of computational time.We evaluate this method by using the Lorenz-96 model on the basis of the ensemble forecasting ability achieved and on the time consumed for computing O-CNOPs.The results demonstrate that the parallel iterative method causes O-CNOPs to yield reliable ensemble members and to achieve ensemble forecasting skills similar to or even slightly higher than those produced by the sequential method.Moreover,the parallel method significantly reduces the computational time for O-CNOPs.Therefore,the parallel iterative method provides a highly effective and efficient approach for calculating O-CNOPs for ensemble forecasts.Expectedly,it can play an important role in the application of the O-CNOPs to realistic ensemble forecasts for high-impact weather and climate events.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
文摘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.
基金supported by the National Natural Science Foundation of China(NSFC,Nos.51772205,51572192,51772208,51472179)the General Program of Municipal Natural Science Foundation of Tianjin(Nos.17JCYBJC17000,17JCYBJC22700)。
文摘Lithium-ion batteries(LIBs)have evolved into the mainstream power source of ene rgy sto rage equipment by reason of their advantages such as high energy density,high power,long cycle life and less pollution.With the expansion of their applications in deep-sea exploration,aerospace and military equipment,special working conditions have placed higher demands on the low-temperature performance of LIBs.However,at low temperatures,the severe polarization and inferior electrochemical activity of electrode materials cause the acute capacity fading upon cycling,which greatly hindered the further development of LIBs.In this review,we summarize the recent important progress of LIBs in low-temperature operations and introduce the key methods and the related action mechanisms for enhancing the capacity of the various cathode and anode materials.It aims to promote the development of high-performance electrode materials and broaden the application range of LIBs.
文摘This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.
文摘Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms （EAs） and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.
基金supported by the National Natural Science Foundation of China (Nos. 40825017, 40576001 and 31100369)
文摘In order to thoroughly investigate the diversity of glacier microorganisms, four DNA extraction methods with differem lysis pat- terns were tested and two screened methods （the Bosshard-Bano method and the Zhou method） were optimized for the most ef- fective form of the filter membrane （cut vs. uncut）, the DNA extraction method, and the precipitation method. The two optimized methods were then compared with the commercial Mo-Bio DNA extraction kit, and the results showed that the kit was generally suitable for extraction of microorganism DNA fi＇om glacier surface snow. Procedurally, it was found that a modified Boss- hard-Bano method （i.e., cutting the filter membrane into pieces, using a specific lysis pattern [lysozyme （5 mg/mL）-protease K （ 1 mg/mL）-CTAB （ 1%）-SDS （ 1%）], performing the extraction only once by chloroform-isoamyl alcohol （24： 1）, and conducting DNA precipitation by pure ethanol） was also an effective and less expensive method for extraction of microorganism DNA from glacier surface snow.
基金This work was supported by the National Natural Science Foundation of China(10071037)
文摘In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.
基金funded by National Nature Science Foundation of China(Nos.51776010 and 91860205)supported by the Academic Excellence Foundation of BUAA for PhD Students,China。
文摘The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.
基金National Natural Science Foundation of China (61071014)National Basic Research Program of China(2009CB72400201)
文摘This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order to reduce experimental cost,the D-optimal criteria with constraints are constructed.The turntable positions and the number of test points are chosen to build D-optimal experimental designs.The D-optimal experimental designs are tested by multi-position calibration experiment for tactical-grade DTGs.Test results show that,with the same cost,the fit uncertainty is reduced by about 50% by using the D-optimal 8-position experimental procedure,compared to using a defacto standard experimental procedure in ANSI/IEEE Std 813-1988.Furthermore,the new experimental procedure almost achieves optimal accuracy with only 12-position which is half the cost of the widely adopted 24-position experimental procedure for achieving optimal accuracy.
基金This study was funded by the National Natural Science Foundation of China(Grant 11272167).
文摘The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.
基金The China Hi-Tech R&D Program(863 Program) Project Number 2001AA602023
文摘The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practice.An improved optimal elemental method is presented that defines a new objective function,and as a byproduct,circumvents the need for mass normalized modal shapes,which are also not readily available in practice.To solve the group of nonlinear equations created by the improved optimal method,the Lagrange multiplier method and Matlab function fmincon are employed.To deal with actual complex structures, the float-encoding genetic algorithm(FGA)is introduced to enhance the capability of the improved method.Two examples,a 7- degree of freedom(DOF)mass-spring system and a 53-DOF planar frame,respectively,are updated using the improved method. The example results demonstrate the advantages of the improved method over existing optimal methods,and show that the genetic algorithm is an effective way to update the models used for actual complex structures.
基金supported in part by the CNRST Morocco,the Volkswagen Foundation:Grant number I/79315Hydromed project
文摘In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.
基金Project (No. ICA4-CT-2001-10039) supported by Manporivers(Management policies for priority water pollutants and their effects onfoods and human health: general methodology and application toChinese river basins)
文摘Based on the widely used DRASTIC method, a fuzzy pattern recognition and optimization method was proposed and applied to the fissured-karstic aquifer of Zhangji area for assessing groundwater vulnerability to pollution. The result is compared with DRASTIC method. It is shown that by taking the fuzziness into consideration, the fuzzy pattern recognition and optimization method reflects more efficiently the fuzzy nature of the groundwater vulnerability to pollution and is more applicable in reality.
基金supported in part by the Teaching Reform Research Foundation of Shengli College in China University of Petroleum(East China)(JG201725)the Natural Science Foundation Shandong Province of China(ZR2018PEE009)the Project of Science and Technology of Shandong Universities in China(J17KA044,J17KB061)。
文摘Many effective optimization algorithms require partial derivatives of objective functions,while some optimization problems'objective functions have no derivatives.According to former research studies,some search directions are obtained using the quadratic hypothesis of objective functions.Based on derivatives,quadratic function assumptions,and directional derivatives,the computational formulas of numerical first-order partial derivatives,second-order partial derivatives,and numerical second-order mixed partial derivatives were constructed.Based on the coordinate transformation relation,a set of orthogonal vectors in the fixed coordinate system was established according to the optimization direction.A numerical algorithm was proposed,taking the second order approximation direction as an example.A large stepsize numerical algorithm based on coordinate transformation was proposed.Several algorithms were validated by an unconstrained optimization of the two-dimensional Rosenbrock objective function.The numerical second order approximation direction with the numerical mixed partial derivatives showed good results.Its calculated amount is 0.2843%of that of without second-order mixed partial derivative.In the process of rotating the local coordinate system 360°,because the objective function is more complex than the quadratic function,if the numerical direction derivative is used instead of the analytic partial derivative,the optimization direction varies with a range of 103.05°.Because theoretical error is in the numerical negative gradient direction,the calculation with the coordinate transformation is 94.71%less than the calculation without coordinate transformation.If there is no theoretical error in the numerical negative gradient direction or in the large-stepsize numerical optimization algorithm based on the coordinate transformation,the sawtooth phenomenon occurs.When each numerical mixed partial derivative takes more than one point,the optimization results cannot be improved.The numerical direction based on the quadratic hypothesis only requires the objective function to be obtained,but does not require derivability and does not take into account truncation error and rounding error.Thus,the application scopes of many optimization methods are extended.
基金supported by the National Natural Science Foundation of China(10472003)Beijing Natural Science(3002002)+1 种基金Beijing Educational Committee Foundations(KM200410005019)Suspensofled by American MSC Company.
文摘ICM （Independent Continuous Mapping） method can solve topological optimization problems with the minimized weight as the objective and subjected to displacement constraints. To get a clearer topological configuration, by introducing the discrete condition of topological variables and integrating with the original objective, an optimal model with multi-objectives is formulated to make the topological variables approach 0 or 1 as near as possible, and the model reduces the effect of deleting rate on the result. The image-filtering method is employed to eliminate the checkerboard patterns and mesh dependence that occurred in the topology optimization of a continuum structure. The computational efficiency is enhanced through selecting quasi-active displacement constraints and a design region. Numerical examples indicate that this algorithm is robust and practicable, though the number of iterations is slightly increased with respect to the original algorithm.
基金supported by the National Natural Science Foundation of China(No.11902335)。
文摘Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verified the performance of the optimized airfoil.Using CRA09-A as the baseline rotor airfoil,the CRA09-B optimized rotor airfoil was designed successfully.Combined with the foundation of high-precision rotor airfoil stationary test technology,the CRA09-B and CRA09-A rotor airfoils were tested in the S3 MA high-speed wind tunnel of ONERA.In order to correct the aerodynamic data,a single parameter linear wall pressure method is used to consider the tunnel effects.The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable,and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.
