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 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.展开更多
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.展开更多
This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the...This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the most efficient methodology for computing exact expressions of sensitivities, of any order, of model responses with respect to features of model parameters and, subsequently, with respect to the model’s uncertain parameters, boundaries, and internal interfaces. The unparalleled efficiency and accuracy of the n<sup>th</sup>-FASAM-N methodology stems from the maximal reduction of the number of adjoint computations (which are considered to be “large-scale” computations) for computing high-order sensitivities. When applying the n<sup>th</sup>-FASAM-N methodology to compute the second- and higher-order sensitivities, the number of large-scale computations is proportional to the number of “model features” as opposed to being proportional to the number of model parameters (which are considerably more than the number of features).When a model has no “feature” functions of parameters, but only comprises primary parameters, the n<sup>th</sup>-FASAM-N methodology becomes identical to the extant n<sup>th</sup> CASAM-N (“n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems”) methodology. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are formulated in linearly increasing higher-dimensional Hilbert spaces as opposed to exponentially increasing parameter-dimensional spaces thus overcoming the curse of dimensionality in sensitivity analysis of nonlinear systems. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N are incomparably more efficient and more accurate than any other methods (statistical, finite differences, etc.) for computing exact expressions of response sensitivities of any order with respect to the model’s features and/or primary uncertain parameters, boundaries, and internal interfaces.展开更多
Sensitivity analysis(SA) has been widely used to screen out a small number of sensitive parameters for model outputs from all adjustable parameters in weather and climate models, helping to improve model predictions b...Sensitivity analysis(SA) has been widely used to screen out a small number of sensitive parameters for model outputs from all adjustable parameters in weather and climate models, helping to improve model predictions by tuning the parameters. However, most parametric SA studies have focused on a single SA method and a single model output evaluation function, which makes the screened sensitive parameters less comprehensive. In addition, qualitative SA methods are often used because simulations using complex weather and climate models are time-consuming. Unlike previous SA studies, this research has systematically evaluated the sensitivity of parameters that affect precipitation and temperature simulations in the Weather Research and Forecasting(WRF) model using both qualitative and quantitative global SA methods. In the SA studies, multiple model output evaluation functions were used to conduct various SA experiments for precipitation and temperature. The results showed that five parameters(P3, P5, P7, P10, and P16) had the greatest effect on precipitation simulation results and that two parameters(P7 and P10) had the greatest effect for temperature. Using quantitative SA, the two-way interactive effect between P7 and P10 was also found to be important, especially for precipitation. The microphysics scheme had more sensitive parameters for precipitation, and P10(the multiplier for saturated soil water content) was the most sensitive parameter for both precipitation and temperature. From the ensemble simulations, preliminary results indicated that the precipitation and temperature simulation accuracies could be improved by tuning the respective sensitive parameter values, especially for simulations of moderate and heavy rain.展开更多
基金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.
基金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.
文摘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.
文摘This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the most efficient methodology for computing exact expressions of sensitivities, of any order, of model responses with respect to features of model parameters and, subsequently, with respect to the model’s uncertain parameters, boundaries, and internal interfaces. The unparalleled efficiency and accuracy of the n<sup>th</sup>-FASAM-N methodology stems from the maximal reduction of the number of adjoint computations (which are considered to be “large-scale” computations) for computing high-order sensitivities. When applying the n<sup>th</sup>-FASAM-N methodology to compute the second- and higher-order sensitivities, the number of large-scale computations is proportional to the number of “model features” as opposed to being proportional to the number of model parameters (which are considerably more than the number of features).When a model has no “feature” functions of parameters, but only comprises primary parameters, the n<sup>th</sup>-FASAM-N methodology becomes identical to the extant n<sup>th</sup> CASAM-N (“n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems”) methodology. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are formulated in linearly increasing higher-dimensional Hilbert spaces as opposed to exponentially increasing parameter-dimensional spaces thus overcoming the curse of dimensionality in sensitivity analysis of nonlinear systems. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N are incomparably more efficient and more accurate than any other methods (statistical, finite differences, etc.) for computing exact expressions of response sensitivities of any order with respect to the model’s features and/or primary uncertain parameters, boundaries, and internal interfaces.
基金supported by the Special Fund for Meteorological Scientific Research in the Public Interest (Grant No. GYHY201506002, CRA40: 40-year CMA global atmospheric reanalysis)the National Basic Research Program of China (Grant No. 2015CB953703)+1 种基金the Intergovernmental Key International S & T Innovation Cooperation Program (Grant No. 2016YFE0102400)the National Natural Science Foundation of China (Grant Nos. 41305052 & 41375139)
文摘Sensitivity analysis(SA) has been widely used to screen out a small number of sensitive parameters for model outputs from all adjustable parameters in weather and climate models, helping to improve model predictions by tuning the parameters. However, most parametric SA studies have focused on a single SA method and a single model output evaluation function, which makes the screened sensitive parameters less comprehensive. In addition, qualitative SA methods are often used because simulations using complex weather and climate models are time-consuming. Unlike previous SA studies, this research has systematically evaluated the sensitivity of parameters that affect precipitation and temperature simulations in the Weather Research and Forecasting(WRF) model using both qualitative and quantitative global SA methods. In the SA studies, multiple model output evaluation functions were used to conduct various SA experiments for precipitation and temperature. The results showed that five parameters(P3, P5, P7, P10, and P16) had the greatest effect on precipitation simulation results and that two parameters(P7 and P10) had the greatest effect for temperature. Using quantitative SA, the two-way interactive effect between P7 and P10 was also found to be important, especially for precipitation. The microphysics scheme had more sensitive parameters for precipitation, and P10(the multiplier for saturated soil water content) was the most sensitive parameter for both precipitation and temperature. From the ensemble simulations, preliminary results indicated that the precipitation and temperature simulation accuracies could be improved by tuning the respective sensitive parameter values, especially for simulations of moderate and heavy rain.
