Implication of community-level ecophysiological parameterization to modelling ecosystem productivity:a case study across nine contrasting forest sites in eastern China

Parameterization is a critical step in modelling ecosystem dynamics.However,assigning parameter values can be a technical challenge for structurally complex natural plant communities;uncertainties in model simulations often arise from inappropriate model parameterization.Here we compared five methods for defining community-level specific leaf area(SLA)and leaf C:N across nine contrasting forest sites along the North-South Transect of Eastern China,including biomass-weighted average for the entire plant community(AP_BW)and four simplified selective sampling(biomass-weighted average over five dominant tree species[5DT_BW],basal area weighted average over five dominant tree species[5DT_AW],biomass-weighted average over all tree species[AT_BW]and basal area weighted average over all tree species[AT_AW]).We found that the default values for SLA and leaf C:N embedded in the Biome-BGC v4.2 were higher than the five computational methods produced across the nine sites,with deviations ranging from 28.0 to 73.3%.In addition,there were only slight deviations(<10%)between the whole plant community sampling(AP_BW)predicted NPP and the four simplified selective sampling methods,and no significant difference between the predictions of AT_BW and AP_BW except the Shennongjia site.The findings in this study highlights the critical importance of computational strategies for community-level parameterization in ecosystem process modelling,and will support the choice of parameterization methods.

Introduction to Mesh Based Generated Lumped Parameter Models for Electromagnetic Problems using Triangular Elements

This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.

Failure Prediction Modeling of Lithium Ion Battery toward Distributed Parameter Estimation

Lithium ion battery has typical character of distributed parameter system, and can be described precisely by partial differential equations and multi-physics theory because lithium ion battery is a complicated electrochemical energy storage system. A novel failure prediction modeling method of lithium ion battery based on distributed parameter estimation and single particle model is proposed in this work. Lithium ion concentration in the anode of lithium ion battery is an unmeasurable distributed variable. Failure prediction system can estimate lithium ion concentration online, track the failure residual which is the difference between the estimated value and the ideal value. The precaution signal will be triggered when the failure residual is beyond the predefined failure precaution threshold, and the failure countdown prediction module will be activated. The remaining time of the severe failure threshold can be estimated by the failure countdown prediction module according to the changing rate of the failure residual. A simulation example verifies that lithium ion concentration in the anode of lithium ion battery can be estimated exactly and effectively by the failure prediction model. The precaution signal can be triggered reliably, and the remaining time of the severe failure can be forecasted accurately by the failure countdown prediction module.

Ensemble Bayesian method for parameter distribution inference:application to reactor physics

The estimation of model parameters is an important subject in engineering.In this area of work,the prevailing approach is to estimate or calculate these as deterministic parameters.In this study,we consider the model parameters from the perspective of random variables and describe the general form of the parameter distribution inference problem.Under this framework,we propose an ensemble Bayesian method by introducing Bayesian inference and the Markov chain Monte Carlo(MCMC)method.Experiments on a finite cylindrical reactor and a 2D IAEA benchmark problem show that the proposed method converges quickly and can estimate parameters effectively,even for several correlated parameters simultaneously.Our experiments include cases of engineering software calls,demonstrating that the method can be applied to engineering,such as nuclear reactor engineering.

INVESTIGATION ON EFFECT OF DRUG DOSING REGIMENTS ON DRUG DELIVERY IN SOLID TUMOR VIA LUMPED PARAMETER MODELING AND ANIMAL EXPERIMENTS

This work aims to investigate the effects of dosing regiments on drug delivery in solid tumors and to validate them with experiments on rats. The lumped parameter models of pharmacokinetics and of drug delivery in tumor were developed to simulate time courses of average drug concentration （Ct） of tumor interstitium in two types of dosing regiments （i.e., single-shot and triple-shot ones）. The two regiments were performed via antitumor drug, hydroxycamptothecin （HCPT）, on rats, to measure the drug concentration in the tumor. The simulations of the drug concentration in the tumor of the two dosing regiments were conducted and compared with the experimental data on rats. The coefficients in the models were investigated. It is concluded that the triple-shot method is more effective than that of single-shot injection. The present lumped-parameter model is quantitatively competent for drug delivery in solid tumor.

A 3D attention U-Net network and its application in geological model parameterization

To solve the problems of convolutional neural network–principal component analysis(CNN-PCA)in fine description and generalization of complex reservoir geological features,a 3D attention U-Net network was proposed not using a trained C3D video motion analysis model to extract the style of a 3D model,and applied to complement the details of geologic model lost in the dimension reduction of PCA method in this study.The 3D attention U-Net network was applied to a complex river channel sandstone reservoir to test its effects.The results show that compared with CNN-PCA method,the 3D attention U-Net network could better complement the details of geological model lost in the PCA dimension reduction,better reflect the fluid flow features in the original geologic model,and improve history matching results.

Introduction to Mesh Based Generated Lumped Parameter Models for Electromagnetic Problems

This paper is an introduction to mesh based generated reluctance network modeling.An overview of scientific works which led to the development of this approach is first presented.Basic concepts of the approach are then presented in the case of electromagnetic devices.A step-by-step procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB and Scilab environments are also included.

Parametric investigation of railway fastenings into the formation and mitigation of short pitch corrugation

Short pitch corrugation has been a problem for railways worldwide over one century.In this paper,a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insights into corrugation mitigation.A three-dimensional finite element vehicle-track dynamic interaction model is employed,which considers the coupling between the structural dynamics and the contact mechanics,while the damage mechanism is assumed to be differential wear.Various fastening models with different configurations,boundary conditions,and parameters of stiffness and damping are built up and analysed.These models may represent different service stages of fastenings in the field.Besides,the effect of train speeds on corrugation features is studied.The results indicate:(1)Fastening parameters and modelling play an important role in corrugation formation.(2)The fastening longitudinal constraint to the rail is the major factor that determines the corrugation formation.The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components.(3)The strengthening of fastening constraints in the longitudinal dimension helps to mitigate corrugation.Meanwhile,the inner fastening constraint in the lateral direction is necessary for corrugation alleviation.(4)The increase in fastening longitudinal stiffness and damping can reduce the vibration amplitudes of longitudinal compression modes and thus reduce the track corrugation propensity.The simulation in this work can well explain the field corrugation in terms of the occurrence possibility and major wavelength components.It can also explain the field data with respect to the small variation between the corrugation wavelength and train speed,which is caused by frequency selection and jump between rail longitudinal compression modes.

Introducing the nth-Order Features Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (nth-FASAM-N): I. Mathematical Framework

This work presents the “nth-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “nth-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 nth-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 nth-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 nth-FASAM-N methodology becomes identical to the extant nth CASAM-N (“nth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems”) methodology. Both the nth-FASAM-N and the nth-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 nth-FASAM-N and the nth-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.

Introducing the nth-Order Features Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (nth-FASAM-N): II. Illustrative Example

This work highlights the unparalleled efficiency of the “nth-Order Function/ Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (nth-FASAM-N) by considering the well-known Nordheim-Fuchs reactor dynamics/safety model. This model describes a short-time self-limiting power excursion in a nuclear reactor system having a negative temperature coefficient in which a large amount of reactivity is suddenly inserted, either intentionally or by accident. This nonlinear paradigm model is sufficiently complex to model realistically self-limiting power excursions for short times yet admits closed-form exact expressions for the time-dependent neutron flux, temperature distribution and energy released during the transient power burst. The nth-FASAM-N methodology is compared to the extant “nth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (nth-CASAM-N) showing that: (i) the 1st-FASAM-N and the 1st-CASAM-N methodologies are equally efficient for computing the first-order sensitivities;each methodology requires a single large-scale computation for solving the “First-Level Adjoint Sensitivity System” (1st-LASS);(ii) the 2nd-FASAM-N methodology is considerably more efficient than the 2nd-CASAM-N methodology for computing the second-order sensitivities since the number of feature-functions is much smaller than the number of primary parameters;specifically for the Nordheim-Fuchs model, the 2nd-FASAM-N methodology requires 2 large-scale computations to obtain all of the exact expressions of the 28 distinct second-order response sensitivities with respect to the model parameters while the 2nd-CASAM-N methodology requires 7 large-scale computations for obtaining these 28 second-order sensitivities;(iii) the 3rd-FASAM-N methodology is even more efficient than the 3rd-CASAM-N methodology: only 2 large-scale computations are needed to obtain the exact expressions of the 84 distinct third-order response sensitivities with respect to the Nordheim-Fuchs model’s parameters when applying the 3rd-FASAM-N methodology, while the application of the 3rd-CASAM-N methodology requires at least 22 large-scale computations for computing the same 84 distinct third-order sensitivities. Together, the nth-FASAM-N and the nth-CASAM-N methodologies are the most practical methodologies for computing response sensitivities of any order comprehensively and accurately, overcoming the curse of dimensionality in sensitivity analysis.

Is Model Parameter Error Related to a Significant Spring Predictability Barrier for El Nio events? Results from a Theoretical Model

Within a theoretical ENSO model, the authors investigated whether or not the errors superimposed on model parameters could cause a significant ＂spring predictability barrier＂ （SPB） for El Nio events. First, sensitivity experiments were respectively performed to the air-sea coupling parameter, α and the thermocline effect coefficient μ. The results showed that the uncertainties superimposed on each of the two parameters did not exhibit an obvious season-dependent evolution; furthermore, the uncertainties caused a very small prediction error and consequently failed to yield a significant SPB. Subsequently, the conditional nonlinear optimal perturbation （CNOP） approach was used to study the effect of the optimal mode （CNOP-P） of the uncertainties of the two parameters on the SPB and to demonstrate that the CNOP-P errors neither presented a unified season-dependent evolution for different El Nio events nor caused a large prediction error, and therefore did not cause a significant SPB. The parameter errors played only a trivial role in yielding a significant SPB. To further validate this conclusion, the authors investigated the effect of the optimal combined mode （i.e. CNOP error） of initial and model errors on SPB. The results illustrated that the CNOP errors tended to have a significant season-dependent evolution, with the largest error growth rate in the spring, and yielded a large prediction error, inducing a significant SPB. The inference, therefore, is that initial errors, rather than model parameter errors, may be the dominant source of uncertainties that cause a significant SPB for El Nio events. These results indicate that the ability to forecast ENSO could be greatly increased by improving the initialization of the forecast model.

Theoretical Investigations of the Local Structure Distortion and EPR Parameter for Ni2＋-doped Perovskite Fluorides

By analyzing the optical spectra and electron paramagnetic resonance parameter D, the local structure distortion of （NiF6）4- clusters in AMF3 （A=K, Rb; M=Zn, Cd, Ca） and K2ZnF4 series are studied using the complete energy matrix based on the double spin-orbit coupling parameter model for configuration ions in a tetragonal ligand field. The results indicate that the contribution of ligand to spin-orbit coupling interaction should be considered for our studied systems. Moreover, the relationships between D and the spin-obit coupling coefficients as well as the average parameter and the divergent parameter are discussed.

Modeling mechanism and extension of GM (1,1)

Firstly, the research progress of grey model GM （1,1） is summarized, which is divided into three development stages： assimilation, alienation and melting stages. Then, the matrix analysis theory is used to study the modeling mechanism of GM （1,1）, which decomposes the modeling data matrix into raw data transformation matrix, accumulated generating operation matrix and background value selection matrix. The changes of these three matrices are the essential reasons affecting the modeling and the accuracy of GM （1,1）. Finally, the paper proposes a generalization grey model GGM （1,1）, which is a extended form of GM （1,1） and also a unified form of model GM （1,1）, model GM （1,1,α）, stage grey model, hopping grey model, generalized accumulated model, strengthening operator model, weakening operator model and unequal interval model. And the theory and practical significance of the extended model is analyzed.

Extraction of temperature dependences of small-signal model parameters in SiGe HBT HICUM model

In this work, temperature dependences of small-signal model parameters in the SiGe HBT HICUM model are presented. Electrical elements in the small-signal equivalent circuit are first extracted at each temperature, then the temperature dependences are determined by the series of extracted temperature coefficients, based on the established temperature for- mulas for corresponding model parameters. The proposed method is validated by a 1x 0.2 x 16 μm2 SiGe HBT over a wide temperature range （from 218 K to 473 K）, and good matching is obtained between the extracted and modeled resuits. Therefore, we believe that the proposed extraction flow of model parameter temperature dependence is reliable for characterizing the transistor performance and guiding the circuit design over a wide temperature range.

A modified lumped parameter model of distribution transformer winding

The modelling of the distribution transformer winding is the starting point and serves as important basis for the transformer characteristics analysis and the lightning pulse response prediction.A distributed parameters model can depict the winding characteristics accurately,but it requires complex calculations.Lumped parameter model requires less calculations,but its applicable frequency range is not wide.This paper studies the amplitude-frequency characteristics of the lightning wave,compares the transformer modelling methods and finally proposes a modified lumped parameter model,based on the above comparison.The proposed model minimizes the errors provoked by the lumped parameter approximation,and the hyperbolic functions of the distributed parameter model.By this modification it becomes possible to accurately describe the winding characteristics and rapidly obtain the node voltage response.The proposed model can provide theoretical and experimental support to lightning protection of the distribution transformer.

Model Parameter Transfer for Gear Fault Diagnosis under Varying Working Conditions

Gear fault diagnosis technologies have received rapid development and been effectively implemented in many engineering applications.However,the various working conditions would degrade the diagnostic performance and make gear fault diagnosis(GFD)more and more challenging.In this paper,a novel model parameter transfer(NMPT)is proposed to boost the performance of GFD under varying working conditions.Based on the previous transfer strategy that controls empirical risk of source domain,this method further integrates the superiorities of multi-task learning with the idea of transfer learning(TL)to acquire transferable knowledge by minimizing the discrepancies of separating hyperplanes between one specific working condition(target domain)and another(source domain),and then transferring both commonality and specialty parameters over tasks to make use of source domain samples to assist target GFD task when sufficient labeled samples from target domain are unavailable.For NMPT implementation,insufficient target domain features and abundant source domain features with supervised information are fed into NMPT model to train a robust classifier for target GFD task.Related experiments prove that NMPT is expected to be a valuable technology to boost practical GFD performance under various working conditions.The proposed methods provides a transfer learning-based framework to handle the problem of insufficient training samples in target task caused by variable operation conditions.

Soft Sensor for Inputs and Parameters Using Nonlinear Singular State Observer in Chemical Processes

Chemical processes are usually nonlinear singular systems.In this study,a soft sensor using nonlinear singular state observer is established for unknown inputs and uncertain model parameters in chemical processes,which are augmented as state variables.Based on the observability of the singular system,this paper presents a simplified observability criterion under certain conditions for unknown inputs and uncertain model parameters.When the observability is satisfied,the unknown inputs and the uncertain model parameters are estimated online by the soft sensor using augmented nonlinear singular state observer.The riser reactor of fluid catalytic cracking unit is used as an example for analysis and simulation.With the catalyst circulation rate as the only unknown input without model error,one temperature sensor at the riser reactor outlet will ensure the correct estimation for the catalyst circulation rate.However,when uncertain model parameters also exist,additional temperature sensors must be used to ensure correct estimation for unknown inputs and uncertain model parameters of chemical processes.

MODELING AND EXPERIMENTAL STUDY OF A PASSIVE HYDRAULIC ENGINE MOUNT WITH INERTIA AND DIRECT-DECOUPLER

To investigate the dynamic characteristics and damping theory of the passive hydraulic engine mount （PHEM）, numerical prediction is performed through lumped parameter model. System parameters, including volume compliance of the decoupler chamber, effective piston area, fluid inertia and resistance of inertia track and direct-decoupler, are identified by means of experiments and finite element method （FEM）. Dynamic behaviors are tested with elastomer test system for purpose of validating PHEM. With incorporation of inertia track and direct-decoupler, PHEM behaves effective and efficient vibration isolation in range of both low and high frequencies. The comparison of the numerical results with the experimental observations shows that the present PHEM achieves fairly good performance for the engine vibration isolation.

Application of Muskingum routing method with variable parameters in ungauged basin

This paper describes a flood routing method applied in an ungauged basin, utilizing the Muskingum model with variable parameters of wave travel time K and weight coefficient of discharge x based on the physical characteristics of the river reach and flood, including the reach slope, length, width, and flood discharge. Three formulas for estimating parameters of wide rectangular, triangular, and parabolic cross sections are proposed. The influence of the flood on channel flow routing parameters is taken into account. The HEC-HMS hydrological model and the geospatial hydrologic analysis module HEC-GeoHMS were used to extract channel or watershed characteristics and to divide sub-basins. In addition, the initial and constant-rate method, user synthetic unit hydrograph method, and exponential recession method were used to estimate runoff volumes, the direct runoff hydrograph, and the baseflow hydrograph, respectively. The Muskingum model with variable parameters was then applied in the Louzigou Basin in Henan Province of China, and of the results, the percentages of flood events with a relative error of peak discharge less than 20% and runoff volume less than 10% are both 100%. They also show that the percentages of flood events with coefficients of determination greater than 0.8 are 83.33%, 91.67%, and 87.5%, respectively, for rectangular, triangular, and parabolic cross sections in 24 flood events. Therefore, this method is applicable to ungauged basins.