The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessar...The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessary information for model predictive control.Despite the many advantages of MHE,long computation time has limited its applications for system-level models of liquid propellant rocket engines.To address this issue,we propose an asynchronous MHE method called advanced-multi-step MHE with Noise Covariance Estimation(amsMHE-NCE).This method computes the MHE problem asynchronously to obtain the states and parameters and can be applied to multi-threaded computations.In the background,the state and covariance estimation optimization problems are computed using multiple sampling times.In real-time,sensitivity is used to quickly approximate state and parameter estimates.A covariance estimation method is developed using sensitivity to avoid redundant MHE problem calculations in case of sensor degradation during engine reuse.The amsMHE-NCE is validated through three cases based on the space shuttle main engine system-level model,and we demonstrate that it can provide more accurate real-time estimates of states and parameters compared to other commonly used estimation methods.展开更多
The existing trajectory clustering （TRACLUS） is sensitive to the input parameters c and MinLns. The parameter value is changed a little, but cluster results are entirely different. Aiming at this vulnerability, a sh...The existing trajectory clustering （TRACLUS） is sensitive to the input parameters c and MinLns. The parameter value is changed a little, but cluster results are entirely different. Aiming at this vulnerability, a shielding parameters sensitivity trajectory cluster （SPSTC） algorithm is proposed which is insensitive to the input parameters. Firstly, some definitions about the core distance and reachable distance of line segment are presented, and then the algorithm generates cluster sorting according to the core dis- tance and reachable distance. Secondly, the reachable plots of line segment sets are constructed according to the cluster sorting and reachable distance. Thirdly, a parameterized sequence is extracted according to the reachable plot, and then the final trajectory cluster based on the parameterized sequence is acquired. The parameterized sequence represents the inner cluster structure of trajectory data. Experiments on real data sets and test data sets show that the SPSTC algorithm effectively reduces the sensitivity to the input parameters, meanwhile it can obtain the better quality of the trajectory cluster.展开更多
Ion parameters in electron cyclotron resonance （ECR） microwave plasma were measured by ion sensitive probe and were compared with the electron parameters obtained by double Langmuir probe. The effects of gas pressur...Ion parameters in electron cyclotron resonance （ECR） microwave plasma were measured by ion sensitive probe and were compared with the electron parameters obtained by double Langmuir probe. The effects of gas pressure and microwave power on the ion temperature and density were analyzed. The spatial distribution of the ion parameters was also investigated by the ion sensitive probes with a tunable radial depth installed on different probe windows along the chamber axis. Results showed that the ion density measured by the ion sensitive probe was in good agreement with the electron density measured by the double Langmuir probe. The influ- ence of gas pressure on the ion parameters was stronger than that of microwave power. With the increase in working pressure, the ion temperature decreased monotonously with a decreasing rate larger than that at higher pressure. The ion density first increased to a peak （42.3~ 101~ cm-3） at 1 Pa and then decreased. The ion temperature and density increased little with the increase in the microwave power from 400 W to 800 W, The plasma far away from the resonant point is found to be radially uniform.展开更多
The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequentl...The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.展开更多
Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity ana...Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity analysis of hydrological model is a key step in model uncertainty quantification, which can identify the dominant parameters, reduce the model calibration uncertainty, and enhance the model optimization efficiency. There are, however, some shortcomings in classical approaches, including the long duration of time and high computation cost required to quantitatively assess the sensitivity of a multiple-parameter hydrological model. For this reason, a two-step statistical evaluation framework using global techniques is presented. It is based on （1） a screening method （Morris） for qualitative ranking of parameters, and （2） a variance-based method integrated with a meta-model for quantitative sensitivity analysis, i.e., the Sobol method integrated with the response surface model （RSMSobol）. First, the Morris screening method was used to qualitatively identify the parameters＇ sensitivity, and then ten parameters were selected to quantify the sensitivity indices. Subsequently, the RSMSobol method was used to quantify the sensitivity, i.e., the first-order and total sensitivity indices based on the response surface model （RSM） were calculated. The RSMSobol method can not only quantify the sensitivity, but also reduce the computational cost, with good accuracy compared to the classical approaches. This approach will be effective and reliable in the global sensitivity analysis of a complex large-scale distributed hydrological model.展开更多
The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which compris...The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.展开更多
Reducing the error of sensitive parameters by studying the parameters sensitivity can reduce the uncertainty of the model,while simulating double-gyre variation in Regional Ocean Modeling System(ROMS).Conditional Nonl...Reducing the error of sensitive parameters by studying the parameters sensitivity can reduce the uncertainty of the model,while simulating double-gyre variation in Regional Ocean Modeling System(ROMS).Conditional Nonlinear Optimal Perturbation related to Parameter(CNOP-P)is an effective method of studying the parameters sensitivity,which represents a type of parameter error with maximum nonlinear development at the prediction time.Intelligent algorithms have been widely applied to solving Conditional Nonlinear Optimal Perturbation(CNOP).In the paper,we proposed an improved simulated annealing(SA)algorithm to solve CNOP-P to get the optimal parameters error,studied the sensitivity of the single parameter and the combination of multiple parameters and verified the effect of reducing the error of sensitive parameters on reducing the uncertainty of model simulation.Specifically,we firstly found the non-period oscillation of kinetic energy time series of double gyre variation,then extracted two transition periods,which are respectively from high energy to low energy and from low energy to high energy.For every transition period,three parameters,respectively wind amplitude(WD),viscosity coefficient(VC)and linear bottom drag coefficient(RDRG),were studied by CNOP-P solved with SA algorithm.Finally,for sensitive parameters,their effect on model simulation is verified.Experiments results showed that the sensitivity order is WD>VC>>RDRG,the effect of the combination of multiple sensitive parameters is greater than that of single parameter superposition and the reduction of error of sensitive parameters can effectively reduce model prediction error which confirmed the importance of sensitive parameters analysis.展开更多
In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The for...In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate（ZooAE_N）,zooplankton basal metabolism rate（ZooBM）,maximum specific growth rate of zooplankton（μ_（20）） and maximum chlorophyll to carbon ratio（Chl2C_m） are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_（20＇）,half-saturation for phytoplankton NO_3 uptake（K_（NO_3）） and remineralization rate of small detritusN（SDeRRN） are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^（-1）.The K_（NO_3） and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5（mmol/m^3）^（-1）,and interval of SEdRRN is from 0.01 and 0.1 d^（-1）.The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.展开更多
In order to develop a seasonal snow model of land surface process as accurately as possible for climatic study. it is necessary to fully understand the effects of important snow internal processes and interaction with...In order to develop a seasonal snow model of land surface process as accurately as possible for climatic study. it is necessary to fully understand the effects of important snow internal processes and interaction with air and to get an insight into influence of several relevant parameterization schemes with parameters' uncertainty to some degree. Using the snow model (SAST) developed by first author and other one and some useful field observation data, this paper has conducted a series of sensitivity studies on the parameterization schemes. They are relative to compaction process, snow thermal conduction, methodology of layering snow pack and to key parameters such as snow albedo, water holding capacity. Then, based on the results from the sensitivity studies, some useful conclusions for snow cover model improvement are obtained from the analysis of the results.展开更多
Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately ...Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately predict the productivity of hydrate-bearing sediments(HBS).The multi-phase seepage parameters of HBS include permeability,porosity,which is closely related to permeability,and hydrate saturation,which has a direct impact on hydrate content.Existing research has shown that these multi-phase seepage parameters have a great impact on HBS productivity.Permeability directly affects the transmission of pressure-drop and discharge of methane gas,porosity and initial hydrate saturation affect the amount of hydrate decomposition and transmission process of pressure-drop,and also indirectly affect temperature variation of the reservoir.Considering the spatial heterogeneity of multi-phase seepage parameters,a depressurization production model with layered heterogeneity is established based on the clayey silt hydrate reservoir at W11 station in the Shenhu Sea area of the South China Sea.Tough+Hydrate software was used to calculate the production model;the process of gas production and seepage parameter evolution under different multi-phase seepage conditions were obtained.A sensitivity analysis of the parameters affecting the reservoir productivity was conducted so that:(a)a HBS model with layered heterogeneity can better describe the transmission process of pressure and thermal compensation mechanism of hydrate reservoir;(b)considering the multi-phase seepage parameter heterogeneity,the influence degrees of the parameters on HBS productivity were permeability,porosity and initial hydrate saturation,in order from large to small,and the influence of permeability was significantly greater than that of other parameters;(c)the production potential of the clayey silt reservoir should not only be determined by hydrate content or seepage capacity,but also by the comprehensive effect of the two;and(d)time scales need to be considered when studying the effects of changes in multi-phase seepage parameters on HBS productivity.展开更多
Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about di...Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.展开更多
Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to ...Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.展开更多
Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the sel...Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the selection and application of higher sensitivity seismic attribute parameters to improve the ability to discriminate fluid and lithology. To better integrate with seismic interpretation, we construct a template to highlight rock physics parameters in sensitivity space, providing guidance for the quantitative seismic interpretation of hydrocarbon-bearing reservoirs.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
Sensitivity analysis of thermal equilibrium parameters in the reservoir module of MIKE 11 model was conducted for the Wuxikou Reservoir in Jiangxi Province of China in order to apply the module to the environmental im...Sensitivity analysis of thermal equilibrium parameters in the reservoir module of MIKE 11 model was conducted for the Wuxikou Reservoir in Jiangxi Province of China in order to apply the module to the environmental impact assessment to accurately predict water temperature of reservoirs.Results showed that radiation parameter A and evaporation-first parameter were much more sensitive than other parameters.The values of the radiation parameter A ranged from 0.10 to 0.34.The values of evaporation-first parameter varied from 0 to 10.The sensitivity of solar absorption parameters was less than that of evaporation parameter,of which light attenuation values ranged from 0.5 to 0.7,and this parameter would not impact model results if it was more than 2.Constants in Beer's law ranged from 0.2 to 0.7.Radiation parameter B was not more sensitive than evaporation parameter and its reasonable range was higher than 0.48.The fitting curves showed consistent changing tendency for these parameters within the reasonable ranges.Additionally,all the thermal equilibrium parameters had much more important effects on surface water temperature than deep water temperature.Moreover,if no observed data could be obtained,the local empirical value would be used to input to the MIKE 11 model to simulate the changes in the discharged outflow-water temperature qualitatively.展开更多
Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration numbe...Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration number, computation time and convergence in a mine ventilation network analysis, were investigated based on 5 mine ventilation systems. The results show that a higher computation accuracy greatly influences the iteration number. When the accuracy reaches 10-6m3·s-1 for solving a complicated mine ventilation network, the running time is too long though a high-speed computer is used. The preliminary value of airflow rate in the range of 1100m3·s-1 has little effects the iteration number. The structure of network also has some effect on the iteration number.展开更多
In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through...In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.展开更多
Accurate measurement of helicopter rotor motion parameters(flap,lead-lag,torsion,and azimuth angles)is essential for rotor blade design,helicopter dynamics modeling,and flight safety and health monitoring.However,the ...Accurate measurement of helicopter rotor motion parameters(flap,lead-lag,torsion,and azimuth angles)is essential for rotor blade design,helicopter dynamics modeling,and flight safety and health monitoring.However,the existing methods face challenges in testing equipment installation,calibration,and data transmission,resulting in limited reports on real-time in-flight measurements of blade motion parameters.This paper proposes a non-contact optoelectronic method based on two-dimensional position-sensitive detectors for in-flight measurement and a ground calibration system to obtain real-time rotor motion parameters during helicopter flight.The proposed method establishes the time evolution relationship of rotor motion parameters and verifies the performance of the in-flight measurement system regarding measurement resolution and accuracy through the construction of a blade motion posture experimental platform.The proposed method has been applied to the flight measurement of a medium-sized single-rotor helicopter,and the obtained results have been compared with theoretical analysis outcomes.Furthermore,this paper examines the characteristics of blade motion parameters during flight and discusses the challenges and potential solutions for measuring rotor motion parameters during helicopter flight using the proposed method.展开更多
That a model has sensitivity responses to parameter uncertainties is a key concept in implementing model parameter es- timation using filtering theory and methodology. Depending on the nature of associated physics and...That a model has sensitivity responses to parameter uncertainties is a key concept in implementing model parameter es- timation using filtering theory and methodology. Depending on the nature of associated physics and characteristic variability of the fluid in a coupled system, the response time scales of a model to parameters can be different, from hourly to decadal. Unlike state estimation, where the update frequency is usually linked with observational frequency, the update frequency for parameter estimation must be associated with the time scale of the model sensitivity response to the parameter being esti- mated. Here, with a simple coupled model, the impact of model sensitivity response time scales on coupled model parameter estimation is studied. The model includes characteristic synoptic to decadal scales by coupling a long-term varying deep ocean with a slow-varying upper ocean forced by a chaotic atmosphere. Results show that, using the update frequency deter- mined by the model sensitivity response time scale, both the reliability and quality of parameter estimation can be improved significantly, and thus the estimated parameters make the model more consistent with the observation. These simple model results provide a guideline for when real observations are used to optimize the parameters in a coupled general circulation model for improving climate analysis and prediction initialization.展开更多
This paper presents an analysis method, based on MacCormack＇s technique, for the evaluation of the time domain sensitivity of distributed parameter elements in high-speed circuit networks. Sensitivities can be calcul...This paper presents an analysis method, based on MacCormack＇s technique, for the evaluation of the time domain sensitivity of distributed parameter elements in high-speed circuit networks. Sensitivities can be calculated from electrical and physical parameters of the distributed parameter elements. The proposed method is a direct numerical method of time-space discretization and does not require complicated mathematical deductive process. Therefore, it is very convenient to program this method. It can be applied to sensitivity analysis of general transmission lines in linear or nonlinear circuit networks. The proposed method is second-order-accurate. Numerical experiment is presented to demonstrate its accuracy and efficiency.展开更多
基金supported by the National Natural Science Foundation of China(Nos.62120106003 and 62173301)。
文摘The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessary information for model predictive control.Despite the many advantages of MHE,long computation time has limited its applications for system-level models of liquid propellant rocket engines.To address this issue,we propose an asynchronous MHE method called advanced-multi-step MHE with Noise Covariance Estimation(amsMHE-NCE).This method computes the MHE problem asynchronously to obtain the states and parameters and can be applied to multi-threaded computations.In the background,the state and covariance estimation optimization problems are computed using multiple sampling times.In real-time,sensitivity is used to quickly approximate state and parameter estimates.A covariance estimation method is developed using sensitivity to avoid redundant MHE problem calculations in case of sensor degradation during engine reuse.The amsMHE-NCE is validated through three cases based on the space shuttle main engine system-level model,and we demonstrate that it can provide more accurate real-time estimates of states and parameters compared to other commonly used estimation methods.
基金supported by the National High Technology Research and Development Program of China(863 Program)(2007AA01Z404)the Funding of Jiangsu Provincial Innovation Program for Graduate Education(CXLX110206)
文摘The existing trajectory clustering （TRACLUS） is sensitive to the input parameters c and MinLns. The parameter value is changed a little, but cluster results are entirely different. Aiming at this vulnerability, a shielding parameters sensitivity trajectory cluster （SPSTC） algorithm is proposed which is insensitive to the input parameters. Firstly, some definitions about the core distance and reachable distance of line segment are presented, and then the algorithm generates cluster sorting according to the core dis- tance and reachable distance. Secondly, the reachable plots of line segment sets are constructed according to the cluster sorting and reachable distance. Thirdly, a parameterized sequence is extracted according to the reachable plot, and then the final trajectory cluster based on the parameterized sequence is acquired. The parameterized sequence represents the inner cluster structure of trajectory data. Experiments on real data sets and test data sets show that the SPSTC algorithm effectively reduces the sensitivity to the input parameters, meanwhile it can obtain the better quality of the trajectory cluster.
基金supported by National Natural Science Foundation of China (No. 10875093)
文摘Ion parameters in electron cyclotron resonance （ECR） microwave plasma were measured by ion sensitive probe and were compared with the electron parameters obtained by double Langmuir probe. The effects of gas pressure and microwave power on the ion temperature and density were analyzed. The spatial distribution of the ion parameters was also investigated by the ion sensitive probes with a tunable radial depth installed on different probe windows along the chamber axis. Results showed that the ion density measured by the ion sensitive probe was in good agreement with the electron density measured by the double Langmuir probe. The influ- ence of gas pressure on the ion parameters was stronger than that of microwave power. With the increase in working pressure, the ion temperature decreased monotonously with a decreasing rate larger than that at higher pressure. The ion density first increased to a peak （42.3~ 101~ cm-3） at 1 Pa and then decreased. The ion temperature and density increased little with the increase in the microwave power from 400 W to 800 W, The plasma far away from the resonant point is found to be radially uniform.
基金The project supported by the Innovative Project of CAS (KJCX-SW-L08)the National Basic Research Program of China(973)
文摘The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.
基金supported by the National Natural Science Foundation of China (Grant No. 41271003)the National Basic Research Program of China (Grants No. 2010CB428403 and 2010CB951103)
文摘Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity analysis of hydrological model is a key step in model uncertainty quantification, which can identify the dominant parameters, reduce the model calibration uncertainty, and enhance the model optimization efficiency. There are, however, some shortcomings in classical approaches, including the long duration of time and high computation cost required to quantitatively assess the sensitivity of a multiple-parameter hydrological model. For this reason, a two-step statistical evaluation framework using global techniques is presented. It is based on （1） a screening method （Morris） for qualitative ranking of parameters, and （2） a variance-based method integrated with a meta-model for quantitative sensitivity analysis, i.e., the Sobol method integrated with the response surface model （RSMSobol）. First, the Morris screening method was used to qualitatively identify the parameters＇ sensitivity, and then ten parameters were selected to quantify the sensitivity indices. Subsequently, the RSMSobol method was used to quantify the sensitivity, i.e., the first-order and total sensitivity indices based on the response surface model （RSM） were calculated. The RSMSobol method can not only quantify the sensitivity, but also reduce the computational cost, with good accuracy compared to the classical approaches. This approach will be effective and reliable in the global sensitivity analysis of a complex large-scale distributed hydrological model.
基金Project(2015B37714)supported by the Fundamental Research Funds for the Central Universities of ChinaProject(51605005)supported by the National Natural Science Foundation of China+1 种基金Project(ZK16-03-03)supported by the Open Foundation of Jiangsu Wind Technology Center,ChinaProject([2013]56)supported by the First Group of 2011 Plan of Jiangsu Province,China
文摘The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.
基金Supported by the National Natural Science Foundation of China(No.41405097)the Fundamental Research Funds for the Central Universities of China in 2017
文摘Reducing the error of sensitive parameters by studying the parameters sensitivity can reduce the uncertainty of the model,while simulating double-gyre variation in Regional Ocean Modeling System(ROMS).Conditional Nonlinear Optimal Perturbation related to Parameter(CNOP-P)is an effective method of studying the parameters sensitivity,which represents a type of parameter error with maximum nonlinear development at the prediction time.Intelligent algorithms have been widely applied to solving Conditional Nonlinear Optimal Perturbation(CNOP).In the paper,we proposed an improved simulated annealing(SA)algorithm to solve CNOP-P to get the optimal parameters error,studied the sensitivity of the single parameter and the combination of multiple parameters and verified the effect of reducing the error of sensitive parameters on reducing the uncertainty of model simulation.Specifically,we firstly found the non-period oscillation of kinetic energy time series of double gyre variation,then extracted two transition periods,which are respectively from high energy to low energy and from low energy to high energy.For every transition period,three parameters,respectively wind amplitude(WD),viscosity coefficient(VC)and linear bottom drag coefficient(RDRG),were studied by CNOP-P solved with SA algorithm.Finally,for sensitive parameters,their effect on model simulation is verified.Experiments results showed that the sensitivity order is WD>VC>>RDRG,the effect of the combination of multiple sensitive parameters is greater than that of single parameter superposition and the reduction of error of sensitive parameters can effectively reduce model prediction error which confirmed the importance of sensitive parameters analysis.
基金The National Natural Science Foundation of China under contract Nos 41206023,41222038 and 41076011the National Basic Research Project(973 Program)of China under contract No.2011CB403606+2 种基金the China-Korea Joint Ocean Research Center"Cooperation on the Development of Basic Technologies for the Yellow Sea and East China Sea Operational Oceanographic System(YOOS)"the Public Science and Technology Research Funds Projects of Ocean under contrcat No.201205018the"Strategic Priority Research Program"of the Chinese Academy of Sciences,under contract No.XDA01020304
文摘In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate（ZooAE_N）,zooplankton basal metabolism rate（ZooBM）,maximum specific growth rate of zooplankton（μ_（20）） and maximum chlorophyll to carbon ratio（Chl2C_m） are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_（20＇）,half-saturation for phytoplankton NO_3 uptake（K_（NO_3）） and remineralization rate of small detritusN（SDeRRN） are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^（-1）.The K_（NO_3） and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5（mmol/m^3）^（-1）,and interval of SEdRRN is from 0.01 and 0.1 d^（-1）.The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.
基金This work is financially supported by 1) National Key Programme for Developing Basic Sciences.G1998040900-Part 1, 2) NSF (key
文摘In order to develop a seasonal snow model of land surface process as accurately as possible for climatic study. it is necessary to fully understand the effects of important snow internal processes and interaction with air and to get an insight into influence of several relevant parameterization schemes with parameters' uncertainty to some degree. Using the snow model (SAST) developed by first author and other one and some useful field observation data, this paper has conducted a series of sensitivity studies on the parameterization schemes. They are relative to compaction process, snow thermal conduction, methodology of layering snow pack and to key parameters such as snow albedo, water holding capacity. Then, based on the results from the sensitivity studies, some useful conclusions for snow cover model improvement are obtained from the analysis of the results.
基金supported by the National Natural Science Foundation of China(Grant Nos.42276224,and 42206230)the Jilin Scientific and Technological Development Program(Grant No.20190303083SF)+2 种基金the International Cooperation Key Laboratory of Underground Energy Development and Geological Restoration(Grant No.YDZJ202102CXJD014)the Interdisciplinary Integration and Innovation Project of JLU(Grant No.JLUXKJC2021ZZ18)the Graduate Innovation Fund of Jilin University(Grant No.2023CX100)。
文摘Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately predict the productivity of hydrate-bearing sediments(HBS).The multi-phase seepage parameters of HBS include permeability,porosity,which is closely related to permeability,and hydrate saturation,which has a direct impact on hydrate content.Existing research has shown that these multi-phase seepage parameters have a great impact on HBS productivity.Permeability directly affects the transmission of pressure-drop and discharge of methane gas,porosity and initial hydrate saturation affect the amount of hydrate decomposition and transmission process of pressure-drop,and also indirectly affect temperature variation of the reservoir.Considering the spatial heterogeneity of multi-phase seepage parameters,a depressurization production model with layered heterogeneity is established based on the clayey silt hydrate reservoir at W11 station in the Shenhu Sea area of the South China Sea.Tough+Hydrate software was used to calculate the production model;the process of gas production and seepage parameter evolution under different multi-phase seepage conditions were obtained.A sensitivity analysis of the parameters affecting the reservoir productivity was conducted so that:(a)a HBS model with layered heterogeneity can better describe the transmission process of pressure and thermal compensation mechanism of hydrate reservoir;(b)considering the multi-phase seepage parameter heterogeneity,the influence degrees of the parameters on HBS productivity were permeability,porosity and initial hydrate saturation,in order from large to small,and the influence of permeability was significantly greater than that of other parameters;(c)the production potential of the clayey silt reservoir should not only be determined by hydrate content or seepage capacity,but also by the comprehensive effect of the two;and(d)time scales need to be considered when studying the effects of changes in multi-phase seepage parameters on HBS productivity.
基金The first author would like to thank the China Scholarship Council(CSC No.201806430001)for the financial support for his study in Germany.The authors are grateful to the anonymous reviewer for the recommendations to improve the paper。
文摘Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.
基金supported by the National Natural Science Foundation of China (11142013, 11172260 and 11072214)the Doctoral Fund of Ministry of Education of China (20110101110016)the Fundamental Research Funds for the Central Universities of China(2011QNA4001)
文摘Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.
文摘Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the selection and application of higher sensitivity seismic attribute parameters to improve the ability to discriminate fluid and lithology. To better integrate with seismic interpretation, we construct a template to highlight rock physics parameters in sensitivity space, providing guidance for the quantitative seismic interpretation of hydrocarbon-bearing reservoirs.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
基金Under the auspices of Research&Development Special Fund for Public Welfare Industry of Ministry of Environmental Protection(No.201309062201309003)
文摘Sensitivity analysis of thermal equilibrium parameters in the reservoir module of MIKE 11 model was conducted for the Wuxikou Reservoir in Jiangxi Province of China in order to apply the module to the environmental impact assessment to accurately predict water temperature of reservoirs.Results showed that radiation parameter A and evaporation-first parameter were much more sensitive than other parameters.The values of the radiation parameter A ranged from 0.10 to 0.34.The values of evaporation-first parameter varied from 0 to 10.The sensitivity of solar absorption parameters was less than that of evaporation parameter,of which light attenuation values ranged from 0.5 to 0.7,and this parameter would not impact model results if it was more than 2.Constants in Beer's law ranged from 0.2 to 0.7.Radiation parameter B was not more sensitive than evaporation parameter and its reasonable range was higher than 0.48.The fitting curves showed consistent changing tendency for these parameters within the reasonable ranges.Additionally,all the thermal equilibrium parameters had much more important effects on surface water temperature than deep water temperature.Moreover,if no observed data could be obtained,the local empirical value would be used to input to the MIKE 11 model to simulate the changes in the discharged outflow-water temperature qualitatively.
基金Project (50474050) supported by the National Natural Science Foundation of China
文摘Depending on the numerical test approach on a computer, the relationships among relevant parameters, eg branch number, node number, mesh number, computation accuracy, preliminary value of airflow rate, iteration number, computation time and convergence in a mine ventilation network analysis, were investigated based on 5 mine ventilation systems. The results show that a higher computation accuracy greatly influences the iteration number. When the accuracy reaches 10-6m3·s-1 for solving a complicated mine ventilation network, the running time is too long though a high-speed computer is used. The preliminary value of airflow rate in the range of 1100m3·s-1 has little effects the iteration number. The structure of network also has some effect on the iteration number.
基金Supported by the National Natural Science Foundation of China(10272109)
文摘In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.
基金the funding provided by the National Helicopter Development Project of China。
文摘Accurate measurement of helicopter rotor motion parameters(flap,lead-lag,torsion,and azimuth angles)is essential for rotor blade design,helicopter dynamics modeling,and flight safety and health monitoring.However,the existing methods face challenges in testing equipment installation,calibration,and data transmission,resulting in limited reports on real-time in-flight measurements of blade motion parameters.This paper proposes a non-contact optoelectronic method based on two-dimensional position-sensitive detectors for in-flight measurement and a ground calibration system to obtain real-time rotor motion parameters during helicopter flight.The proposed method establishes the time evolution relationship of rotor motion parameters and verifies the performance of the in-flight measurement system regarding measurement resolution and accuracy through the construction of a blade motion posture experimental platform.The proposed method has been applied to the flight measurement of a medium-sized single-rotor helicopter,and the obtained results have been compared with theoretical analysis outcomes.Furthermore,this paper examines the characteristics of blade motion parameters during flight and discusses the challenges and potential solutions for measuring rotor motion parameters during helicopter flight using the proposed method.
基金funded by the National Natural Science Foundation of China (Grant No.41676088)the National Key Research and Development Project of China (2016YFC1401800,2017YFC1404100,2017YFC1404102)+1 种基金the Fundamental Research Funds for the Central Universities (HEUCF 041705)the Foundation of the Key Laboratory of Marine Environmental Information Technology
文摘That a model has sensitivity responses to parameter uncertainties is a key concept in implementing model parameter es- timation using filtering theory and methodology. Depending on the nature of associated physics and characteristic variability of the fluid in a coupled system, the response time scales of a model to parameters can be different, from hourly to decadal. Unlike state estimation, where the update frequency is usually linked with observational frequency, the update frequency for parameter estimation must be associated with the time scale of the model sensitivity response to the parameter being esti- mated. Here, with a simple coupled model, the impact of model sensitivity response time scales on coupled model parameter estimation is studied. The model includes characteristic synoptic to decadal scales by coupling a long-term varying deep ocean with a slow-varying upper ocean forced by a chaotic atmosphere. Results show that, using the update frequency deter- mined by the model sensitivity response time scale, both the reliability and quality of parameter estimation can be improved significantly, and thus the estimated parameters make the model more consistent with the observation. These simple model results provide a guideline for when real observations are used to optimize the parameters in a coupled general circulation model for improving climate analysis and prediction initialization.
文摘This paper presents an analysis method, based on MacCormack＇s technique, for the evaluation of the time domain sensitivity of distributed parameter elements in high-speed circuit networks. Sensitivities can be calculated from electrical and physical parameters of the distributed parameter elements. The proposed method is a direct numerical method of time-space discretization and does not require complicated mathematical deductive process. Therefore, it is very convenient to program this method. It can be applied to sensitivity analysis of general transmission lines in linear or nonlinear circuit networks. The proposed method is second-order-accurate. Numerical experiment is presented to demonstrate its accuracy and efficiency.
