A method for model identification and parameter estimation

International Nuclear Information System (INIS)

Bambach, M; Heinkenschloss, M; Herty, M

2013-01-01

We propose and analyze a new method for the identification of a parameter-dependent model that best describes a given system. This problem arises, for example, in the mathematical modeling of material behavior where several competing constitutive equations are available to describe a given material. In this case, the models are differential equations that arise from the different constitutive equations, and the unknown parameters are coefficients in the constitutive equations. One has to determine the best-suited constitutive equations for a given material and application from experiments. We assume that the true model is one of the N possible parameter-dependent models. To identify the correct model and the corresponding parameters, we can perform experiments, where for each experiment we prescribe an input to the system and observe a part of the system state. Our approach consists of two stages. In the first stage, for each pair of models we determine the experiment, i.e. system input and observation, that best differentiates between the two models, and measure the distance between the two models. Then we conduct N(N − 1) or, depending on the approach taken, N(N − 1)/2 experiments and use the result of the experiments as well as the previously computed model distances to determine the true model. We provide sufficient conditions on the model distances and measurement errors which guarantee that our approach identifies the correct model. Given the model, we identify the corresponding model parameters in the second stage. The problem in the second stage is a standard parameter estimation problem and we use a method suitable for the given application. We illustrate our approach on three examples, including one where the models are elliptic partial differential equations with different parameterized right-hand sides and an example where we identify the constitutive equation in a problem from computational viscoplasticity. (paper)

Automated parameter estimation for biological models using Bayesian statistical model checking.

Science.gov (United States)

Hussain, Faraz; Langmead, Christopher J; Mi, Qi; Dutta-Moscato, Joyeeta; Vodovotz, Yoram; Jha, Sumit K

2015-01-01

Probabilistic models have gained widespread acceptance in the systems biology community as a useful way to represent complex biological systems. Such models are developed using existing knowledge of the structure and dynamics of the system, experimental observations, and inferences drawn from statistical analysis of empirical data. A key bottleneck in building such models is that some system variables cannot be measured experimentally. These variables are incorporated into the model as numerical parameters. Determining values of these parameters that justify existing experiments and provide reliable predictions when model simulations are performed is a key research problem. Using an agent-based model of the dynamics of acute inflammation, we demonstrate a novel parameter estimation algorithm by discovering the amount and schedule of doses of bacterial lipopolysaccharide that guarantee a set of observed clinical outcomes with high probability. We synthesized values of twenty-eight unknown parameters such that the parameterized model instantiated with these parameter values satisfies four specifications describing the dynamic behavior of the model. We have developed a new algorithmic technique for discovering parameters in complex stochastic models of biological systems given behavioral specifications written in a formal mathematical logic. Our algorithm uses Bayesian model checking, sequential hypothesis testing, and stochastic optimization to automatically synthesize parameters of probabilistic biological models.

Seven-parameter statistical model for BRDF in the UV band.

Science.gov (United States)

Bai, Lu; Wu, Zhensen; Zou, Xiren; Cao, Yunhua

2012-05-21

A new semi-empirical seven-parameter BRDF model is developed in the UV band using experimentally measured data. The model is based on the five-parameter model of Wu and the fourteen-parameter model of Renhorn and Boreman. Surface scatter, bulk scatter and retro-reflection scatter are considered. An optimizing modeling method, the artificial immune network genetic algorithm, is used to fit the BRDF measurement data over a wide range of incident angles. The calculation time and accuracy of the five- and seven-parameter models are compared. After fixing the seven parameters, the model can well describe scattering data in the UV band.

Seasonal and spatial variation in broadleaf forest model parameters

Science.gov (United States)

Groenendijk, M.; van der Molen, M. K.; Dolman, A. J.

2009-04-01

Process based, coupled ecosystem carbon, energy and water cycle models are used with the ultimate goal to project the effect of future climate change on the terrestrial carbon cycle. A typical dilemma in such exercises is how much detail the model must be given to describe the observations reasonably realistic while also be general. We use a simple vegetation model (5PM) with five model parameters to study the variability of the parameters. These parameters are derived from the observed carbon and water fluxes from the FLUXNET database. For 15 broadleaf forests the model parameters were derived for different time resolutions. It appears that in general for all forests, the correlation coefficient between observed and simulated carbon and water fluxes improves with a higher parameter time resolution. The quality of the simulations is thus always better when a higher time resolution is used. These results show that annual parameters are not capable of properly describing weather effects on ecosystem fluxes, and that two day time resolution yields the best results. A first indication of the climate constraints can be found by the seasonal variation of the covariance between Jm, which describes the maximum electron transport for photosynthesis, and climate variables. A general seasonality we found is that during winter the covariance with all climate variables is zero. Jm increases rapidly after initial spring warming, resulting in a large covariance with air temperature and global radiation. During summer Jm is less variable, but co-varies negatively with air temperature and vapour pressure deficit and positively with soil water content. A temperature response appears during spring and autumn for broadleaf forests. This shows that an annual model parameter cannot be representative for the entire year. And relations with mean annual temperature are not possible. During summer the photosynthesis parameters are constrained by water availability, soil water content and

Lumped-parameter models

Energy Technology Data Exchange (ETDEWEB)

Ibsen, Lars Bo; Liingaard, M.

2006-12-15

A lumped-parameter model represents the frequency dependent soil-structure interaction of a massless foundation placed on or embedded into an unbounded soil domain. In this technical report the steps of establishing a lumped-parameter model are presented. Following sections are included in this report: Static and dynamic formulation, Simple lumped-parameter models and Advanced lumped-parameter models. (au)

Using the MWC model to describe heterotropic interactions in hemoglobin

Science.gov (United States)

Rapp, Olga

2017-01-01

Hemoglobin is a classical model allosteric protein. Research on hemoglobin parallels the development of key cooperativity and allostery concepts, such as the ‘all-or-none’ Hill formalism, the stepwise Adair binding formulation and the concerted Monod-Wymann-Changuex (MWC) allosteric model. While it is clear that the MWC model adequately describes the cooperative binding of oxygen to hemoglobin, rationalizing the effects of H+, CO2 or organophosphate ligands on hemoglobin-oxygen saturation using the same model remains controversial. According to the MWC model, allosteric ligands exert their effect on protein function by modulating the quaternary conformational transition of the protein. However, data fitting analysis of hemoglobin oxygen saturation curves in the presence or absence of inhibitory ligands persistently revealed effects on both relative oxygen affinity (c) and conformational changes (L), elementary MWC parameters. The recent realization that data fitting analysis using the traditional MWC model equation may not provide reliable estimates for L and c thus calls for a re-examination of previous data using alternative fitting strategies. In the current manuscript, we present two simple strategies for obtaining reliable estimates for MWC mechanistic parameters of hemoglobin steady-state saturation curves in cases of both evolutionary and physiological variations. Our results suggest that the simple MWC model provides a reasonable description that can also account for heterotropic interactions in hemoglobin. The results, moreover, offer a general roadmap for successful data fitting analysis using the MWC model. PMID:28793329

Modelling hydrodynamic parameters to predict flow assisted corrosion

International Nuclear Information System (INIS)

Poulson, B.; Greenwell, B.; Chexal, B.; Horowitz, J.

1992-01-01

During the past 15 years, flow assisted corrosion has been a worldwide problem in the power generating industry. The phenomena is complex and depends on environment, material composition, and hydrodynamic factors. Recently, modeling of flow assisted corrosion has become a subject of great importance. A key part of this effort is modeling the hydrodynamic aspects of this issue. This paper examines which hydrodynamic parameter should be used to correlate the occurrence and rate of flow assisted corrosion with physically meaningful parameters, discusses ways of measuring the relevant hydrodynamic parameter, and describes how the hydrodynamic data is incorporated into the predictive model

Statistical Analysis of Input Parameters Impact on the Modelling of Underground Structures

Directory of Open Access Journals (Sweden)

M. Hilar

2008-01-01

Full Text Available The behaviour of a geomechanical model and its final results are strongly affected by the input parameters. As the inherent variability of rock mass is difficult to model, engineers are frequently forced to face the question “Which input values should be used for analyses?” The correct answer to such a question requires a probabilistic approach, considering the uncertainty of site investigations and variation in the ground. This paper describes the statistical analysis of input parameters for FEM calculations of traffic tunnels in the city of Prague. At the beginning of the paper, the inaccuracy in the geotechnical modelling is discussed. In the following part the Fuzzy techniques are summarized, including information about an application of the Fuzzy arithmetic on the shotcrete parameters. The next part of the paper is focused on the stochastic simulation – Monte Carlo Simulation is briefly described, Latin Hypercubes method is described more in details. At the end several practical examples are described: statistical analysis of the input parameters on the numerical modelling of the completed Mrázovka tunnel (profile West Tunnel Tube km 5.160 and modelling of the constructed tunnel Špejchar – Pelc Tyrolka. 

Determining Rheological Parameters of Generalized Yield-Power-Law Fluid Model

Directory of Open Access Journals (Sweden)

Stryczek Stanislaw

2004-09-01

Full Text Available The principles of determining rheological parameters of drilling muds described by a generalized yield-power-law are presented in the paper. Functions between tangent stresses and shear rate are given. The conditions of laboratory measurements of rheological parameters of generalized yield-power-law fluids are described and necessary mathematical relations for rheological model parameters given. With the block diagrams, the methodics of numerical solution of these relations has been presented. Rheological parameters of an exemplary drilling mud have been calculated with the use of this numerical program.

Microarchitecture Parameters Describe Bone Structure and Its Strength Better Than BMD

Directory of Open Access Journals (Sweden)

Tomasz Topoliński

2012-01-01

Full Text Available Introduction and Hypothesis. Some papers have shown that bone mineral density (BMD may not be accurate in predicting fracture risk. Recently microarchitecture parameters have been reported to give information on bone characteristics. The aim of this study was to find out if the values of volume, fractal dimension, and bone mineral density are correlated with bone strength. Methods. Forty-two human bone samples harvested during total hip replacement surgery were cut to cylindrical samples. The geometrical mesh of layers of bone mass obtained from microCT investigation and the volumes of each layer and fractal dimension were calculated. The finite element method was applied to calculate the compression force F causing ε=0.8% strain. Results. There were stronger correlations for microarchitecture parameters with strength than those for bone mineral density. The values of determination coefficient R2 for mean volume and force were 0.88 and 0.90 for mean fractal dimension and force, while for BMD and force the value was 0.53. The samples with bigger mean bone volume of layers and bigger mean fractal dimension of layers (more complex structure presented higher strength. Conclusion. The volumetric and fractal dimension parameters better describe bone structure and strength than BMD.

Environmental Transport Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

M. A. Wasiolek

2003-06-27

This analysis report is one of the technical reports documenting the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN), a biosphere model supporting the total system performance assessment (TSPA) for the geologic repository at Yucca Mountain. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows relationships among the reports developed for biosphere modeling and biosphere abstraction products for the TSPA, as identified in the ''Technical Work Plan: for Biosphere Modeling and Expert Support'' (TWP) (BSC 2003 [163602]). Some documents in Figure 1-1 may be under development and not available when this report is issued. This figure provides an understanding of how this report contributes to biosphere modeling in support of the license application (LA), but access to the listed documents is not required to understand the contents of this report. This report is one of the reports that develops input parameter values for the biosphere model. The ''Biosphere Model Report'' (BSC 2003 [160699]) describes the conceptual model, the mathematical model, and the input parameters. The purpose of this analysis is to develop biosphere model parameter values related to radionuclide transport and accumulation in the environment. These parameters support calculations of radionuclide concentrations in the environmental media (e.g., soil, crops, animal products, and air) resulting from a given radionuclide concentration at the source of contamination (i.e., either in groundwater or volcanic ash). The analysis was performed in accordance with the TWP (BSC 2003 [163602]). This analysis develops values of parameters associated with many features, events, and processes (FEPs) applicable to the reference biosphere (DTN: M00303SEPFEPS2.000 [162452]), which are addressed in the biosphere model (BSC 2003 [160699]). The treatment of these FEPs is described in BSC (2003 [160699

Environmental Transport Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

Wasiolek, M. A.

2003-01-01

This analysis report is one of the technical reports documenting the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN), a biosphere model supporting the total system performance assessment (TSPA) for the geologic repository at Yucca Mountain. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows relationships among the reports developed for biosphere modeling and biosphere abstraction products for the TSPA, as identified in the ''Technical Work Plan: for Biosphere Modeling and Expert Support'' (TWP) (BSC 2003 [163602]). Some documents in Figure 1-1 may be under development and not available when this report is issued. This figure provides an understanding of how this report contributes to biosphere modeling in support of the license application (LA), but access to the listed documents is not required to understand the contents of this report. This report is one of the reports that develops input parameter values for the biosphere model. The ''Biosphere Model Report'' (BSC 2003 [160699]) describes the conceptual model, the mathematical model, and the input parameters. The purpose of this analysis is to develop biosphere model parameter values related to radionuclide transport and accumulation in the environment. These parameters support calculations of radionuclide concentrations in the environmental media (e.g., soil, crops, animal products, and air) resulting from a given radionuclide concentration at the source of contamination (i.e., either in groundwater or volcanic ash). The analysis was performed in accordance with the TWP (BSC 2003 [163602]). This analysis develops values of parameters associated with many features, events, and processes (FEPs) applicable to the reference biosphere (DTN: M00303SEPFEPS2.000 [162452]), which are addressed in the biosphere model (BSC 2003 [160699]). The treatment of these FEPs is described in BSC (2003 [160699], Section 6.2). Parameter values

Distribution-centric 3-parameter thermodynamic models of partition gas chromatography.

Science.gov (United States)

Blumberg, Leonid M

2017-03-31

If both parameters (the entropy, ΔS, and the enthalpy, ΔH) of the classic van't Hoff model of dependence of distribution coefficients (K) of analytes on temperature (T) are treated as the temperature-independent constants then the accuracy of the model is known to be insufficient for the needed accuracy of retention time prediction. A more accurate 3-parameter Clarke-Glew model offers a way to treat ΔS and ΔH as functions, ΔS(T) and ΔH(T), of T. A known T-centric construction of these functions is based on relating them to the reference values (ΔS ref and ΔH ref ) corresponding to a predetermined reference temperature (T ref ). Choosing a single T ref for all analytes in a complex sample or in a large database might lead to practically irrelevant values of ΔS ref and ΔH ref for those analytes that have too small or too large retention factors at T ref . Breaking all analytes in several subsets each with its own T ref leads to discontinuities in the analyte parameters. These problems are avoided in the K-centric modeling where ΔS(T) and ΔS(T) and other analyte parameters are described in relation to their values corresponding to a predetermined reference distribution coefficient (K Ref ) - the same for all analytes. In this report, the mathematics of the K-centric modeling are described and the properties of several types of K-centric parameters are discussed. It has been shown that the earlier introduced characteristic parameters of the analyte-column interaction (the characteristic temperature, T char , and the characteristic thermal constant, θ char ) are a special chromatographically convenient case of the K-centric parameters. Transformations of T-centric parameters into K-centric ones and vice-versa as well as the transformations of one set of K-centric parameters into another set and vice-versa are described. Copyright © 2017 Elsevier B.V. All rights reserved.

Constant-parameter capture-recapture models

Science.gov (United States)

Brownie, C.; Hines, J.E.; Nichols, J.D.

1986-01-01

Jolly (1982, Biometrics 38, 301-321) presented modifications of the Jolly-Seber model for capture-recapture data, which assume constant survival and/or capture rates. Where appropriate, because of the reduced number of parameters, these models lead to more efficient estimators than the Jolly-Seber model. The tests to compare models given by Jolly do not make complete use of the data, and we present here the appropriate modifications, and also indicate how to carry out goodness-of-fit tests which utilize individual capture history information. We also describe analogous models for the case where young and adult animals are tagged. The availability of computer programs to perform the analysis is noted, and examples are given using output from these programs.

Radon decay product in-door behaviour - parameter, measurement method, and model review

International Nuclear Information System (INIS)

Scofield, P.

1988-01-01

This report reviews parameters used to characterize indoor radon daughter behavior and concentrations. Certain parameters that affect indoor radon daughter concentrations are described and the values obtained experimentally or theoretically are summarized. Radon daughter measurement methods are reviewed, such as, PAEC, unattached daughters, particle size distributions, and plateout measurement methods. In addition, certain radon pressure driven/diffusion models and indoor radon daughter models are briefly described. (orig.)

Describing Growth Pattern of Bali Cows Using Non-linear Regression Models

Directory of Open Access Journals (Sweden)

Mohd. Hafiz A.W

2016-12-01

Full Text Available The objective of this study was to evaluate the best fit non-linear regression model to describe the growth pattern of Bali cows. Estimates of asymptotic mature weight, rate of maturing and constant of integration were derived from Brody, von Bertalanffy, Gompertz and Logistic models which were fitted to cross-sectional data of body weight taken from 74 Bali cows raised in MARDI Research Station Muadzam Shah Pahang. Coefficient of determination (R2 and residual mean squares (MSE were used to determine the best fit model in describing the growth pattern of Bali cows. Von Bertalanffy model was the best model among the four growth functions evaluated to determine the mature weight of Bali cattle as shown by the highest R2 and lowest MSE values (0.973 and 601.9, respectively, followed by Gompertz (0.972 and 621.2, respectively, Logistic (0.971 and 648.4, respectively and Brody (0.932 and 660.5, respectively models. The correlation between rate of maturing and mature weight was found to be negative in the range of -0.170 to -0.929 for all models, indicating that animals of heavier mature weight had lower rate of maturing. The use of non-linear model could summarize the weight-age relationship into several biologically interpreted parameters compared to the entire lifespan weight-age data points that are difficult and time consuming to interpret.

Zener Diode Compact Model Parameter Extraction Using Xyce-Dakota Optimization.

Energy Technology Data Exchange (ETDEWEB)

Buchheit, Thomas E. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Wilcox, Ian Zachary [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Sandoval, Andrew J [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Reza, Shahed [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-12-01

This report presents a detailed process for compact model parameter extraction for DC circuit Zener diodes. Following the traditional approach of Zener diode parameter extraction, circuit model representation is defined and then used to capture the different operational regions of a real diode's electrical behavior. The circuit model contains 9 parameters represented by resistors and characteristic diodes as circuit model elements. The process of initial parameter extraction, the identification of parameter values for the circuit model elements, is presented in a way that isolates the dependencies between certain electrical parameters and highlights both the empirical nature of the extraction and portions of the real diode physical behavior which of the parameters are intended to represent. Optimization of the parameters, a necessary part of a robost parameter extraction process, is demonstrated using a 'Xyce-Dakota' workflow, discussed in more detail in the report. Among other realizations during this systematic approach of electrical model parameter extraction, non-physical solutions are possible and can be difficult to avoid because of the interdependencies between the different parameters. The process steps described are fairly general and can be leveraged for other types of semiconductor device model extractions. Also included in the report are recommendations for experiment setups for generating optimum dataset for model extraction and the Parameter Identification and Ranking Table (PIRT) for Zener diodes.

Describing pediatric dysphonia with nonlinear dynamic parameters

Science.gov (United States)

Meredith, Morgan L.; Theis, Shannon M.; McMurray, J. Scott; Zhang, Yu; Jiang, Jack J.

2008-01-01

Objective Nonlinear dynamic analysis has emerged as a reliable and objective tool for assessing voice disorders. However, it has only been tested on adult populations. In the present study, nonlinear dynamic analysis was applied to normal and dysphonic pediatric populations with the goal of collecting normative data. Jitter analysis was also applied in order to compare nonlinear dynamic and perturbation measures. This study’s findings will be useful in creating standards for the use of nonlinear dynamic analysis as a tool to describe dysphonia in the pediatric population. Methods The study included 38 pediatric subjects (23 children with dysphonia and 15 without). Recordings of sustained vowels were obtained from each subject and underwent nonlinear dynamic analysis and percent jitter analysis. The resulting correlation dimension (D2) and percent jitter values were compared across the two groups using t-tests set at a significance level of p = 0.05. Results It was shown that D2 values covary with the presence of pathology in children. D2 values were significantly higher in dysphonic children than in normal children (p = 0.002). Standard deviations indicated a higher level of variation in normal children’s D2 values than in dysphonic children’s D2 values. Jitter analysis showed markedly higher percent jitter in dysphonic children than in normal children (p = 0.025) and large standard deviations for both groups. Conclusion This study indicates that nonlinear dynamic analysis could be a viable tool for the detection and assessment of dysphonia in children. Further investigations and more normative data are needed to create standards for using nonlinear dynamic parameters for the clinical evaluation of pediatric dysphonia. PMID:18947887

Applicability of Different Hydraulic Parameters to Describe Soil Detachment in Eroding Rills

Science.gov (United States)

Wirtz, Stefan; Seeger, Manuel; Zell, Andreas; Wagner, Christian; Wagner, Jean-Frank; Ries, Johannes B.

2013-01-01

This study presents the comparison of experimental results with assumptions used in numerical models. The aim of the field experiments is to test the linear relationship between different hydraulic parameters and soil detachment. For example correlations between shear stress, unit length shear force, stream power, unit stream power and effective stream power and the detachment rate does not reveal a single parameter which consistently displays the best correlation. More importantly, the best fit does not only vary from one experiment to another, but even between distinct measurement points. Different processes in rill erosion are responsible for the changing correlations. However, not all these procedures are considered in soil erosion models. Hence, hydraulic parameters alone are not sufficient to predict detachment rates. They predict the fluvial incising in the rill's bottom, but the main sediment sources are not considered sufficiently in its equations. The results of this study show that there is still a lack of understanding of the physical processes underlying soil erosion. Exerted forces, soil stability and its expression, the abstraction of the detachment and transport processes in shallow flowing water remain still subject of unclear description and dependence. PMID:23717669

Agricultural and Environmental Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

Kaylie Rasmuson; Kurt Rautenstrauch

2003-01-01

This analysis is one of nine technical reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) biosphere model. It documents input parameters for the biosphere model, and supports the use of the model to develop Biosphere Dose Conversion Factors (BDCF). The biosphere model is one of a series of process models supporting the Total System Performance Assessment (TSPA) for the repository at Yucca Mountain. The ERMYN provides the TSPA with the capability to perform dose assessments. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships between the major activities and their products (the analysis and model reports) that were planned in the biosphere Technical Work Plan (TWP, BSC 2003a). It should be noted that some documents identified in Figure 1-1 may be under development and therefore not available at the time this document is issued. The ''Biosphere Model Report'' (BSC 2003b) describes the ERMYN and its input parameters. This analysis report, ANL-MGR-MD-000006, ''Agricultural and Environmental Input Parameters for the Biosphere Model'', is one of the five reports that develop input parameters for the biosphere model. This report defines and justifies values for twelve parameters required in the biosphere model. These parameters are related to use of contaminated groundwater to grow crops. The parameter values recommended in this report are used in the soil, plant, and carbon-14 submodels of the ERMYN

Lumped-Parameter Models for Windturbine Footings on Layered Ground

DEFF Research Database (Denmark)

Andersen, Lars

The design of modern wind turbines is typically based on lifetime analyses using aeroelastic codes. In this regard, the impedance of the foundations must be described accurately without increasing the overall size of the computationalmodel significantly. This may be obtained by the fitting...... of a lumped-parameter model to the results of a rigorous model or experimental results. In this paper, guidelines are given for the formulation of such lumped-parameter models and examples are given in which the models are utilised for the analysis of a wind turbine supported by a surface footing on a layered...

Cosmological model with viscosity media (dark fluid) described by an effective equation of state

International Nuclear Information System (INIS)

Ren Jie; Meng Xinhe

2006-01-01

A generally parameterized equation of state (EOS) is investigated in the cosmological evolution with bulk viscosity media modelled as dark fluid, which can be regarded as a unification of dark energy and dark matter. Compared with the case of the perfect fluid, this EOS has possessed four additional parameters, which can be interpreted as the case of the non-perfect fluid with time-dependent viscosity or the model with variable cosmological constant. From this general EOS, a completely integrable dynamical equation to the scale factor is obtained with its solution explicitly given out. (i) In this parameterized model of cosmology, for a special choice of the parameters we can explain the late-time accelerating expansion universe in a new view. The early inflation, the median (relatively late time) deceleration, and the recently cosmic acceleration may be unified in a single equation. (ii) A generalized relation of the Hubble parameter scaling with the redshift is obtained for some cosmology interests. (iii) By using the SNe Ia data to fit the effective viscosity model we show that the case of matter described by p=0 plus with effective viscosity contributions can fit the observational gold data in an acceptable level

Suitability of parametric models to describe the hydraulic properties of an unsaturated coarse sand and gravel

Science.gov (United States)

Mace, Andy; Rudolph, David L.; Kachanoski , R. Gary

1998-01-01

The performance of parametric models used to describe soil water retention (SWR) properties and predict unsaturated hydraulic conductivity (K) as a function of volumetric water content (θ) is examined using SWR and K(θ) data for coarse sand and gravel sediments. Six 70 cm long, 10 cm diameter cores of glacial outwash were instrumented at eight depths with porous cup ten-siometers and time domain reflectometry probes to measure soil water pressure head (h) and θ, respectively, for seven unsaturated and one saturated steady-state flow conditions. Forty-two θ(h) and K(θ) relationships were measured from the infiltration tests on the cores. Of the four SWR models compared in the analysis, the van Genuchten (1980) equation with parameters m and n restricted according to the Mualem (m = 1 - 1/n) criterion is best suited to describe the θ(h) relationships. The accuracy of two models that predict K(θ) using parameter values derived from the SWR models was also evaluated. The model developed by van Genuchten (1980) based on the theoretical expression of Mualem (1976) predicted K(θ) more accurately than the van Genuchten (1980) model based on the theory of Burdine (1953). A sensitivity analysis shows that more accurate predictions of K(θ) are achieved using SWR model parameters derived with residual water content (θr) specified according to independent measurements of θ at values of h where θ/h ∼ 0 rather than model-fit θr values. The accuracy of the model K(θ) function improves markedly when at least one value of unsaturated K is used to scale the K(θ) function predicted using the saturated K. The results of this investigation indicate that the hydraulic properties of coarse-grained sediments can be accurately described using the parametric models. In addition, data collection efforts should focus on measuring at least one value of unsaturated hydraulic conductivity and as complete a set of SWR data as possible, particularly in the dry range.

Importance of hydrological parameters in contaminant transport modeling in a terrestrial environment

International Nuclear Information System (INIS)

Tsuduki, Katsunori; Matsunaga, Takeshi

2007-01-01

A grid type multi-layered distributed parameter model for calculating discharge in a watershed was described. Model verification with our field observation resulted in different sets of hydrological parameter values, all of which reproduced the observed discharge. The effect of those varied hydrological parameters on contaminant transport calculation was examined and discussed by simulation of event water transfer. (author)

An improved robust model predictive control for linear parameter-varying input-output models

NARCIS (Netherlands)

Abbas, H.S.; Hanema, J.; Tóth, R.; Mohammadpour, J.; Meskin, N.

2018-01-01

This paper describes a new robust model predictive control (MPC) scheme to control the discrete-time linear parameter-varying input-output models subject to input and output constraints. Closed-loop asymptotic stability is guaranteed by including a quadratic terminal cost and an ellipsoidal terminal

An extended car-following model to describe connected traffic dynamics under cyberattacks

Science.gov (United States)

Wang, Pengcheng; Yu, Guizhen; Wu, Xinkai; Qin, Hongmao; Wang, Yunpeng

2018-04-01

In this paper, the impacts of the potential cyberattacks on vehicles are modeled through an extended car-following model. To better understand the mechanism of traffic disturbance under cyberattacks, the linear and nonlinear stability analysis are conducted respectively. Particularly, linear stability analysis is performed to obtain different neutral stability conditions with various parameters; and nonlinear stability analysis is carried out by using reductive perturbation method to derive the soliton solution of the modified Korteweg de Vries equation (mKdV) near the critical point, which is used to draw coexisting stability lines. Furthermore, by applying linear and nonlinear stability analysis, traffic flow state can be divided into three states, i.e., stable, metastable and unstable states which are useful to describe shockwave dynamics and driving behaviors under cyberattacks. The theoretical results show that the proposed car-following model is capable of successfully describing the car-following behavior of connected vehicles with cyberattacks. Finally, numerical simulation using real values has confirmed the validity of theoretical analysis. The results further demonstrate our model can be used to help avoid collisions and relieve traffic congestion with cybersecurity threats.

Agricultural and Environmental Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

Kaylie Rasmuson; Kurt Rautenstrauch

2003-06-20

This analysis is one of nine technical reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) biosphere model. It documents input parameters for the biosphere model, and supports the use of the model to develop Biosphere Dose Conversion Factors (BDCF). The biosphere model is one of a series of process models supporting the Total System Performance Assessment (TSPA) for the repository at Yucca Mountain. The ERMYN provides the TSPA with the capability to perform dose assessments. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships between the major activities and their products (the analysis and model reports) that were planned in the biosphere Technical Work Plan (TWP, BSC 2003a). It should be noted that some documents identified in Figure 1-1 may be under development and therefore not available at the time this document is issued. The ''Biosphere Model Report'' (BSC 2003b) describes the ERMYN and its input parameters. This analysis report, ANL-MGR-MD-000006, ''Agricultural and Environmental Input Parameters for the Biosphere Model'', is one of the five reports that develop input parameters for the biosphere model. This report defines and justifies values for twelve parameters required in the biosphere model. These parameters are related to use of contaminated groundwater to grow crops. The parameter values recommended in this report are used in the soil, plant, and carbon-14 submodels of the ERMYN.

Cognitive Models of Risky Choice: Parameter Stability and Predictive Accuracy of Prospect Theory

Science.gov (United States)

Glockner, Andreas; Pachur, Thorsten

2012-01-01

In the behavioral sciences, a popular approach to describe and predict behavior is cognitive modeling with adjustable parameters (i.e., which can be fitted to data). Modeling with adjustable parameters allows, among other things, measuring differences between people. At the same time, parameter estimation also bears the risk of overfitting. Are…

Retrospective forecast of ETAS model with daily parameters estimate

Science.gov (United States)

Falcone, Giuseppe; Murru, Maura; Console, Rodolfo; Marzocchi, Warner; Zhuang, Jiancang

2016-04-01

We present a retrospective ETAS (Epidemic Type of Aftershock Sequence) model based on the daily updating of free parameters during the background, the learning and the test phase of a seismic sequence. The idea was born after the 2011 Tohoku-Oki earthquake. The CSEP (Collaboratory for the Study of Earthquake Predictability) Center in Japan provided an appropriate testing benchmark for the five 1-day submitted models. Of all the models, only one was able to successfully predict the number of events that really happened. This result was verified using both the real time and the revised catalogs. The main cause of the failure was in the underestimation of the forecasted events, due to model parameters maintained fixed during the test. Moreover, the absence in the learning catalog of an event similar to the magnitude of the mainshock (M9.0), which drastically changed the seismicity in the area, made the learning parameters not suitable to describe the real seismicity. As an example of this methodological development we show the evolution of the model parameters during the last two strong seismic sequences in Italy: the 2009 L'Aquila and the 2012 Reggio Emilia episodes. The achievement of the model with daily updated parameters is compared with that of same model where the parameters remain fixed during the test time.

Lumped-Parameter Models for Wind-Turbine Footings on Layered Ground

DEFF Research Database (Denmark)

Andersen, Lars; Liingaard, Morten

2007-01-01

The design of modern wind turbines is typically based on lifetime analyses using aeroelastic codes. In this regard, the impedance of the foundations must be described accurately without increasing the overall size of the computational model significantly. This may be obtained by the fitting...... of a lumped-parameter model to the results of a rigorous model or experimental results. In this paper, guidelines are given for the formulation of such lumped-parameter models and examples are given in which the models are utilised for the analysis of a wind turbine supported by a surface footing on a layered...

Evidence in Support of the Independent Channel Model Describing the Sensorimotor Control of Human Stance Using a Humanoid Robot.

Science.gov (United States)

Pasma, Jantsje H; Assländer, Lorenz; van Kordelaar, Joost; de Kam, Digna; Mergner, Thomas; Schouten, Alfred C

2018-01-01

The Independent Channel (IC) model is a commonly used linear balance control model in the frequency domain to analyze human balance control using system identification and parameter estimation. The IC model is a rudimentary and noise-free description of balance behavior in the frequency domain, where a stable model representation is not guaranteed. In this study, we conducted firstly time-domain simulations with added noise, and secondly robot experiments by implementing the IC model in a real-world robot (PostuRob II) to test the validity and stability of the model in the time domain and for real world situations. Balance behavior of seven healthy participants was measured during upright stance by applying pseudorandom continuous support surface rotations. System identification and parameter estimation were used to describe the balance behavior with the IC model in the frequency domain. The IC model with the estimated parameters from human experiments was implemented in Simulink for computer simulations including noise in the time domain and robot experiments using the humanoid robot PostuRob II. Again, system identification and parameter estimation were used to describe the simulated balance behavior. Time series, Frequency Response Functions, and estimated parameters from human experiments, computer simulations, and robot experiments were compared with each other. The computer simulations showed similar balance behavior and estimated control parameters compared to the human experiments, in the time and frequency domain. Also, the IC model was able to control the humanoid robot by keeping it upright, but showed small differences compared to the human experiments in the time and frequency domain, especially at high frequencies. We conclude that the IC model, a descriptive model in the frequency domain, can imitate human balance behavior also in the time domain, both in computer simulations with added noise and real world situations with a humanoid robot. This

Evidence in Support of the Independent Channel Model Describing the Sensorimotor Control of Human Stance Using a Humanoid Robot

Directory of Open Access Journals (Sweden)

Jantsje H. Pasma

2018-03-01

Full Text Available The Independent Channel (IC model is a commonly used linear balance control model in the frequency domain to analyze human balance control using system identification and parameter estimation. The IC model is a rudimentary and noise-free description of balance behavior in the frequency domain, where a stable model representation is not guaranteed. In this study, we conducted firstly time-domain simulations with added noise, and secondly robot experiments by implementing the IC model in a real-world robot (PostuRob II to test the validity and stability of the model in the time domain and for real world situations. Balance behavior of seven healthy participants was measured during upright stance by applying pseudorandom continuous support surface rotations. System identification and parameter estimation were used to describe the balance behavior with the IC model in the frequency domain. The IC model with the estimated parameters from human experiments was implemented in Simulink for computer simulations including noise in the time domain and robot experiments using the humanoid robot PostuRob II. Again, system identification and parameter estimation were used to describe the simulated balance behavior. Time series, Frequency Response Functions, and estimated parameters from human experiments, computer simulations, and robot experiments were compared with each other. The computer simulations showed similar balance behavior and estimated control parameters compared to the human experiments, in the time and frequency domain. Also, the IC model was able to control the humanoid robot by keeping it upright, but showed small differences compared to the human experiments in the time and frequency domain, especially at high frequencies. We conclude that the IC model, a descriptive model in the frequency domain, can imitate human balance behavior also in the time domain, both in computer simulations with added noise and real world situations with a

Environmental Transport Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

M. Wasiolek

2004-01-01

This analysis report is one of the technical reports documenting the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN), a biosphere model supporting the total system performance assessment for the license application (TSPA-LA) for the geologic repository at Yucca Mountain. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows relationships among the reports developed for biosphere modeling and biosphere abstraction products for the TSPA-LA, as identified in the ''Technical Work Plan for Biosphere Modeling and Expert Support'' (BSC 2004 [DIRS 169573]) (TWP). This figure provides an understanding of how this report contributes to biosphere modeling in support of the license application (LA). This report is one of the five reports that develop input parameter values for the biosphere model. The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the conceptual model and the mathematical model. The input parameter reports, shown to the right of the Biosphere Model Report in Figure 1-1, contain detailed description of the model input parameters. The output of this report is used as direct input in the ''Nominal Performance Biosphere Dose Conversion Factor Analysis'' and in the ''Disruptive Event Biosphere Dose Conversion Factor Analysis'' that calculate the values of biosphere dose conversion factors (BDCFs) for the groundwater and volcanic ash exposure scenarios, respectively. The purpose of this analysis was to develop biosphere model parameter values related to radionuclide transport and accumulation in the environment. These parameters support calculations of radionuclide concentrations in the environmental media (e.g., soil, crops, animal products, and air) resulting from a given radionuclide concentration at the source of contamination (i.e., either in groundwater or in volcanic ash). The analysis was performed in accordance with the TWP (BSC 2004 [DIRS 169573])

Environmental Transport Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

M. Wasiolek

2004-09-10

This analysis report is one of the technical reports documenting the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN), a biosphere model supporting the total system performance assessment for the license application (TSPA-LA) for the geologic repository at Yucca Mountain. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows relationships among the reports developed for biosphere modeling and biosphere abstraction products for the TSPA-LA, as identified in the ''Technical Work Plan for Biosphere Modeling and Expert Support'' (BSC 2004 [DIRS 169573]) (TWP). This figure provides an understanding of how this report contributes to biosphere modeling in support of the license application (LA). This report is one of the five reports that develop input parameter values for the biosphere model. The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the conceptual model and the mathematical model. The input parameter reports, shown to the right of the Biosphere Model Report in Figure 1-1, contain detailed description of the model input parameters. The output of this report is used as direct input in the ''Nominal Performance Biosphere Dose Conversion Factor Analysis'' and in the ''Disruptive Event Biosphere Dose Conversion Factor Analysis'' that calculate the values of biosphere dose conversion factors (BDCFs) for the groundwater and volcanic ash exposure scenarios, respectively. The purpose of this analysis was to develop biosphere model parameter values related to radionuclide transport and accumulation in the environment. These parameters support calculations of radionuclide concentrations in the environmental media (e.g., soil, crops, animal products, and air) resulting from a given radionuclide concentration at the source of contamination (i.e., either in groundwater or in volcanic ash). The analysis

Four-parameter model for polarization-resolved rough-surface BRDF.

Science.gov (United States)

Renhorn, Ingmar G E; Hallberg, Tomas; Bergström, David; Boreman, Glenn D

2011-01-17

A modeling procedure is demonstrated, which allows representation of polarization-resolved BRDF data using only four parameters: the real and imaginary parts of an effective refractive index with an added parameter taking grazing incidence absorption into account and an angular-scattering parameter determined from the BRDF measurement of a chosen angle of incidence, preferably close to normal incidence. These parameters allow accurate predictions of s- and p-polarized BRDF for a painted rough surface, over three decades of variation in BRDF magnitude. To characterize any particular surface of interest, the measurements required to determine these four parameters are the directional hemispherical reflectance (DHR) for s- and p-polarized input radiation and the BRDF at a selected angle of incidence. The DHR data describes the angular and polarization dependence, as well as providing the overall normalization constraint. The resulting model conserves energy and fulfills the reciprocity criteria.

Soil-Related Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

Smith, A. J.

2004-01-01

This report presents one of the analyses that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN). The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the details of the conceptual model as well as the mathematical model and the required input parameters. The biosphere model is one of a series of process models supporting the postclosure Total System Performance Assessment (TSPA) for the Yucca Mountain repository. A schematic representation of the documentation flow for the Biosphere input to TSPA is presented in Figure 1-1. This figure shows the evolutionary relationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the biosphere abstraction products for TSPA, as identified in the ''Technical Work Plan for Biosphere Modeling and Expert Support'' (TWP) (BSC 2004 [DIRS 169573]). This figure is included to provide an understanding of how this analysis report contributes to biosphere modeling in support of the license application, and is not intended to imply that access to the listed documents is required to understand the contents of this report. This report, ''Soil-Related Input Parameters for the Biosphere Model'', is one of the five analysis reports that develop input parameters for use in the ERMYN model. This report is the source documentation for the six biosphere parameters identified in Table 1-1. The purpose of this analysis was to develop the biosphere model parameters associated with the accumulation and depletion of radionuclides in the soil. These parameters support the calculation of radionuclide concentrations in soil from on-going irrigation or ash deposition and, as a direct consequence, radionuclide concentration in other environmental media that are affected by radionuclide concentrations in soil. The analysis was performed in accordance with the TWP (BSC 2004 [DIRS 169573]) where the governing procedure was defined as AP-SIII.9Q, ''Scientific Analyses''. This

Metamodel-based inverse method for parameter identification: elastic-plastic damage model

Science.gov (United States)

Huang, Changwu; El Hami, Abdelkhalak; Radi, Bouchaïb

2017-04-01

This article proposed a metamodel-based inverse method for material parameter identification and applies it to elastic-plastic damage model parameter identification. An elastic-plastic damage model is presented and implemented in numerical simulation. The metamodel-based inverse method is proposed in order to overcome the disadvantage in computational cost of the inverse method. In the metamodel-based inverse method, a Kriging metamodel is constructed based on the experimental design in order to model the relationship between material parameters and the objective function values in the inverse problem, and then the optimization procedure is executed by the use of a metamodel. The applications of the presented material model and proposed parameter identification method in the standard A 2017-T4 tensile test prove that the presented elastic-plastic damage model is adequate to describe the material's mechanical behaviour and that the proposed metamodel-based inverse method not only enhances the efficiency of parameter identification but also gives reliable results.

Applications of the solvation parameter model in reversed-phase liquid chromatography.

Science.gov (United States)

Poole, Colin F; Lenca, Nicole

2017-02-24

The solvation parameter model is widely used to provide insight into the retention mechanism in reversed-phase liquid chromatography, for column characterization, and in the development of surrogate chromatographic models for biopartitioning processes. The properties of the separation system are described by five system constants representing all possible intermolecular interactions for neutral molecules. The general model can be extended to include ions and enantiomers by adding new descriptors to encode the specific properties of these compounds. System maps provide a comprehensive overview of the separation system as a function of mobile phase composition and/or temperature for method development. The solvation parameter model has been applied to gradient elution separations but here theory and practice suggest a cautious approach since the interpretation of system and compound properties derived from its use are approximate. A growing application of the solvation parameter model in reversed-phase liquid chromatography is the screening of surrogate chromatographic systems for estimating biopartitioning properties. Throughout the discussion of the above topics success as well as known and likely deficiencies of the solvation parameter model are described with an emphasis on the role of the heterogeneous properties of the interphase region on the interpretation and understanding of the general retention mechanism in reversed-phase liquid chromatography for porous chemically bonded sorbents. Copyright © 2016 Elsevier B.V. All rights reserved.

Description of the Hexadecapole Deformation Parameter in the sdg Interacting Boson Model

Science.gov (United States)

Liu, Yu-xin; Sun, Di; Wang, Jia-jun; Han, Qi-zhi

1998-04-01

The hexadecapole deformation parameter β4 of the rare-earth and actinide nuclei is investigated in the framework of the sdg interacing boson model. An explicit relation between the geometric hexadecapole deformation parameter β4 and the intrinsic deformation parameters epsilon4, epsilon2 are obtained. The deformation parameters β4 of the rare-earths and actinides are determined without any free parameter. The calculated results agree with experimental data well. It also shows that the SU(5) limit of the sdg interacting boson model can describe the β4 systematics as well as the SU(3) limit.

Description of the hexadecapole deformation parameter in the sdg interacting boson model

International Nuclear Information System (INIS)

Liu Yuxin; Sun Di; Wang Jiajun; Han Qizhi

1998-01-01

The hexadecapole deformation parameter β 4 of the rare-earth and actinide nuclei is investigated in the framework of the sdg interacting boson model. An explicit relation between the geometric hexadecapole deformation parameter β 4 and the intrinsic deformation parameters ε 4 , ε 2 are obtained. The deformation parameters β 4 of the rare-earths and actinides are determined without any free parameter. The calculated results agree with experimental data well. It also shows that the SU(5) limit of the sdg interacting boson model can describe the β 4 systematics as well as the SU(3) limit

Model checking biological systems described using ambient calculus

DEFF Research Database (Denmark)

Mardare, Radu Iulian; Priami, Corrado; Qualia, Paola

2005-01-01

Model checking biological systems described using ambient calculus. In Proc. of the second International Workshop on Computational Methods in Systems Biology (CMSB04), Lecture Notes in Bioinformatics 3082:85-103, Springer, 2005.......Model checking biological systems described using ambient calculus. In Proc. of the second International Workshop on Computational Methods in Systems Biology (CMSB04), Lecture Notes in Bioinformatics 3082:85-103, Springer, 2005....

q-deformed Einstein's model to describe specific heat of solid

Science.gov (United States)

Guha, Atanu; Das, Prasanta Kumar

2018-04-01

Realistic phenomena can be described more appropriately using generalized canonical ensemble, with proper parameter sets involved. We have generalized the Einstein's theory for specific heat of solid in Tsallis statistics, where the temperature fluctuation is introduced into the theory via the fluctuation parameter q. At low temperature the Einstein's curve of the specific heat in the nonextensive Tsallis scenario exactly lies on the experimental data points. Consequently this q-modified Einstein's curve is found to be overlapping with the one predicted by Debye. Considering only the temperature fluctuation effect(even without considering more than one mode of vibration is being triggered) we found that the CV vs T curve is as good as obtained by considering the different modes of vibration as suggested by Debye. Generalizing the Einstein's theory in Tsallis statistics we found that a unique value of the Einstein temperature θE along with a temperature dependent deformation parameter q(T) , can well describe the phenomena of specific heat of solid i.e. the theory is equivalent to Debye's theory with a temperature dependent θD.

Estimating Parameters in Physical Models through Bayesian Inversion: A Complete Example

KAUST Repository

Allmaras, Moritz

2013-02-07

All mathematical models of real-world phenomena contain parameters that need to be estimated from measurements, either for realistic predictions or simply to understand the characteristics of the model. Bayesian statistics provides a framework for parameter estimation in which uncertainties about models and measurements are translated into uncertainties in estimates of parameters. This paper provides a simple, step-by-step example-starting from a physical experiment and going through all of the mathematics-to explain the use of Bayesian techniques for estimating the coefficients of gravity and air friction in the equations describing a falling body. In the experiment we dropped an object from a known height and recorded the free fall using a video camera. The video recording was analyzed frame by frame to obtain the distance the body had fallen as a function of time, including measures of uncertainty in our data that we describe as probability densities. We explain the decisions behind the various choices of probability distributions and relate them to observed phenomena. Our measured data are then combined with a mathematical model of a falling body to obtain probability densities on the space of parameters we seek to estimate. We interpret these results and discuss sources of errors in our estimation procedure. © 2013 Society for Industrial and Applied Mathematics.

Parameter Estimation of Spacecraft Fuel Slosh Model

Science.gov (United States)

Gangadharan, Sathya; Sudermann, James; Marlowe, Andrea; Njengam Charles

2004-01-01

Fuel slosh in the upper stages of a spinning spacecraft during launch has been a long standing concern for the success of a space mission. Energy loss through the movement of the liquid fuel in the fuel tank affects the gyroscopic stability of the spacecraft and leads to nutation (wobble) which can cause devastating control issues. The rate at which nutation develops (defined by Nutation Time Constant (NTC can be tedious to calculate and largely inaccurate if done during the early stages of spacecraft design. Pure analytical means of predicting the influence of onboard liquids have generally failed. A strong need exists to identify and model the conditions of resonance between nutation motion and liquid modes and to understand the general characteristics of the liquid motion that causes the problem in spinning spacecraft. A 3-D computerized model of the fuel slosh that accounts for any resonant modes found in the experimental testing will allow for increased accuracy in the overall modeling process. Development of a more accurate model of the fuel slosh currently lies in a more generalized 3-D computerized model incorporating masses, springs and dampers. Parameters describing the model include the inertia tensor of the fuel, spring constants, and damper coefficients. Refinement and understanding the effects of these parameters allow for a more accurate simulation of fuel slosh. The current research will focus on developing models of different complexity and estimating the model parameters that will ultimately provide a more realistic prediction of Nutation Time Constant obtained through simulation.

Inhalation Exposure Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

K. Rautenstrauch

2004-09-10

This analysis is one of 10 reports that support the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN) biosphere model. The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents development of input parameters for the biosphere model that are related to atmospheric mass loading and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for a Yucca Mountain repository. Inhalation Exposure Input Parameters for the Biosphere Model is one of five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the plan for development of the biosphere abstraction products for TSPA, as identified in the Technical Work Plan for Biosphere Modeling and Expert Support (BSC 2004 [DIRS 169573]). This analysis report defines and justifies values of mass loading for the biosphere model. Mass loading is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Mass loading values are used in the air submodel of ERMYN to calculate concentrations of radionuclides in air inhaled by a receptor and concentrations in air surrounding crops. Concentrations in air to which the receptor is exposed are then used in the inhalation submodel to calculate the dose contribution to the receptor from inhalation of contaminated airborne particles. Concentrations in air surrounding plants are used in the plant submodel to calculate the concentrations of radionuclides in foodstuffs contributed from uptake by foliar interception.

Inhalation Exposure Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

K. Rautenstrauch

2004-01-01

This analysis is one of 10 reports that support the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN) biosphere model. The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents development of input parameters for the biosphere model that are related to atmospheric mass loading and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for a Yucca Mountain repository. Inhalation Exposure Input Parameters for the Biosphere Model is one of five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the plan for development of the biosphere abstraction products for TSPA, as identified in the Technical Work Plan for Biosphere Modeling and Expert Support (BSC 2004 [DIRS 169573]). This analysis report defines and justifies values of mass loading for the biosphere model. Mass loading is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Mass loading values are used in the air submodel of ERMYN to calculate concentrations of radionuclides in air inhaled by a receptor and concentrations in air surrounding crops. Concentrations in air to which the receptor is exposed are then used in the inhalation submodel to calculate the dose contribution to the receptor from inhalation of contaminated airborne particles. Concentrations in air surrounding plants are used in the plant submodel to calculate the concentrations of radionuclides in foodstuffs contributed from uptake by foliar interception

Constraining statistical-model parameters using fusion and spallation reactions

Directory of Open Access Journals (Sweden)

Charity Robert J.

2011-10-01

Full Text Available The de-excitation of compound nuclei has been successfully described for several decades by means of statistical models. However, such models involve a large number of free parameters and ingredients that are often underconstrained by experimental data. We show how the degeneracy of the model ingredients can be partially lifted by studying different entrance channels for de-excitation, which populate different regions of the parameter space of the compound nucleus. Fusion reactions, in particular, play an important role in this strategy because they fix three out of four of the compound-nucleus parameters (mass, charge and total excitation energy. The present work focuses on fission and intermediate-mass-fragment emission cross sections. We prove how equivalent parameter sets for fusion-fission reactions can be resolved using another entrance channel, namely spallation reactions. Intermediate-mass-fragment emission can be constrained in a similar way. An interpretation of the best-fit IMF barriers in terms of the Wigner energies of the nascent fragments is discussed.

Agricultural and Environmental Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

K. Rasmuson; K. Rautenstrauch

2004-09-14

This analysis is one of 10 technical reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) (i.e., the biosphere model). It documents development of agricultural and environmental input parameters for the biosphere model, and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for the repository at Yucca Mountain. The ERMYN provides the TSPA with the capability to perform dose assessments. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships between the major activities and their products (the analysis and model reports) that were planned in ''Technical Work Plan for Biosphere Modeling and Expert Support'' (BSC 2004 [DIRS 169573]). The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the ERMYN and its input parameters.

Agricultural and Environmental Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

K. Rasmuson; K. Rautenstrauch

2004-01-01

This analysis is one of 10 technical reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) (i.e., the biosphere model). It documents development of agricultural and environmental input parameters for the biosphere model, and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for the repository at Yucca Mountain. The ERMYN provides the TSPA with the capability to perform dose assessments. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships between the major activities and their products (the analysis and model reports) that were planned in ''Technical Work Plan for Biosphere Modeling and Expert Support'' (BSC 2004 [DIRS 169573]). The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the ERMYN and its input parameters

Soil-Related Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

A. J. Smith

2004-09-09

This report presents one of the analyses that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN). The ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes the details of the conceptual model as well as the mathematical model and the required input parameters. The biosphere model is one of a series of process models supporting the postclosure Total System Performance Assessment (TSPA) for the Yucca Mountain repository. A schematic representation of the documentation flow for the Biosphere input to TSPA is presented in Figure 1-1. This figure shows the evolutionary relationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the biosphere abstraction products for TSPA, as identified in the ''Technical Work Plan for Biosphere Modeling and Expert Support'' (TWP) (BSC 2004 [DIRS 169573]). This figure is included to provide an understanding of how this analysis report contributes to biosphere modeling in support of the license application, and is not intended to imply that access to the listed documents is required to understand the contents of this report. This report, ''Soil-Related Input Parameters for the Biosphere Model'', is one of the five analysis reports that develop input parameters for use in the ERMYN model. This report is the source documentation for the six biosphere parameters identified in Table 1-1. The purpose of this analysis was to develop the biosphere model parameters associated with the accumulation and depletion of radionuclides in the soil. These parameters support the calculation of radionuclide concentrations in soil from on-going irrigation or ash deposition and, as a direct consequence, radionuclide concentration in other environmental media that are affected by radionuclide concentrations in soil. The analysis was performed in accordance with the TWP (BSC 2004 [DIRS 169573]) where the governing procedure

Estimating Parameters for the PVsyst Version 6 Photovoltaic Module Performance Model

Energy Technology Data Exchange (ETDEWEB)

Hansen, Clifford [Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)

2015-10-01

We present an algorithm to determine parameters for the photovoltaic module perf ormance model encoded in the software package PVsyst(TM) version 6. Our method operates on current - voltage (I - V) measured over a range of irradiance and temperature conditions. We describe the method and illustrate its steps using data for a 36 cell crystalli ne silicon module. We qualitatively compare our method with one other technique for estimating parameters for the PVsyst(TM) version 6 model .

An Application of a Multidimensional Extension of the Two-Parameter Logistic Latent Trait Model.

Science.gov (United States)

McKinley, Robert L.; Reckase, Mark D.

A latent trait model is described that is appropriate for use with tests that measure more than one dimension, and its application to both real and simulated test data is demonstrated. Procedures for estimating the parameters of the model are presented. The research objectives are to determine whether the two-parameter logistic model more…

Model description and kinetic parameter analysis of MTBE biodegradation in a packed bed reactor

DEFF Research Database (Denmark)

Waul, Christopher Kevin; Arvin, Erik; Schmidt, Jens Ejbye

2008-01-01

A dynamic modeling approach was used to estimate in-situ model parameters, which describe the degradation of methyl tert-butyl ether (MTBE) in a laboratory packed bed reactor. The measured dynamic response of MTBE pulses injected at the reactor's inlet was analyzed by least squares and parameter...

Experimental kinetic parameters in the thermo-fluid-dynamic modelling of coal combustion

International Nuclear Information System (INIS)

Migliavacca, G.; Perini, M.; Parodi, E.

2001-01-01

The designing and the optimisation of modern and efficient combustion systems are nowadays frequently based on calculation tools for mathematical modelling, which are able to predict the evolution of the process starting from the first principles of physics. Otherwise, in many cases, specific experimental parameters are needed to describe the specific nature of the materials considered in the calculations. It is especially true in the modelling of coal combustion, which is a complex process strongly dependent on the chemical and physical features of the fuel. This paper describes some experimental techniques used to estimate the fundamental kinetic parameters of coal combustion and shows how this data may be introduced in a model calculation to predict the pollutant emissions from a real scale combustion plant [it

Modeling the outflow of liquid with initial supercritical parameters using the relaxation model for condensation

Directory of Open Access Journals (Sweden)

Lezhnin Sergey

2017-01-01

Full Text Available The two-temperature model of the outflow from a vessel with initial supercritical parameters of medium has been realized. The model uses thermodynamic non-equilibrium relaxation approach to describe phase transitions. Based on a new asymptotic model for computing the relaxation time, the outflow of water with supercritical initial pressure and super- and subcritical temperatures has been calculated.

Parameter optimization for surface flux transport models

Science.gov (United States)

Whitbread, T.; Yeates, A. R.; Muñoz-Jaramillo, A.; Petrie, G. J. D.

2017-11-01

Accurate prediction of solar activity calls for precise calibration of solar cycle models. Consequently we aim to find optimal parameters for models which describe the physical processes on the solar surface, which in turn act as proxies for what occurs in the interior and provide source terms for coronal models. We use a genetic algorithm to optimize surface flux transport models using National Solar Observatory (NSO) magnetogram data for Solar Cycle 23. This is applied to both a 1D model that inserts new magnetic flux in the form of idealized bipolar magnetic regions, and also to a 2D model that assimilates specific shapes of real active regions. The genetic algorithm searches for parameter sets (meridional flow speed and profile, supergranular diffusivity, initial magnetic field, and radial decay time) that produce the best fit between observed and simulated butterfly diagrams, weighted by a latitude-dependent error structure which reflects uncertainty in observations. Due to the easily adaptable nature of the 2D model, the optimization process is repeated for Cycles 21, 22, and 24 in order to analyse cycle-to-cycle variation of the optimal solution. We find that the ranges and optimal solutions for the various regimes are in reasonable agreement with results from the literature, both theoretical and observational. The optimal meridional flow profiles for each regime are almost entirely within observational bounds determined by magnetic feature tracking, with the 2D model being able to accommodate the mean observed profile more successfully. Differences between models appear to be important in deciding values for the diffusive and decay terms. In like fashion, differences in the behaviours of different solar cycles lead to contrasts in parameters defining the meridional flow and initial field strength.

Model parameter updating using Bayesian networks

International Nuclear Information System (INIS)

Treml, C.A.; Ross, Timothy J.

2004-01-01

This paper outlines a model parameter updating technique for a new method of model validation using a modified model reference adaptive control (MRAC) framework with Bayesian Networks (BNs). The model parameter updating within this method is generic in the sense that the model/simulation to be validated is treated as a black box. It must have updateable parameters to which its outputs are sensitive, and those outputs must have metrics that can be compared to that of the model reference, i.e., experimental data. Furthermore, no assumptions are made about the statistics of the model parameter uncertainty, only upper and lower bounds need to be specified. This method is designed for situations where a model is not intended to predict a complete point-by-point time domain description of the item/system behavior; rather, there are specific points, features, or events of interest that need to be predicted. These specific points are compared to the model reference derived from actual experimental data. The logic for updating the model parameters to match the model reference is formed via a BN. The nodes of this BN consist of updateable model input parameters and the specific output values or features of interest. Each time the model is executed, the input/output pairs are used to adapt the conditional probabilities of the BN. Each iteration further refines the inferred model parameters to produce the desired model output. After parameter updating is complete and model inputs are inferred, reliabilities for the model output are supplied. Finally, this method is applied to a simulation of a resonance control cooling system for a prototype coupled cavity linac. The results are compared to experimental data.

A New Six-Parameter Model Based on Chebyshev Polynomials for Solar Cells

Directory of Open Access Journals (Sweden)

Shu-xian Lun

2015-01-01

Full Text Available This paper presents a new current-voltage (I-V model for solar cells. It has been proved that series resistance of a solar cell is related to temperature. However, the existing five-parameter model ignores the temperature dependence of series resistance and then only accurately predicts the performance of monocrystalline silicon solar cells. Therefore, this paper uses Chebyshev polynomials to describe the relationship between series resistance and temperature. This makes a new parameter called temperature coefficient for series resistance introduced into the single-diode model. Then, a new six-parameter model for solar cells is established in this paper. This new model can improve the accuracy of the traditional single-diode model and reflect the temperature dependence of series resistance. To validate the accuracy of the six-parameter model in this paper, five kinds of silicon solar cells with different technology types, that is, monocrystalline silicon, polycrystalline silicon, thin film silicon, and tripe-junction amorphous silicon, are tested at different irradiance and temperature conditions. Experiment results show that the six-parameter model proposed in this paper is an I-V model with moderate computational complexity and high precision.

PEM fuel cell model and simulation in Matlab–Simulink based on physical parameters

International Nuclear Information System (INIS)

Abdin, Z.; Webb, C.J.; Gray, E.MacA.

2016-01-01

An advanced PEM fuel cell mathematical model is described and realised in four ancillaries in the Matlab–Simulink environment. Where possible, the model is based on parameters with direct physical meaning, with the aim of going beyond empirically describing the characteristics of the fuel cell. The model can therefore be used to predict enhanced performance owing to, for instance, improved electrode materials, and to relate changes in the measured performance to internal changes affecting influential physical parameters. Some simplifying assumptions make the model fairly light in computational demand and therefore amenable to extension to simulate an entire fuel-cell stack as part of an energy system. Despite these assumptions, the model emulates experimental data well, especially at high current density. The influences of pressure, temperature, humidification and reactant partial pressure on cell performance are explored. The dominating effect of membrane hydration is clearly revealed. - Highlights: • Model based on physical parameters where possible. • Effective binary diffusion modelled in detail on an atomistic basis. • The dominating effect of membrane hydration is clearly revealed. • Documented Simulink model so others can use it. • Conceived as a research tool for exploring enhanced fuel cell performance and diagnosing problems.

Model parameter learning using Kullback-Leibler divergence

Science.gov (United States)

Lin, Chungwei; Marks, Tim K.; Pajovic, Milutin; Watanabe, Shinji; Tung, Chih-kuan

2018-02-01

In this paper, we address the following problem: For a given set of spin configurations whose probability distribution is of the Boltzmann type, how do we determine the model coupling parameters? We demonstrate that directly minimizing the Kullback-Leibler divergence is an efficient method. We test this method against the Ising and XY models on the one-dimensional (1D) and two-dimensional (2D) lattices, and provide two estimators to quantify the model quality. We apply this method to two types of problems. First, we apply it to the real-space renormalization group (RG). We find that the obtained RG flow is sufficiently good for determining the phase boundary (within 1% of the exact result) and the critical point, but not accurate enough for critical exponents. The proposed method provides a simple way to numerically estimate amplitudes of the interactions typically truncated in the real-space RG procedure. Second, we apply this method to the dynamical system composed of self-propelled particles, where we extract the parameter of a statistical model (a generalized XY model) from a dynamical system described by the Viscek model. We are able to obtain reasonable coupling values corresponding to different noise strengths of the Viscek model. Our method is thus able to provide quantitative analysis of dynamical systems composed of self-propelled particles.

PARAMETER ESTIMATION IN BREAD BAKING MODEL

Directory of Open Access Journals (Sweden)

Hadiyanto Hadiyanto

2012-05-01

Full Text Available Bread product quality is highly dependent to the baking process. A model for the development of product quality, which was obtained by using quantitative and qualitative relationships, was calibrated by experiments at a fixed baking temperature of 200°C alone and in combination with 100 W microwave powers. The model parameters were estimated in a stepwise procedure i.e. first, heat and mass transfer related parameters, then the parameters related to product transformations and finally product quality parameters. There was a fair agreement between the calibrated model results and the experimental data. The results showed that the applied simple qualitative relationships for quality performed above expectation. Furthermore, it was confirmed that the microwave input is most meaningful for the internal product properties and not for the surface properties as crispness and color. The model with adjusted parameters was applied in a quality driven food process design procedure to derive a dynamic operation pattern, which was subsequently tested experimentally to calibrate the model. Despite the limited calibration with fixed operation settings, the model predicted well on the behavior under dynamic convective operation and on combined convective and microwave operation. It was expected that the suitability between model and baking system could be improved further by performing calibration experiments at higher temperature and various microwave power levels.  Abstrak  PERKIRAAN PARAMETER DALAM MODEL UNTUK PROSES BAKING ROTI. Kualitas produk roti sangat tergantung pada proses baking yang digunakan. Suatu model yang telah dikembangkan dengan metode kualitatif dan kuantitaif telah dikalibrasi dengan percobaan pada temperatur 200oC dan dengan kombinasi dengan mikrowave pada 100 Watt. Parameter-parameter model diestimasi dengan prosedur bertahap yaitu pertama, parameter pada model perpindahan masa dan panas, parameter pada model transformasi, dan

Frameworks for understanding and describing business models

DEFF Research Database (Denmark)

Nielsen, Christian; Roslender, Robin

2014-01-01

This chapter provides in a chronological fashion an introduction to six frameworks that one can apply to describing, understanding and also potentially innovating business models. These six frameworks have been chosen carefully as they represent six very different perspectives on business models...... and in this manner “complement” each other. There are a multitude of varying frameworks that could be chosen from and we urge the reader to search and trial these for themselves. The six chosen models (year of release in parenthesis) are: • Service-Profit Chain (1994) • Strategic Systems Auditing (1997) • Strategy...... Maps (2001) • Intellectual Capital Statements (2003) • Chesbrough’s framework for Open Business Models (2006) • Business Model Canvas (2008)...

Soil-related Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

A. J. Smith

2003-01-01

This analysis is one of the technical reports containing documentation of the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN), a biosphere model supporting the Total System Performance Assessment (TSPA) for the geologic repository at Yucca Mountain. The biosphere model is one of a series of process models supporting the Total System Performance Assessment (TSPA) for the Yucca Mountain repository. A graphical representation of the documentation hierarchy for the ERMYN biosphere model is presented in Figure 1-1. This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the plan for development of the biosphere abstraction products for TSPA, as identified in the ''Technical Work Plan: for Biosphere Modeling and Expert Support'' (BSC 2003 [163602]). It should be noted that some documents identified in Figure 1-1 may be under development at the time this report is issued and therefore not available. This figure is included to provide an understanding of how this analysis report contributes to biosphere modeling in support of the license application, and is not intended to imply that access to the listed documents is required to understand the contents of this report. This report, ''Soil Related Input Parameters for the Biosphere Model'', is one of the five analysis reports that develop input parameters for use in the ERMYN model. This report is the source documentation for the six biosphere parameters identified in Table 1-1. ''The Biosphere Model Report'' (BSC 2003 [160699]) describes in detail the conceptual model as well as the mathematical model and its input parameters. The purpose of this analysis was to develop the biosphere model parameters needed to evaluate doses from pathways associated with the accumulation and depletion of radionuclides in the soil. These parameters support the calculation of radionuclide concentrations in soil from on-going irrigation and ash

Three-dimensional FEM model of FBGs in PANDA fibers with experimentally determined model parameters

Science.gov (United States)

Lindner, Markus; Hopf, Barbara; Koch, Alexander W.; Roths, Johannes

2017-04-01

A 3D-FEM model has been developed to improve the understanding of multi-parameter sensing with Bragg gratings in attached or embedded polarization maintaining fibers. The material properties of the fiber, especially Young's modulus and Poisson's ratio of the fiber's stress applying parts, are crucial for accurate simulations, but are usually not provided by the manufacturers. A methodology is presented to determine the unknown parameters by using experimental characterizations of the fiber and iterative FEM simulations. The resulting 3D-Model is capable of describing the change in birefringence of the free fiber when exposed to longitudinal strain. In future studies the 3D-FEM model will be employed to study the interaction of PANDA fibers with the surrounding materials in which they are embedded.

Impact of various parameters in detecting chromosomal aberrations by FISH to describe radiosensitivity

International Nuclear Information System (INIS)

Keller, U.; Mueller, E.; Grabenbauer, G.; Sauer, R.; Distel, L.; Kuechler, A.; Liehr, T.

2004-01-01

Background and purpose: analysis of radiation-induced chromosomal aberrations is regarded as the ''gold standard'' for classifying individual radiosensitivity. A variety of different parameters can be used. The crucial question, however, is to explore which parameter is suited best to describe the differences between patients with increased radiosensitivity and healthy individuals. Patients and methods: in this study, five patients with severe radiation-induced late effects of at least grade 3, classified according to the Radiation Therapy Oncology Group (RTOG), and eleven healthy individuals were examined retrospectively. Peripheral blood lymphocytes were irradiated in vitro with 0.7 Gy and 2.0 Gy prior to cultivation and stained by means of three-color fluorescence in situ hybridization (FISH). The detailed analysis was focused on the number of breaks per metaphase, on breaks from complex chromosomal rearrangements per metaphase, as well as on the percentage of translocations, dicentric chromosomes, breaks, and excess acentric fragments - each in comparison with the total number of mitoses analyzed. Results: using the number of breaks from complex chromosomal rearrangements after 2.0 Gy, radiosensitive patients as endpoint were clearly to be distinguished (p = 0.001) from healthy individuals. Translocations (p = 0.001) as well as breaks per metaphase (p = 0.002) were also suitable indicators for detecting differences between patients and healthy individuals. The parameters ''percentage of dicentric chromosomes'', ''breaks'', and ''excess acentric fragments'' in comparison to the total number of mitoses analyzed could neither serve as meaningful nor as significant criteria, since they showed a strong interindividual variability. Conclusion: to detect a difference in chromosomal aberrations between healthy and radiosensitive individuals, the parameters ''frequency of breaks per metaphase'', ''complex chromosomal rearrangements'', and ''translocations'' are most

Parameter identification in a generalized time-harmonic Rayleigh damping model for elastography.

Directory of Open Access Journals (Sweden)

Elijah E W Van Houten

Full Text Available The identifiability of the two damping components of a Generalized Rayleigh Damping model is investigated through analysis of the continuum equilibrium equations as well as a simple spring-mass system. Generalized Rayleigh Damping provides a more diversified attenuation model than pure Viscoelasticity, with two parameters to describe attenuation effects and account for the complex damping behavior found in biological tissue. For heterogeneous Rayleigh Damped materials, there is no equivalent Viscoelastic system to describe the observed motions. For homogeneous systems, the inverse problem to determine the two Rayleigh Damping components is seen to be uniquely posed, in the sense that the inverse matrix for parameter identification is full rank, with certain conditions: when either multi-frequency data is available or when both shear and dilatational wave propagation is taken into account. For the multi-frequency case, the frequency dependency of the elastic parameters adds a level of complexity to the reconstruction problem that must be addressed for reasonable solutions. For the dilatational wave case, the accuracy of compressional wave measurement in fluid saturated soft tissues becomes an issue for qualitative parameter identification. These issues can be addressed with reasonable assumptions on the negligible damping levels of dilatational waves in soft tissue. In general, the parameters of a Generalized Rayleigh Damping model are identifiable for the elastography inverse problem, although with more complex conditions than the simpler Viscoelastic damping model. The value of this approach is the additional structural information provided by the Generalized Rayleigh Damping model, which can be linked to tissue composition as well as rheological interpretations.

Combustion Model and Control Parameter Optimization Methods for Single Cylinder Diesel Engine

Directory of Open Access Journals (Sweden)

Bambang Wahono

2014-01-01

Full Text Available This research presents a method to construct a combustion model and a method to optimize some control parameters of diesel engine in order to develop a model-based control system. The construction purpose of the model is to appropriately manage some control parameters to obtain the values of fuel consumption and emission as the engine output objectives. Stepwise method considering multicollinearity was applied to construct combustion model with the polynomial model. Using the experimental data of a single cylinder diesel engine, the model of power, BSFC, NOx, and soot on multiple injection diesel engines was built. The proposed method succesfully developed the model that describes control parameters in relation to the engine outputs. Although many control devices can be mounted to diesel engine, optimization technique is required to utilize this method in finding optimal engine operating conditions efficiently beside the existing development of individual emission control methods. Particle swarm optimization (PSO was used to calculate control parameters to optimize fuel consumption and emission based on the model. The proposed method is able to calculate control parameters efficiently to optimize evaluation item based on the model. Finally, the model which added PSO then was compiled in a microcontroller.

A Consistent Methodology Based Parameter Estimation for a Lactic Acid Bacteria Fermentation Model

DEFF Research Database (Denmark)

Spann, Robert; Roca, Christophe; Kold, David

2017-01-01

Lactic acid bacteria are used in many industrial applications, e.g. as starter cultures in the dairy industry or as probiotics, and research on their cell production is highly required. A first principles kinetic model was developed to describe and understand the biological, physical, and chemical...... mechanisms in a lactic acid bacteria fermentation. We present here a consistent approach for a methodology based parameter estimation for a lactic acid fermentation. In the beginning, just an initial knowledge based guess of parameters was available and an initial parameter estimation of the complete set...... of parameters was performed in order to get a good model fit to the data. However, not all parameters are identifiable with the given data set and model structure. Sensitivity, identifiability, and uncertainty analysis were completed and a relevant identifiable subset of parameters was determined for a new...

Robust estimation of hydrological model parameters

Directory of Open Access Journals (Sweden)

A. Bárdossy

2008-11-01

Full Text Available The estimation of hydrological model parameters is a challenging task. With increasing capacity of computational power several complex optimization algorithms have emerged, but none of the algorithms gives a unique and very best parameter vector. The parameters of fitted hydrological models depend upon the input data. The quality of input data cannot be assured as there may be measurement errors for both input and state variables. In this study a methodology has been developed to find a set of robust parameter vectors for a hydrological model. To see the effect of observational error on parameters, stochastically generated synthetic measurement errors were applied to observed discharge and temperature data. With this modified data, the model was calibrated and the effect of measurement errors on parameters was analysed. It was found that the measurement errors have a significant effect on the best performing parameter vector. The erroneous data led to very different optimal parameter vectors. To overcome this problem and to find a set of robust parameter vectors, a geometrical approach based on Tukey's half space depth was used. The depth of the set of N randomly generated parameters was calculated with respect to the set with the best model performance (Nash-Sutclife efficiency was used for this study for each parameter vector. Based on the depth of parameter vectors, one can find a set of robust parameter vectors. The results show that the parameters chosen according to the above criteria have low sensitivity and perform well when transfered to a different time period. The method is demonstrated on the upper Neckar catchment in Germany. The conceptual HBV model was used for this study.

Verification Techniques for Parameter Selection and Bayesian Model Calibration Presented for an HIV Model

Science.gov (United States)

Wentworth, Mami Tonoe

Uncertainty quantification plays an important role when making predictive estimates of model responses. In this context, uncertainty quantification is defined as quantifying and reducing uncertainties, and the objective is to quantify uncertainties in parameter, model and measurements, and propagate the uncertainties through the model, so that one can make a predictive estimate with quantified uncertainties. Two of the aspects of uncertainty quantification that must be performed prior to propagating uncertainties are model calibration and parameter selection. There are several efficient techniques for these processes; however, the accuracy of these methods are often not verified. This is the motivation for our work, and in this dissertation, we present and illustrate verification frameworks for model calibration and parameter selection in the context of biological and physical models. First, HIV models, developed and improved by [2, 3, 8], describe the viral infection dynamics of an HIV disease. These are also used to make predictive estimates of viral loads and T-cell counts and to construct an optimal control for drug therapy. Estimating input parameters is an essential step prior to uncertainty quantification. However, not all the parameters are identifiable, implying that they cannot be uniquely determined by the observations. These unidentifiable parameters can be partially removed by performing parameter selection, a process in which parameters that have minimal impacts on the model response are determined. We provide verification techniques for Bayesian model calibration and parameter selection for an HIV model. As an example of a physical model, we employ a heat model with experimental measurements presented in [10]. A steady-state heat model represents a prototypical behavior for heat conduction and diffusion process involved in a thermal-hydraulic model, which is a part of nuclear reactor models. We employ this simple heat model to illustrate verification

Photovoltaic module parameters acquisition model

Energy Technology Data Exchange (ETDEWEB)

Cibira, Gabriel, E-mail: cibira@lm.uniza.sk; Koščová, Marcela, E-mail: mkoscova@lm.uniza.sk

2014-09-01

Highlights: • Photovoltaic five-parameter model is proposed using Matlab{sup ®} and Simulink. • The model acquisits input sparse data matrix from stigmatic measurement. • Computer simulations lead to continuous I–V and P–V characteristics. • Extrapolated I–V and P–V characteristics are in hand. • The model allows us to predict photovoltaics exploitation in different conditions. - Abstract: This paper presents basic procedures for photovoltaic (PV) module parameters acquisition using MATLAB and Simulink modelling. In first step, MATLAB and Simulink theoretical model are set to calculate I–V and P–V characteristics for PV module based on equivalent electrical circuit. Then, limited I–V data string is obtained from examined PV module using standard measurement equipment at standard irradiation and temperature conditions and stated into MATLAB data matrix as a reference model. Next, the theoretical model is optimized to keep-up with the reference model and to learn its basic parameters relations, over sparse data matrix. Finally, PV module parameters are deliverable for acquisition at different realistic irradiation, temperature conditions as well as series resistance. Besides of output power characteristics and efficiency calculation for PV module or system, proposed model validates computing statistical deviation compared to reference model.

Photovoltaic module parameters acquisition model

International Nuclear Information System (INIS)

Cibira, Gabriel; Koščová, Marcela

2014-01-01

Highlights: • Photovoltaic five-parameter model is proposed using Matlab ® and Simulink. • The model acquisits input sparse data matrix from stigmatic measurement. • Computer simulations lead to continuous I–V and P–V characteristics. • Extrapolated I–V and P–V characteristics are in hand. • The model allows us to predict photovoltaics exploitation in different conditions. - Abstract: This paper presents basic procedures for photovoltaic (PV) module parameters acquisition using MATLAB and Simulink modelling. In first step, MATLAB and Simulink theoretical model are set to calculate I–V and P–V characteristics for PV module based on equivalent electrical circuit. Then, limited I–V data string is obtained from examined PV module using standard measurement equipment at standard irradiation and temperature conditions and stated into MATLAB data matrix as a reference model. Next, the theoretical model is optimized to keep-up with the reference model and to learn its basic parameters relations, over sparse data matrix. Finally, PV module parameters are deliverable for acquisition at different realistic irradiation, temperature conditions as well as series resistance. Besides of output power characteristics and efficiency calculation for PV module or system, proposed model validates computing statistical deviation compared to reference model

Inhalation Exposure Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

M. Wasiolek

2006-06-05

This analysis is one of the technical reports that support the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN), referred to in this report as the biosphere model. ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents development of input parameters for the biosphere model that are related to atmospheric mass loading and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for a Yucca Mountain repository. ''Inhalation Exposure Input Parameters for the Biosphere Model'' is one of five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the biosphere model is presented in Figure 1-1 (based on BSC 2006 [DIRS 176938]). This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling and how this analysis report contributes to biosphere modeling. This analysis report defines and justifies values of atmospheric mass loading for the biosphere model. Mass loading is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Mass loading values are used in the air submodel of the biosphere model to calculate concentrations of radionuclides in air inhaled by a receptor and concentrations in air surrounding crops. Concentrations in air to which the receptor is exposed are then used in the inhalation submodel to calculate the dose contribution to the receptor from inhalation of contaminated airborne particles. Concentrations in air surrounding plants are used in the plant submodel to calculate the concentrations of radionuclides in foodstuffs contributed from uptake by foliar interception. This

Inhalation Exposure Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

M. Wasiolek

2006-01-01

This analysis is one of the technical reports that support the Environmental Radiation Model for Yucca Mountain, Nevada (ERMYN), referred to in this report as the biosphere model. ''Biosphere Model Report'' (BSC 2004 [DIRS 169460]) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents development of input parameters for the biosphere model that are related to atmospheric mass loading and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the total system performance assessment (TSPA) for a Yucca Mountain repository. ''Inhalation Exposure Input Parameters for the Biosphere Model'' is one of five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the biosphere model is presented in Figure 1-1 (based on BSC 2006 [DIRS 176938]). This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling and how this analysis report contributes to biosphere modeling. This analysis report defines and justifies values of atmospheric mass loading for the biosphere model. Mass loading is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Mass loading values are used in the air submodel of the biosphere model to calculate concentrations of radionuclides in air inhaled by a receptor and concentrations in air surrounding crops. Concentrations in air to which the receptor is exposed are then used in the inhalation submodel to calculate the dose contribution to the receptor from inhalation of contaminated airborne particles. Concentrations in air surrounding plants are used in the plant submodel to calculate the concentrations of radionuclides in foodstuffs contributed from uptake by foliar interception. This report is concerned primarily with the

Parameters of Models of Structural Transformations in Alloy Steel Under Welding Thermal Cycle

Science.gov (United States)

Kurkin, A. S.; Makarov, E. L.; Kurkin, A. B.; Rubtsov, D. E.; Rubtsov, M. E.

2017-05-01

A mathematical model of structural transformations in an alloy steel under the thermal cycle of multipass welding is suggested for computer implementation. The minimum necessary set of parameters for describing the transformations under heating and cooling is determined. Ferritic-pearlitic, bainitic and martensitic transformations under cooling of a steel are considered. A method for deriving the necessary temperature and time parameters of the model from the chemical composition of the steel is described. Published data are used to derive regression models of the temperature ranges and parameters of transformation kinetics in alloy steels. It is shown that the disadvantages of the active visual methods of analysis of the final phase composition of steels are responsible for inaccuracy and mismatch of published data. The hardness of a specimen, which correlates with some other mechanical properties of the material, is chosen as the most objective and reproducible criterion of the final phase composition. The models developed are checked by a comparative analysis of computational results and experimental data on the hardness of 140 alloy steels after cooling at various rates.

Inhalation Exposure Input Parameters for the Biosphere Model

International Nuclear Information System (INIS)

M. A. Wasiolek

2003-01-01

This analysis is one of the nine reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) biosphere model. The ''Biosphere Model Report'' (BSC 2003a) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents a set of input parameters for the biosphere model, and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the Total System Performance Assessment (TSPA) for a Yucca Mountain repository. This report, ''Inhalation Exposure Input Parameters for the Biosphere Model'', is one of the five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the plan for development of the biosphere abstraction products for TSPA, as identified in the ''Technical Work Plan: for Biosphere Modeling and Expert Support'' (BSC 2003b). It should be noted that some documents identified in Figure 1-1 may be under development at the time this report is issued and therefore not available at that time. This figure is included to provide an understanding of how this analysis report contributes to biosphere modeling in support of the license application, and is not intended to imply that access to the listed documents is required to understand the contents of this analysis report. This analysis report defines and justifies values of mass loading, which is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Measurements of mass loading are used in the air submodel of ERMYN to calculate concentrations of radionuclides in air surrounding crops and concentrations in air inhaled by a receptor. Concentrations in air to which the

A note on modeling of tumor regression for estimation of radiobiological parameters

International Nuclear Information System (INIS)

Zhong, Hualiang; Chetty, Indrin

2014-01-01

Purpose: Accurate calculation of radiobiological parameters is crucial to predicting radiation treatment response. Modeling differences may have a significant impact on derived parameters. In this study, the authors have integrated two existing models with kinetic differential equations to formulate a new tumor regression model for estimation of radiobiological parameters for individual patients. Methods: A system of differential equations that characterizes the birth-and-death process of tumor cells in radiation treatment was analytically solved. The solution of this system was used to construct an iterative model (Z-model). The model consists of three parameters: tumor doubling time T d , half-life of dead cells T r , and cell survival fraction SF D under dose D. The Jacobian determinant of this model was proposed as a constraint to optimize the three parameters for six head and neck cancer patients. The derived parameters were compared with those generated from the two existing models: Chvetsov's model (C-model) and Lim's model (L-model). The C-model and L-model were optimized with the parameter T d fixed. Results: With the Jacobian-constrained Z-model, the mean of the optimized cell survival fractions is 0.43 ± 0.08, and the half-life of dead cells averaged over the six patients is 17.5 ± 3.2 days. The parameters T r and SF D optimized with the Z-model differ by 1.2% and 20.3% from those optimized with the T d -fixed C-model, and by 32.1% and 112.3% from those optimized with the T d -fixed L-model, respectively. Conclusions: The Z-model was analytically constructed from the differential equations of cell populations that describe changes in the number of different tumor cells during the course of radiation treatment. The Jacobian constraints were proposed to optimize the three radiobiological parameters. The generated model and its optimization method may help develop high-quality treatment regimens for individual patients

A simple but accurate procedure for solving the five-parameter model

International Nuclear Information System (INIS)

Mares, Oana; Paulescu, Marius; Badescu, Viorel

2015-01-01

Highlights: • A new procedure for extracting the parameters of the one-diode model is proposed. • Only the basic information listed in the datasheet of PV modules are required. • Results demonstrate a simple, robust and accurate procedure. - Abstract: The current–voltage characteristic of a photovoltaic module is typically evaluated by using a model based on the solar cell equivalent circuit. The complexity of the procedure applied for extracting the model parameters depends on data available in manufacture’s datasheet. Since the datasheet is not detailed enough, simplified models have to be used in many cases. This paper proposes a new procedure for extracting the parameters of the one-diode model in standard test conditions, using only the basic data listed by all manufactures in datasheet (short circuit current, open circuit voltage and maximum power point). The procedure is validated by using manufacturers’ data for six commercially crystalline silicon photovoltaic modules. Comparing the computed and measured current–voltage characteristics the determination coefficient is in the range 0.976–0.998. Thus, the proposed procedure represents a feasible tool for solving the five-parameter model applied to crystalline silicon photovoltaic modules. The procedure is described in detail, to guide potential users to derive similar models for other types of photovoltaic modules.

Elastic-plastic analysis of AS4/PEEK composite laminate using a one-parameter plasticity model

Science.gov (United States)

Sun, C. T.; Yoon, K. J.

1992-01-01

A one-parameter plasticity model was shown to adequately describe the plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The elastic-plastic stress-strain relations of coupon specimens were measured and compared with those predicted by the finite element analysis using the one-parameter plasticity model. The results show that the one-parameter plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Nonlinear analysis of AS4/PEEK thermoplastic composite laminate using a one parameter plasticity model

Science.gov (United States)

Sun, C. T.; Yoon, K. J.

1990-01-01

A one-parameter plasticity model was shown to adequately describe the orthotropic plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The nonlinear stress-strain relations were measured and compared with those predicted by the finite element analysis using the one-parameter elastic-plastic constitutive model. The results show that the one-parameter orthotropic plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Model calibration and parameter estimation for environmental and water resource systems

CERN Document Server

Sun, Ne-Zheng

2015-01-01

This three-part book provides a comprehensive and systematic introduction to the development of useful models for complex systems. Part 1 covers the classical inverse problem for parameter estimation in both deterministic and statistical frameworks, Part 2 is dedicated to system identification, hyperparameter estimation, and model dimension reduction, and Part 3 considers how to collect data and construct reliable models for prediction and decision-making. For the first time, topics such as multiscale inversion, stochastic field parameterization, level set method, machine learning, global sensitivity analysis, data assimilation, model uncertainty quantification, robust design, and goal-oriented modeling, are systematically described and summarized in a single book from the perspective of model inversion, and elucidated with numerical examples from environmental and water resources modeling. Readers of this book will not only learn basic concepts and methods for simple parameter estimation, but also get famili...

Multi-objective genetic algorithm parameter estimation in a reduced nuclear reactor model

Energy Technology Data Exchange (ETDEWEB)

Marseguerra, M.; Zio, E.; Canetta, R. [Polytechnic of Milan, Dept. of Nuclear Engineering, Milano (Italy)

2005-07-01

The fast increase in computing power has rendered, and will continue to render, more and more feasible the incorporation of dynamics in the safety and reliability models of complex engineering systems. In particular, the Monte Carlo simulation framework offers a natural environment for estimating the reliability of systems with dynamic features. However, the time-integration of the dynamic processes may render the Monte Carlo simulation quite burdensome so that it becomes mandatory to resort to validated, simplified models of process evolution. Such models are typically based on lumped effective parameters whose values need to be suitably estimated so as to best fit to the available plant data. In this paper we propose a multi-objective genetic algorithm approach for the estimation of the effective parameters of a simplified model of nuclear reactor dynamics. The calibration of the effective parameters is achieved by best fitting the model responses of the quantities of interest to the actual evolution profiles. A case study is reported in which the real reactor is simulated by the QUAndry based Reactor Kinetics (Quark) code available from the Nuclear Energy Agency and the simplified model is based on the point kinetics approximation to describe the neutron balance in the core and on thermal equilibrium relations to describe the energy exchange between the different loops. (authors)

Multi-objective genetic algorithm parameter estimation in a reduced nuclear reactor model

International Nuclear Information System (INIS)

Marseguerra, M.; Zio, E.; Canetta, R.

2005-01-01

The fast increase in computing power has rendered, and will continue to render, more and more feasible the incorporation of dynamics in the safety and reliability models of complex engineering systems. In particular, the Monte Carlo simulation framework offers a natural environment for estimating the reliability of systems with dynamic features. However, the time-integration of the dynamic processes may render the Monte Carlo simulation quite burdensome so that it becomes mandatory to resort to validated, simplified models of process evolution. Such models are typically based on lumped effective parameters whose values need to be suitably estimated so as to best fit to the available plant data. In this paper we propose a multi-objective genetic algorithm approach for the estimation of the effective parameters of a simplified model of nuclear reactor dynamics. The calibration of the effective parameters is achieved by best fitting the model responses of the quantities of interest to the actual evolution profiles. A case study is reported in which the real reactor is simulated by the QUAndry based Reactor Kinetics (Quark) code available from the Nuclear Energy Agency and the simplified model is based on the point kinetics approximation to describe the neutron balance in the core and on thermal equilibrium relations to describe the energy exchange between the different loops. (authors)

The lumped parameter model for fuel pins

Energy Technology Data Exchange (ETDEWEB)

Liu, W S [Ontario Hydro, Toronto, ON (Canada)

1996-12-31

The use of a lumped fuel-pin model in a thermal-hydraulic code is advantageous because of computational simplicity and efficiency. The model uses an averaging approach over the fuel cross section and makes some simplifying assumptions to describe the transient equations for the averaged fuel, fuel centerline and sheath temperatures. It is shown that by introducing a factor in the effective fuel conductivity, the analytical solution of the mean fuel temperature can be modified to simulate the effects of the flux depression in the heat generation rate and the variation in fuel thermal conductivity. The simplified analytical method used in the transient equation is presented. The accuracy of the lumped parameter model has been compared with the results from the finite difference method. (author). 4 refs., 2 tabs., 4 figs.

Prediction of interest rate using CKLS model with stochastic parameters

International Nuclear Information System (INIS)

Ying, Khor Chia; Hin, Pooi Ah

2014-01-01

The Chan, Karolyi, Longstaff and Sanders (CKLS) model is a popular one-factor model for describing the spot interest rates. In this paper, the four parameters in the CKLS model are regarded as stochastic. The parameter vector φ (j) of four parameters at the (J+n)-th time point is estimated by the j-th window which is defined as the set consisting of the observed interest rates at the j′-th time point where j≤j′≤j+n. To model the variation of φ (j) , we assume that φ (j) depends on φ (j−m) , φ (j−m+1) ,…, φ (j−1) and the interest rate r j+n at the (j+n)-th time point via a four-dimensional conditional distribution which is derived from a [4(m+1)+1]-dimensional power-normal distribution. Treating the (j+n)-th time point as the present time point, we find a prediction interval for the future value r j+n+1 of the interest rate at the next time point when the value r j+n of the interest rate is given. From the above four-dimensional conditional distribution, we also find a prediction interval for the future interest rate r j+n+d at the next d-th (d≥2) time point. The prediction intervals based on the CKLS model with stochastic parameters are found to have better ability of covering the observed future interest rates when compared with those based on the model with fixed parameters

Prediction of interest rate using CKLS model with stochastic parameters

Energy Technology Data Exchange (ETDEWEB)

Ying, Khor Chia [Faculty of Computing and Informatics, Multimedia University, Jalan Multimedia, 63100 Cyberjaya, Selangor (Malaysia); Hin, Pooi Ah [Sunway University Business School, No. 5, Jalan Universiti, Bandar Sunway, 47500 Subang Jaya, Selangor (Malaysia)

2014-06-19

The Chan, Karolyi, Longstaff and Sanders (CKLS) model is a popular one-factor model for describing the spot interest rates. In this paper, the four parameters in the CKLS model are regarded as stochastic. The parameter vector φ{sup (j)} of four parameters at the (J+n)-th time point is estimated by the j-th window which is defined as the set consisting of the observed interest rates at the j′-th time point where j≤j′≤j+n. To model the variation of φ{sup (j)}, we assume that φ{sup (j)} depends on φ{sup (j−m)}, φ{sup (j−m+1)},…, φ{sup (j−1)} and the interest rate r{sub j+n} at the (j+n)-th time point via a four-dimensional conditional distribution which is derived from a [4(m+1)+1]-dimensional power-normal distribution. Treating the (j+n)-th time point as the present time point, we find a prediction interval for the future value r{sub j+n+1} of the interest rate at the next time point when the value r{sub j+n} of the interest rate is given. From the above four-dimensional conditional distribution, we also find a prediction interval for the future interest rate r{sub j+n+d} at the next d-th (d≥2) time point. The prediction intervals based on the CKLS model with stochastic parameters are found to have better ability of covering the observed future interest rates when compared with those based on the model with fixed parameters.

A morphing technique for signal modelling in a multidimensional space of coupling parameters

CERN Document Server

The ATLAS collaboration

2015-01-01

This note describes a morphing method that produces signal models for fits to data in which both the affected event yields and kinematic distributions are simultaneously taken into account. The signal model is morphed in a continuous manner through the available multi-dimensional parameter space. Searches for deviations from Standard Model predictions for Higgs boson properties have so far used information either from event yields or kinematic distributions. The combined approach described here is expected to substantially enhance the sensitivity to beyond the Standard Model contributions.

A termination criterion for parameter estimation in stochastic models in systems biology.

Science.gov (United States)

Zimmer, Christoph; Sahle, Sven

2015-11-01

Parameter estimation procedures are a central aspect of modeling approaches in systems biology. They are often computationally expensive, especially when the models take stochasticity into account. Typically parameter estimation involves the iterative optimization of an objective function that describes how well the model fits some measured data with a certain set of parameter values. In order to limit the computational expenses it is therefore important to apply an adequate stopping criterion for the optimization process, so that the optimization continues at least until a reasonable fit is obtained, but not much longer. In the case of stochastic modeling, at least some parameter estimation schemes involve an objective function that is itself a random variable. This means that plain convergence tests are not a priori suitable as stopping criteria. This article suggests a termination criterion suited to optimization problems in parameter estimation arising from stochastic models in systems biology. The termination criterion is developed for optimization algorithms that involve populations of parameter sets, such as particle swarm or evolutionary algorithms. It is based on comparing the variance of the objective function over the whole population of parameter sets with the variance of repeated evaluations of the objective function at the best parameter set. The performance is demonstrated for several different algorithms. To test the termination criterion we choose polynomial test functions as well as systems biology models such as an Immigration-Death model and a bistable genetic toggle switch. The genetic toggle switch is an especially challenging test case as it shows a stochastic switching between two steady states which is qualitatively different from the model behavior in a deterministic model. Copyright © 2015. Published by Elsevier Ireland Ltd.

The sensitivity of flowline models of tidewater glaciers to parameter uncertainty

Directory of Open Access Journals (Sweden)

E. M. Enderlin

2013-10-01

Full Text Available Depth-integrated (1-D flowline models have been widely used to simulate fast-flowing tidewater glaciers and predict change because the continuous grounding line tracking, high horizontal resolution, and physically based calving criterion that are essential to realistic modeling of tidewater glaciers can easily be incorporated into the models while maintaining high computational efficiency. As with all models, the values for parameters describing ice rheology and basal friction must be assumed and/or tuned based on observations. For prognostic studies, these parameters are typically tuned so that the glacier matches observed thickness and speeds at an initial state, to which a perturbation is applied. While it is well know that ice flow models are sensitive to these parameters, the sensitivity of tidewater glacier models has not been systematically investigated. Here we investigate the sensitivity of such flowline models of outlet glacier dynamics to uncertainty in three key parameters that influence a glacier's resistive stress components. We find that, within typical observational uncertainty, similar initial (i.e., steady-state glacier configurations can be produced with substantially different combinations of parameter values, leading to differing transient responses after a perturbation is applied. In cases where the glacier is initially grounded near flotation across a basal over-deepening, as typically observed for rapidly changing glaciers, these differences can be dramatic owing to the threshold of stability imposed by the flotation criterion. The simulated transient response is particularly sensitive to the parameterization of ice rheology: differences in ice temperature of ~ 2 °C can determine whether the glaciers thin to flotation and retreat unstably or remain grounded on a marine shoal. Due to the highly non-linear dependence of tidewater glaciers on model parameters, we recommend that their predictions are accompanied by

Parameter estimation and sensitivity analysis for a mathematical model with time delays of leukemia

Science.gov (United States)

Cândea, Doina; Halanay, Andrei; Rǎdulescu, Rodica; Tǎlmaci, Rodica

2017-01-01

We consider a system of nonlinear delay differential equations that describes the interaction between three competing cell populations: healthy, leukemic and anti-leukemia T cells involved in Chronic Myeloid Leukemia (CML) under treatment with Imatinib. The aim of this work is to establish which model parameters are the most important in the success or failure of leukemia remission under treatment using a sensitivity analysis of the model parameters. For the most significant parameters of the model which affect the evolution of CML disease during Imatinib treatment we try to estimate the realistic values using some experimental data. For these parameters, steady states are calculated and their stability is analyzed and biologically interpreted.

Equation-free analysis of agent-based models and systematic parameter determination

Science.gov (United States)

Thomas, Spencer A.; Lloyd, David J. B.; Skeldon, Anne C.

2016-12-01

Agent based models (ABM)s are increasingly used in social science, economics, mathematics, biology and computer science to describe time dependent systems in circumstances where a description in terms of equations is difficult. Yet few tools are currently available for the systematic analysis of ABM behaviour. Numerical continuation and bifurcation analysis is a well-established tool for the study of deterministic systems. Recently, equation-free (EF) methods have been developed to extend numerical continuation techniques to systems where the dynamics are described at a microscopic scale and continuation of a macroscopic property of the system is considered. To date, the practical use of EF methods has been limited by; (1) the over-head of application-specific implementation; (2) the laborious configuration of problem-specific parameters; and (3) large ensemble sizes (potentially) leading to computationally restrictive run-times. In this paper we address these issues with our tool for the EF continuation of stochastic systems, which includes algorithms to systematically configuration problem specific parameters and enhance robustness to noise. Our tool is generic and can be applied to any 'black-box' simulator and determines the essential EF parameters prior to EF analysis. Robustness is significantly improved using our convergence-constraint with a corrector-repeat (C3R) method. This algorithm automatically detects outliers based on the dynamics of the underlying system enabling both an order of magnitude reduction in ensemble size and continuation of systems at much higher levels of noise than classical approaches. We demonstrate our method with application to several ABM models, revealing parameter dependence, bifurcation and stability analysis of these complex systems giving a deep understanding of the dynamical behaviour of the models in a way that is not otherwise easily obtainable. In each case we demonstrate our systematic parameter determination stage for

A critical revisit of the key parameters used to describe microbial electrochemical systems

International Nuclear Information System (INIS)

Sharma, Mohita; Bajracharya, Suman; Gildemyn, Sylvia; Patil, Sunil A.; Alvarez-Gallego, Yolanda; Pant, Deepak; Rabaey, Korneel; Dominguez-Benetton, Xochitl

2014-01-01

Graphical abstract: - Abstract: Many microorganisms have the innate capability to discharge and/or receive electrons to and from solid state materials such as electrodes. This ability is now used towards innovative processes in wastewater treatment, power generation, production of fuels and biochemicals, bioremediation, desalination and resource recovery, among others. Despite being a dynamic field in science and technology, significant challenges remain towards industrial implementation which include representation of judicious performance indicators. This critical review outlines the progress in current density evaluated per projected surface area of electrodes, the most wide-spread performance indicator. It also proposes guidelines to correct current and exchange current per porous surface area, biofilm covered area, electrochemically- or bioelectrochemically- active surface area, of the electrodes. Recommendations for indicators to describe the environmental and electrochemical robustness of electrochemically-active biofilms are portrayed, including preservation of the predominant functionality as well as electrochemical mechanistic and phenomenological features. A few additional key elements for industrial processing are depicted. Whereas Microbial Fuel Cells (MFCs) are the main focus, some important parameters for reporting on cathodic bioproduction performance are also discussed. This critical revision aims to provide key parameters to compare the whole spectrum of microbial electrochemical systems in a consistent way

Parameters Tuning of Model Free Adaptive Control Based on Minimum Entropy

Institute of Scientific and Technical Information of China (English)

Chao Ji; Jing Wang; Liulin Cao; Qibing Jin

2014-01-01

Dynamic linearization based model free adaptive control(MFAC) algorithm has been widely used in practical systems, in which some parameters should be tuned before it is successfully applied to process industries. Considering the random noise existing in real processes, a parameter tuning method based on minimum entropy optimization is proposed,and the feature of entropy is used to accurately describe the system uncertainty. For cases of Gaussian stochastic noise and non-Gaussian stochastic noise, an entropy recursive optimization algorithm is derived based on approximate model or identified model. The extensive simulation results show the effectiveness of the minimum entropy optimization for the partial form dynamic linearization based MFAC. The parameters tuned by the minimum entropy optimization index shows stronger stability and more robustness than these tuned by other traditional index,such as integral of the squared error(ISE) or integral of timeweighted absolute error(ITAE), when the system stochastic noise exists.

Parameter Estimation of Partial Differential Equation Models.

Science.gov (United States)

Xun, Xiaolei; Cao, Jiguo; Mallick, Bani; Carroll, Raymond J; Maity, Arnab

2013-01-01

Partial differential equation (PDE) models are commonly used to model complex dynamic systems in applied sciences such as biology and finance. The forms of these PDE models are usually proposed by experts based on their prior knowledge and understanding of the dynamic system. Parameters in PDE models often have interesting scientific interpretations, but their values are often unknown, and need to be estimated from the measurements of the dynamic system in the present of measurement errors. Most PDEs used in practice have no analytic solutions, and can only be solved with numerical methods. Currently, methods for estimating PDE parameters require repeatedly solving PDEs numerically under thousands of candidate parameter values, and thus the computational load is high. In this article, we propose two methods to estimate parameters in PDE models: a parameter cascading method and a Bayesian approach. In both methods, the underlying dynamic process modeled with the PDE model is represented via basis function expansion. For the parameter cascading method, we develop two nested levels of optimization to estimate the PDE parameters. For the Bayesian method, we develop a joint model for data and the PDE, and develop a novel hierarchical model allowing us to employ Markov chain Monte Carlo (MCMC) techniques to make posterior inference. Simulation studies show that the Bayesian method and parameter cascading method are comparable, and both outperform other available methods in terms of estimation accuracy. The two methods are demonstrated by estimating parameters in a PDE model from LIDAR data.

A mathematical model for describing the mechanical behaviour of root canal instruments.

Science.gov (United States)

Zhang, E W; Cheung, G S P; Zheng, Y F

2011-01-01

The purpose of this study was to establish a general mathematical model for describing the mechanical behaviour of root canal instruments by combining a theoretical analytical approach with a numerical finite-element method. Mathematical formulas representing the longitudinal (taper, helical angle and pitch) and cross-sectional configurations and area, the bending and torsional inertia, the curvature of the boundary point and the (geometry of) loading condition were derived. Torsional and bending stresses and the resultant deformation were expressed mathematically as a function of these geometric parameters, modulus of elasticity of the material and the applied load. As illustrations, three brands of NiTi endodontic files of different cross-sectional configurations (ProTaper, Hero 642, and Mani NRT) were analysed under pure torsion and pure bending situation by entering the model into a finite-element analysis package (ANSYS). Numerical results confirmed that mathematical models were a feasible method to analyse the mechanical properties and predict the stress and deformation for root canal instruments during root canal preparation. Mathematical and numerical model can be a suitable way to examine mechanical behaviours as a criterion of the instrument design and to predict the stress and strain experienced by the endodontic instruments during root canal preparation. © 2010 International Endodontic Journal.

Large-scale parameter extraction in electrocardiology models through Born approximation

KAUST Repository

He, Yuan

2012-12-04

One of the main objectives in electrocardiology is to extract physical properties of cardiac tissues from measured information on electrical activity of the heart. Mathematically, this is an inverse problem for reconstructing coefficients in electrocardiology models from partial knowledge of the solutions of the models. In this work, we consider such parameter extraction problems for two well-studied electrocardiology models: the bidomain model and the FitzHugh-Nagumo model. We propose a systematic reconstruction method based on the Born approximation of the original nonlinear inverse problem. We describe a two-step procedure that allows us to reconstruct not only perturbations of the unknowns, but also the backgrounds around which the linearization is performed. We show some numerical simulations under various conditions to demonstrate the performance of our method. We also introduce a parameterization strategy using eigenfunctions of the Laplacian operator to reduce the number of unknowns in the parameter extraction problem. © 2013 IOP Publishing Ltd.

Synchronous Generator Model Parameter Estimation Based on Noisy Dynamic Waveforms

Science.gov (United States)

Berhausen, Sebastian; Paszek, Stefan

2016-01-01

In recent years, there have occurred system failures in many power systems all over the world. They have resulted in a lack of power supply to a large number of recipients. To minimize the risk of occurrence of power failures, it is necessary to perform multivariate investigations, including simulations, of power system operating conditions. To conduct reliable simulations, the current base of parameters of the models of generating units, containing the models of synchronous generators, is necessary. In the paper, there is presented a method for parameter estimation of a synchronous generator nonlinear model based on the analysis of selected transient waveforms caused by introducing a disturbance (in the form of a pseudorandom signal) in the generator voltage regulation channel. The parameter estimation was performed by minimizing the objective function defined as a mean square error for deviations between the measurement waveforms and the waveforms calculated based on the generator mathematical model. A hybrid algorithm was used for the minimization of the objective function. In the paper, there is described a filter system used for filtering the noisy measurement waveforms. The calculation results of the model of a 44 kW synchronous generator installed on a laboratory stand of the Institute of Electrical Engineering and Computer Science of the Silesian University of Technology are also given. The presented estimation method can be successfully applied to parameter estimation of different models of high-power synchronous generators operating in a power system.

Reuse-centric Requirements Analysis with Task Models, Scenarios, and Critical Parameters

Directory of Open Access Journals (Sweden)

Cyril Montabert

2007-02-01

Full Text Available This paper outlines a requirements-analysis process that unites task models, scenarios, and critical parameters to exploit and generate reusable knowledge at the requirements phase. Through the deployment of a critical-parameter-based approach to task modeling, the process yields the establishment of an integrative and formalized model issued from scenarios that can be used for requirements characterization. Furthermore, not only can this entity serve as interface to a knowledge repository relying on a critical-parameter-based taxonomy to support reuse but its characterization in terms of critical parameters also allows the model to constitute a broader reuse solution. We discuss our vision for a user-centric and reuse-centric approach to requirements analysis, present previous efforts implicated with this line of work, and state the revisions brought to extend the reuse potential and effectiveness of a previous iteration of a requirements tool implementing such process. Finally, the paper describes the sequence and nature of the activities involved with the conduct of our proposed requirements-analysis technique, concluding by previewing ongoing work in the field that will explore the feasibility for designers to use our approach.

Digital Materials - Evaluation of the Possibilities of using Selected Hyperelastic Models to Describe Constitutive Relations

Science.gov (United States)

Mańkowski, J.; Lipnicki, J.

2017-08-01

The authors tried to identify the parameters of numerical models of digital materials, which are a kind of composite resulting from the manufacture of the product in 3D printers. With the arrangement of several heads of the printer, the new material can result from mixing of materials with radically different properties, during the process of producing single layer of the product. The new material has properties dependent on the base materials properties and their proportions. Digital materials tensile characteristics are often non-linear and qualify to be described by hyperelastic materials models. The identification was conducted based on the results of tensile tests models, its various degrees coefficients of the polynomials to various degrees coefficients of the polynomials. The Drucker's stability criterion was also examined. Fourteen different materials were analyzed.

The generalized model of polypeptide chain describing the helix-coil transition in biopolymers

International Nuclear Information System (INIS)

Mamasakhlisov, E.S.; Badasyan, A.V.; Tsarukyan, A.V.; Grigoryan, A.V.; Morozov, V.F.

2005-07-01

In this paper we summarize some results of our theoretical investigations of helix-coil transition both in single-strand (polypeptides) and two-strand (polynucleotides) macromolecules. The Hamiltonian of the Generalized Model of Polypeptide Chain (GMPC) is introduced to describe the system in which the conformations are correlated over some dimensional range Δ (it equals 3 for polypeptide, because one H-bond fixes three pairs of rotation, for double strand DNA it equals to one chain rigidity because of impossibility of loop formation on the scale less than Δ). The Hamiltonian does not contain any parameter designed especially for helix-coil transition and uses pure molecular microscopic parameters (the energy of hydrogen bond formation, reduced partition function of repeated unit, the number of repeated units fixed by one hydrogen bond, the energies of interaction between the repeated units and the solvent molecules). To calculate averages we evaluate the partition function using the transfer-matrix approach. The GMPC allowed to describe the influence of a number of factors, affecting the transition, basing on a unified microscopic approach. Thus we obtained, that solvents change transition temperature and interval in different ways, depending on type of solvent and on energy of solvent- macromolecule interaction; stacking on the background of H-bonding increases stability and decreases cooperativity of melting. For heterogeneous DNA we could analytically derive well known formulae for transition temperature and interval. In the framework of GMPC we calculate and show the difference of two order parameters of helix-coil transition - the helicity degree, and the average fraction of repeated units in helical conformation. Given article has the aim to review the results obtained during twenty years in the context of GMPC. (author)

Quality assessment for radiological model parameters

International Nuclear Information System (INIS)

Funtowicz, S.O.

1989-01-01

A prototype framework for representing uncertainties in radiological model parameters is introduced. This follows earlier development in this journal of a corresponding framework for representing uncertainties in radiological data. Refinements and extensions to the earlier framework are needed in order to take account of the additional contextual factors consequent on using data entries to quantify model parameters. The parameter coding can in turn feed in to methods for evaluating uncertainties in calculated model outputs. (author)

Simplified distributed parameters BWR dynamic model for transient and stability analysis

International Nuclear Information System (INIS)

Espinosa-Paredes, Gilberto; Nunez-Carrera, Alejandro; Vazquez-Rodriguez, Alejandro

2006-01-01

This paper describes a simplified model to perform transient and linear stability analysis for a typical boiling water reactor (BWR). The simplified transient model was based in lumped and distributed parameters approximations, which includes vessel dome and the downcomer, recirculation loops, neutron process, fuel pin temperature distribution, lower and upper plenums reactor core and pressure and level controls. The stability was determined by studying the linearized versions of the equations representing the BWR system in the frequency domain. Numerical examples are used to illustrate the wide application of the simplified BWR model. We concluded that this simplified model describes properly the dynamic of a BWR and can be used for safety analysis or as a first approach in the design of an advanced BWR

Inhalation Exposure Input Parameters for the Biosphere Model

Energy Technology Data Exchange (ETDEWEB)

M. A. Wasiolek

2003-09-24

This analysis is one of the nine reports that support the Environmental Radiation Model for Yucca Mountain Nevada (ERMYN) biosphere model. The ''Biosphere Model Report'' (BSC 2003a) describes in detail the conceptual model as well as the mathematical model and its input parameters. This report documents a set of input parameters for the biosphere model, and supports the use of the model to develop biosphere dose conversion factors (BDCFs). The biosphere model is one of a series of process models supporting the Total System Performance Assessment (TSPA) for a Yucca Mountain repository. This report, ''Inhalation Exposure Input Parameters for the Biosphere Model'', is one of the five reports that develop input parameters for the biosphere model. A graphical representation of the documentation hierarchy for the ERMYN is presented in Figure 1-1. This figure shows the interrelationships among the products (i.e., analysis and model reports) developed for biosphere modeling, and the plan for development of the biosphere abstraction products for TSPA, as identified in the ''Technical Work Plan: for Biosphere Modeling and Expert Support'' (BSC 2003b). It should be noted that some documents identified in Figure 1-1 may be under development at the time this report is issued and therefore not available at that time. This figure is included to provide an understanding of how this analysis report contributes to biosphere modeling in support of the license application, and is not intended to imply that access to the listed documents is required to understand the contents of this analysis report. This analysis report defines and justifies values of mass loading, which is the total mass concentration of resuspended particles (e.g., dust, ash) in a volume of air. Measurements of mass loading are used in the air submodel of ERMYN to calculate concentrations of radionuclides in air surrounding crops and concentrations in air

Determination of Kinetic and Thermodynamic Parameters that Describe Isothermal Seed Germination: A Student Research Project.

Science.gov (United States)

Hageseth, Gaylord T.

1982-01-01

Describes a project for students to collect and fit data to a theoretical mathematical model that describes the rate of isothermal seed germination, including activation energy for substrate and produce and the autocatalytic reaction, and changes in enthalpy, entropy, and the Gibb's free energy. (Author/SK)

A distributed parameter wire model for transient electrical discharges

International Nuclear Information System (INIS)

Maier, W.B. II; Kadish, A.; Sutherland, C.D.; Robiscoe, R.T.

1990-01-01

A model for freely propagating transient electrical discharges, such as lightning and punch-through arcs, is developed in this paper. We describe the electromagnetic fields by Maxwell's equations and we represent the interaction of electric fields with the medium to produce current by ∂J/∂t=ω 2 (E-E*J)/4π, where ω and E* are parameters characteristic of the medium, J≡current density, and J≡J/|J|. We illustrate the properties of this model for small-diameter, guided, cylindrically symmetric discharges. Analytic, numerical, and approximate solutions are given for special cases. The model describes, in a new and comprehensive fashion, certain macroscopic discharge properties, such as threshold behavior, quenching and reignition, path tortuosity, discharge termination with nonzero charge density remaining along the discharge path, and other experimentally observed discharge phenomena. Fields, current densities, and charge densities are quantitatively determined from given boundary and initial conditions. We suggest that many macroscopic discharge properties are properly explained by the model as electromagnetic phenomena, and we discuss extensions of the model to include chemistry, principally ionization and recombination

Simultaneous fitting of statistical-model parameters to symmetric and asymmetric fission cross sections

International Nuclear Information System (INIS)

Mancusi, D; Charity, R J; Cugnon, J

2013-01-01

The de-excitation of compound nuclei has been successfully described for several decades by means of statistical models. However, accurate predictions require some fine-tuning of the model parameters. This task can be simplified by studying several entrance channels, which populate different regions of the parameter space of the compound nucleus. Fusion reactions play an important role in this strategy because they minimise the uncertainty on the entrance channel by fixing mass, charge and excitation energy of the compound nucleus. If incomplete fusion is negligible, the only uncertainty on the compound nucleus comes from the spin distribution. However, some de-excitation channels, such as fission, are quite sensitive to spin. Other entrance channels can then be used to discriminate between equivalent parameter sets. The focus of this work is on fission and intermediate-mass-fragment emission cross sections of compound nuclei with 70 70 ≲ A ≲ 240. 240. The statistical de-excitation model is GEMINI++. The choice of the observables is natural in the framework of GEMINI++, which describes fragment emission using a fissionlike formalism. Equivalent parameter sets for fusion reactions can be resolved using the spallation entrance channel. This promising strategy can lead to the identification of a minimal set of physical ingredients necessary for a unified quantitative description of nuclear de-excitation.

An empirical model describing the postnatal growth of organs in ICRP reference humans: Pt. 1

International Nuclear Information System (INIS)

Walker, J.T.

1991-01-01

An empirical model is presented for describing the postnatal mass growth of lungs in ICRP reference humans. A combined exponential and logistic function containing six parameters is fitted to ICRP 23 lung data using a weighted non-linear least squares technique. The results indicate that the model delineates the data well. Further analysis shows that reference male lungs attain a higher pubertal peak velocity (PPV) and adult mass size than female lungs, although the latter reach their PPV and adult mass size first. Furthermore, the model shows that lung growth rates in infants are two to three orders of magnitude higher than those in mature adults. This finding is important because of the possible association between higher radiation risks in infants' organs that have faster cell turnover rates compared to mature adult organs. The significance of the model for ICRP dosimetric purposes will be discussed. (author)

Parameter Estimation of Partial Differential Equation Models

KAUST Repository

Xun, Xiaolei

2013-09-01

Partial differential equation (PDE) models are commonly used to model complex dynamic systems in applied sciences such as biology and finance. The forms of these PDE models are usually proposed by experts based on their prior knowledge and understanding of the dynamic system. Parameters in PDE models often have interesting scientific interpretations, but their values are often unknown and need to be estimated from the measurements of the dynamic system in the presence of measurement errors. Most PDEs used in practice have no analytic solutions, and can only be solved with numerical methods. Currently, methods for estimating PDE parameters require repeatedly solving PDEs numerically under thousands of candidate parameter values, and thus the computational load is high. In this article, we propose two methods to estimate parameters in PDE models: a parameter cascading method and a Bayesian approach. In both methods, the underlying dynamic process modeled with the PDE model is represented via basis function expansion. For the parameter cascading method, we develop two nested levels of optimization to estimate the PDE parameters. For the Bayesian method, we develop a joint model for data and the PDE and develop a novel hierarchical model allowing us to employ Markov chain Monte Carlo (MCMC) techniques to make posterior inference. Simulation studies show that the Bayesian method and parameter cascading method are comparable, and both outperform other available methods in terms of estimation accuracy. The two methods are demonstrated by estimating parameters in a PDE model from long-range infrared light detection and ranging data. Supplementary materials for this article are available online. © 2013 American Statistical Association.

Parameter Estimation for Thurstone Choice Models

Energy Technology Data Exchange (ETDEWEB)

Vojnovic, Milan [London School of Economics (United Kingdom); Yun, Seyoung [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2017-04-24

We consider the estimation accuracy of individual strength parameters of a Thurstone choice model when each input observation consists of a choice of one item from a set of two or more items (so called top-1 lists). This model accommodates the well-known choice models such as the Luce choice model for comparison sets of two or more items and the Bradley-Terry model for pair comparisons. We provide a tight characterization of the mean squared error of the maximum likelihood parameter estimator. We also provide similar characterizations for parameter estimators defined by a rank-breaking method, which amounts to deducing one or more pair comparisons from a comparison of two or more items, assuming independence of these pair comparisons, and maximizing a likelihood function derived under these assumptions. We also consider a related binary classification problem where each individual parameter takes value from a set of two possible values and the goal is to correctly classify all items within a prescribed classification error. The results of this paper shed light on how the parameter estimation accuracy depends on given Thurstone choice model and the structure of comparison sets. In particular, we found that for unbiased input comparison sets of a given cardinality, when in expectation each comparison set of given cardinality occurs the same number of times, for a broad class of Thurstone choice models, the mean squared error decreases with the cardinality of comparison sets, but only marginally according to a diminishing returns relation. On the other hand, we found that there exist Thurstone choice models for which the mean squared error of the maximum likelihood parameter estimator can decrease much faster with the cardinality of comparison sets. We report empirical evaluation of some claims and key parameters revealed by theory using both synthetic and real-world input data from some popular sport competitions and online labor platforms.

BAYESIAN PARAMETER ESTIMATION IN A MIXED-ORDER MODEL OF BOD DECAY. (U915590)

Science.gov (United States)

We describe a generalized version of the BOD decay model in which the reaction is allowed to assume an order other than one. This is accomplished by making the exponent on BOD concentration a free parameter to be determined by the data. This "mixed-order" model may be ...

Rate-equation modelling and ensemble approach to extraction of parameters for viral infection-induced cell apoptosis and necrosis

Energy Technology Data Exchange (ETDEWEB)

Domanskyi, Sergii; Schilling, Joshua E.; Privman, Vladimir, E-mail: privman@clarkson.edu [Department of Physics, Clarkson University, Potsdam, New York 13676 (United States); Gorshkov, Vyacheslav [National Technical University of Ukraine — KPI, Kiev 03056 (Ukraine); Libert, Sergiy, E-mail: libert@cornell.edu [Department of Biomedical Sciences, Cornell University, Ithaca, New York 14853 (United States)

2016-09-07

We develop a theoretical approach that uses physiochemical kinetics modelling to describe cell population dynamics upon progression of viral infection in cell culture, which results in cell apoptosis (programmed cell death) and necrosis (direct cell death). Several model parameters necessary for computer simulation were determined by reviewing and analyzing available published experimental data. By comparing experimental data to computer modelling results, we identify the parameters that are the most sensitive to the measured system properties and allow for the best data fitting. Our model allows extraction of parameters from experimental data and also has predictive power. Using the model we describe interesting time-dependent quantities that were not directly measured in the experiment and identify correlations among the fitted parameter values. Numerical simulation of viral infection progression is done by a rate-equation approach resulting in a system of “stiff” equations, which are solved by using a novel variant of the stochastic ensemble modelling approach. The latter was originally developed for coupled chemical reactions.

Determining the Walker exponent and developing a modified Smith-Watson-Topper parameter model

Energy Technology Data Exchange (ETDEWEB)

Lv, Zhiqiang; Huang, Hong Zhong; Wang, Hai Kun; Gao, Huiying; Zuo, Fang Jun [University of Electronic Science and Technology of China, Chengdu (China)

2016-03-15

Mean stress effects significantly influence the fatigue life of components. In general, tensile mean stresses are known to reduce the fatigue life of components, whereas compressive mean stresses are known to increase it. To date, various methods that account for mean stress effects have been studied. In this research, considering the high accuracy of mean stress correction and the difficulty in obtaining the material parameter of the Walker method, a practical method is proposed to describe the material parameter of this method. The test data of various materials are then used to verify the proposed practical method. Furthermore, by applying the Walker material parameter and the Smith-Watson-Topper (SWT) parameter, a modified strain-life model is developed to consider sensitivity to mean stress of materials. In addition, three sets of experimental fatigue data from super alloy GH4133, aluminum alloy 7075-T651, and carbon steel are used to estimate the accuracy of the proposed model. A comparison is also made between the SWT parameter method and the proposed strainlife model. The proposed strain-life model provides more accurate life prediction results than the SWT parameter method.

Modeling and Parameter Identification Involving 3-Hydroxypropionaldehyde Inhibitory Effects in Glycerol Continuous Fermentation

Directory of Open Access Journals (Sweden)

Zhaohua Gong

2012-01-01

Full Text Available Mathematical modeling and parameter estimation are critical steps in the optimization of biotechnological processes. In the 1,3-propanediol (1,3-PD production by glycerol fermentation process under anaerobic conditions, 3-hydroxypropionaldehyde (3-HPA accumulation would arouse an irreversible cessation of the fermentation process. Considering 3-HPA inhibitions to cells growth and to activities of enzymes, we propose a novel mathematical model to describe glycerol continuous cultures. Some properties of the above model are discussed. On the basis of the concentrations of extracellular substances, a parameter identification model is established to determine the kinetic parameters in the presented system. Through the penalty function technique combined with an extension of the state space method, an improved genetic algorithm is then constructed to solve the parameter identification model. An illustrative numerical example shows the appropriateness of the proposed model and the validity of optimization algorithm. Since it is difficult to measure the concentrations of intracellular substances, a quantitative robustness analysis method is given to infer whether the model is plausible for the intracellular substances. Numerical results show that the proposed model is of good robustness.

A model describing intra-granular fission gas behaviour in oxide fuel for advanced engineering tools

Science.gov (United States)

Pizzocri, D.; Pastore, G.; Barani, T.; Magni, A.; Luzzi, L.; Van Uffelen, P.; Pitts, S. A.; Alfonsi, A.; Hales, J. D.

2018-04-01

The description of intra-granular fission gas behaviour is a fundamental part of any model for the prediction of fission gas release and swelling in nuclear fuel. In this work we present a model describing the evolution of intra-granular fission gas bubbles in terms of bubble number density and average size, coupled to gas release to grain boundaries. The model considers the fundamental processes of single gas atom diffusion, gas bubble nucleation, re-solution and gas atom trapping at bubbles. The model is derived from a detailed cluster dynamics formulation, yet it consists of only three differential equations in its final form; hence, it can be efficiently applied in engineering fuel performance codes while retaining a physical basis. We discuss improvements relative to previous single-size models for intra-granular bubble evolution. We validate the model against experimental data, both in terms of bubble number density and average bubble radius. Lastly, we perform an uncertainty and sensitivity analysis by propagating the uncertainties in the parameters to model results.

A flowing plasma model to describe drift waves in a cylindrical helicon discharge

International Nuclear Information System (INIS)

Chang, L.; Hole, M. J.; Corr, C. S.

2011-01-01

A two-fluid model developed originally to describe wave oscillations in the vacuum arc centrifuge, a cylindrical, rapidly rotating, low temperature, and confined plasma column, is applied to interpret plasma oscillations in a RF generated linear magnetized plasma [WOMBAT (waves on magnetized beams and turbulence)], with similar density and field strength. Compared to typical centrifuge plasmas, WOMBAT plasmas have slower normalized rotation frequency, lower temperature, and lower axial velocity. Despite these differences, the two-fluid model provides a consistent description of the WOMBAT plasma configuration and yields qualitative agreement between measured and predicted wave oscillation frequencies with axial field strength. In addition, the radial profile of the density perturbation predicted by this model is consistent with the data. Parameter scans show that the dispersion curve is sensitive to the axial field strength and the electron temperature, and the dependence of oscillation frequency with electron temperature matches the experiment. These results consolidate earlier claims that the density and floating potential oscillations are a resistive drift mode, driven by the density gradient. To our knowledge, this is the first detailed physics model of flowing plasmas in the diffusion region away from the RF source. Possible extensions to the model, including temperature nonuniformity and magnetic field oscillations, are also discussed.

Emulation of a complex global aerosol model to quantify sensitivity to uncertain parameters

Directory of Open Access Journals (Sweden)

L. A. Lee

2011-12-01

Full Text Available Sensitivity analysis of atmospheric models is necessary to identify the processes that lead to uncertainty in model predictions, to help understand model diversity through comparison of driving processes, and to prioritise research. Assessing the effect of parameter uncertainty in complex models is challenging and often limited by CPU constraints. Here we present a cost-effective application of variance-based sensitivity analysis to quantify the sensitivity of a 3-D global aerosol model to uncertain parameters. A Gaussian process emulator is used to estimate the model output across multi-dimensional parameter space, using information from a small number of model runs at points chosen using a Latin hypercube space-filling design. Gaussian process emulation is a Bayesian approach that uses information from the model runs along with some prior assumptions about the model behaviour to predict model output everywhere in the uncertainty space. We use the Gaussian process emulator to calculate the percentage of expected output variance explained by uncertainty in global aerosol model parameters and their interactions. To demonstrate the technique, we show examples of cloud condensation nuclei (CCN sensitivity to 8 model parameters in polluted and remote marine environments as a function of altitude. In the polluted environment 95 % of the variance of CCN concentration is described by uncertainty in the 8 parameters (excluding their interaction effects and is dominated by the uncertainty in the sulphur emissions, which explains 80 % of the variance. However, in the remote region parameter interaction effects become important, accounting for up to 40 % of the total variance. Some parameters are shown to have a negligible individual effect but a substantial interaction effect. Such sensitivities would not be detected in the commonly used single parameter perturbation experiments, which would therefore underpredict total uncertainty. Gaussian process

A standard protocol for describing individual-based and agent-based models

Science.gov (United States)

Grimm, Volker; Berger, Uta; Bastiansen, Finn; Eliassen, Sigrunn; Ginot, Vincent; Giske, Jarl; Goss-Custard, John; Grand, Tamara; Heinz, Simone K.; Huse, Geir; Huth, Andreas; Jepsen, Jane U.; Jorgensen, Christian; Mooij, Wolf M.; Muller, Birgit; Pe'er, Guy; Piou, Cyril; Railsback, Steven F.; Robbins, Andrew M.; Robbins, Martha M.; Rossmanith, Eva; Ruger, Nadja; Strand, Espen; Souissi, Sami; Stillman, Richard A.; Vabo, Rune; Visser, Ute; DeAngelis, Donald L.

2006-01-01

Simulation models that describe autonomous individual organisms (individual based models, IBM) or agents (agent-based models, ABM) have become a widely used tool, not only in ecology, but also in many other disciplines dealing with complex systems made up of autonomous entities. However, there is no standard protocol for describing such simulation models, which can make them difficult to understand and to duplicate. This paper presents a proposed standard protocol, ODD, for describing IBMs and ABMs, developed and tested by 28 modellers who cover a wide range of fields within ecology. This protocol consists of three blocks (Overview, Design concepts, and Details), which are subdivided into seven elements: Purpose, State variables and scales, Process overview and scheduling, Design concepts, Initialization, Input, and Submodels. We explain which aspects of a model should be described in each element, and we present an example to illustrate the protocol in use. In addition, 19 examples are available in an Online Appendix. We consider ODD as a first step for establishing a more detailed common format of the description of IBMs and ABMs. Once initiated, the protocol will hopefully evolve as it becomes used by a sufficiently large proportion of modellers.

Using multistage models to describe radiation-induced leukaemia

International Nuclear Information System (INIS)

Little, M.P.; Muirhead, C.R.; Boice, J.D. Jr.; Kleinerman, R.A.

1995-01-01

The Armitage-Doll model of carcinogenesis is fitted to data on leukaemia mortality among the Japanese atomic bomb survivors with the DS86 dosimetry and on leukaemia incidence in the International Radiation Study of Cervical Cancer patients. Two different forms of model are fitted: the first postulates up to two radiation-affected stages and the second additionally allows for the presence at birth of a non-trivial population of cells which have already accumulated the first of the mutations leading to malignancy. Among models of the first form, a model with two adjacent radiation-affected stages appears to fit the data better than other models of the first form, including both models with two affected stages in any order and models with only one affected stage. The best fitting model predicts a linear-quadratic dose-response and reductions of relative risk with increasing time after exposure and age at exposure, in agreement with what has previously been observed in the Japanese and cervical cancer data. However, on the whole it does not provide an adequate fit to either dataset. The second form of model appears to provide a rather better fit, but the optimal models have biologically implausible parameters (the number of initiated cells at birth is negative) so that this model must also be regarded as providing an unsatisfactory description of the data. (author)

DaMoScope and its internet graphics for the visual control of adjusting mathematical models describing experimental data

Science.gov (United States)

Belousov, V. I.; Ezhela, V. V.; Kuyanov, Yu. V.; Tkachenko, N. P.

2015-12-01

The experience of using the dynamic atlas of the experimental data and mathematical models of their description in the problems of adjusting parametric models of observable values depending on kinematic variables is presented. The functional possibilities of an image of a large number of experimental data and the models describing them are shown by examples of data and models of observable values determined by the amplitudes of elastic scattering of hadrons. The Internet implementation of an interactive tool DaMoScope and its interface with the experimental data and codes of adjusted parametric models with the parameters of the best description of data are schematically shown. The DaMoScope codes are freely available.

DaMoScope and its internet graphics for the visual control of adjusting mathematical models describing experimental data

Energy Technology Data Exchange (ETDEWEB)

Belousov, V. I.; Ezhela, V. V.; Kuyanov, Yu. V., E-mail: Yu.Kuyanov@gmail.com; Tkachenko, N. P. [Institute for High Energy Physics, National Research Center Kurchatov Institute, COMPAS Group (Russian Federation)

2015-12-15

The experience of using the dynamic atlas of the experimental data and mathematical models of their description in the problems of adjusting parametric models of observable values depending on kinematic variables is presented. The functional possibilities of an image of a large number of experimental data and the models describing them are shown by examples of data and models of observable values determined by the amplitudes of elastic scattering of hadrons. The Internet implementation of an interactive tool DaMoScope and its interface with the experimental data and codes of adjusted parametric models with the parameters of the best description of data are schematically shown. The DaMoScope codes are freely available.

DaMoScope and its internet graphics for the visual control of adjusting mathematical models describing experimental data

International Nuclear Information System (INIS)

Belousov, V. I.; Ezhela, V. V.; Kuyanov, Yu. V.; Tkachenko, N. P.

2015-01-01

The experience of using the dynamic atlas of the experimental data and mathematical models of their description in the problems of adjusting parametric models of observable values depending on kinematic variables is presented. The functional possibilities of an image of a large number of experimental data and the models describing them are shown by examples of data and models of observable values determined by the amplitudes of elastic scattering of hadrons. The Internet implementation of an interactive tool DaMoScope and its interface with the experimental data and codes of adjusted parametric models with the parameters of the best description of data are schematically shown. The DaMoScope codes are freely available

Calibration and validation of a model describing complete autotrophic nitrogen removal in a granular SBR system

DEFF Research Database (Denmark)

Vangsgaard, Anna Katrine; Mutlu, Ayten Gizem; Gernaey, Krist

2013-01-01

BACKGROUND: A validated model describing the nitritation-anammox process in a granular sequencing batch reactor (SBR) system is an important tool for: a) design of future experiments and b) prediction of process performance during optimization, while applying process control, or during system scale......-up. RESULTS: A model was calibrated using a step-wise procedure customized for the specific needs of the system. The important steps in the procedure were initialization, steady-state and dynamic calibration, and validation. A fast and effective initialization approach was developed to approximate pseudo...... screening of the parameter space proposed by Sin et al. (2008) - to find the best fit of the model to dynamic data. Finally, the calibrated model was validated with an independent data set. CONCLUSION: The presented calibration procedure is the first customized procedure for this type of system...

The influence of model parameters on catchment-response

International Nuclear Information System (INIS)

Shah, S.M.S.; Gabriel, H.F.; Khan, A.A.

2002-01-01

This paper deals with the study of influence of influence of conceptual rainfall-runoff model parameters on catchment response (runoff). A conceptual modified watershed yield model is employed to study the effects of model-parameters on catchment-response, i.e. runoff. The model is calibrated, using manual parameter-fitting approach, also known as trial and error parameter-fitting. In all, there are twenty one (21) parameters that control the functioning of the model. A lumped parametric approach is used. The detailed analysis was performed on Ling River near Kahuta, having catchment area of 56 sq. miles. The model includes physical parameters like GWSM, PETS, PGWRO, etc. fitting coefficients like CINF, CGWS, etc. and initial estimates of the surface-water and groundwater storages i.e. srosp and gwsp. Sensitivity analysis offers a good way, without repetititious computations, the proper weight and consideration that must be taken when each of the influencing factor is evaluated. Sensitivity-analysis was performed to evaluate the influence of model-parameters on runoff. The sensitivity and relative contributions of model parameters influencing catchment-response are studied. (author)

SU-E-T-399: Determination of the Radiobiological Parameters That Describe the Dose-Response Relations of Xerostomia and Disgeusia From Head and Neck Radiotherapy

Energy Technology Data Exchange (ETDEWEB)

Mavroidis, P; Stathakis, S; Papanikolaou, N [University of Texas Health Science Center, UTHSCSA, San Antonio, TX (United States); Peixoto Xavier, C [University of Coimbra, Coimbra, Coimbra (Portugal); Costa Ferreira, B [University of Aveiro, Coimbra, Coimbra (Portugal); Khouri, L; Carmo Lopes, M do [IPOCFG, EPE, Coimbra, Coimbra (Portugal)

2014-06-01

Purpose: To estimate the radiobiological parameters that describe the doseresponse relations of xerostomia and disgeusia from head and neck cancer radiotherapy. To identify the organs that are best correlated with the manifestation of those clinical endpoints. Finally, to evaluate the goodnessof- fit by comparing the model predictions against the actual clinical results. Methods: In this study, 349 head and neck cancer patients were included. For each patient the dose volume histograms (DVH) of parotids (separate and combined), mandible, submandibular glands (separate and combined) and salivary glands were calculated. The follow-up of those patients was recorded at different times after the completion of the treatment (7 weeks, 3, 7, 12, 18 and 24 months). Acute and late xerostomia and acute disgeusia were the clinical endpoints examined. A maximum likelihood fitting was performed to calculate the best estimates of the parameters used by the relative seriality model. The statistical methods of the error distribution, the receiver operating characteristic (ROC) curve, the Pearson's test and the Akaike's information criterion were utilized to assess the goodness-of-fit and the agreement between the pattern of the radiobiological predictions with that of the clinical records. Results: The estimated values of the radiobiological parameters of salivary glands are D50 = 25.2 Gy, γ = 0.52, s = 0.001. The statistical analysis confirmed the clinical validity of those parameters (area under the ROC curve = 0.65 and AIC = 38.3). Conclusion: The analysis proved that the treatment outcome pattern of the patient material can be reproduced by the relative seriality model and the estimated radiobiological parameters. Salivary glands were found to have strong volume dependence (low relative seriality). Diminishing the biologically effective uniform dose to salivary glands below 30 Gy may significantly reduce the risk of complications to the patients irradiated for

SU-E-T-399: Determination of the Radiobiological Parameters That Describe the Dose-Response Relations of Xerostomia and Disgeusia From Head and Neck Radiotherapy

International Nuclear Information System (INIS)

Mavroidis, P; Stathakis, S; Papanikolaou, N; Peixoto Xavier, C; Costa Ferreira, B; Khouri, L; Carmo Lopes, M do

2014-01-01

Purpose: To estimate the radiobiological parameters that describe the doseresponse relations of xerostomia and disgeusia from head and neck cancer radiotherapy. To identify the organs that are best correlated with the manifestation of those clinical endpoints. Finally, to evaluate the goodnessof- fit by comparing the model predictions against the actual clinical results. Methods: In this study, 349 head and neck cancer patients were included. For each patient the dose volume histograms (DVH) of parotids (separate and combined), mandible, submandibular glands (separate and combined) and salivary glands were calculated. The follow-up of those patients was recorded at different times after the completion of the treatment (7 weeks, 3, 7, 12, 18 and 24 months). Acute and late xerostomia and acute disgeusia were the clinical endpoints examined. A maximum likelihood fitting was performed to calculate the best estimates of the parameters used by the relative seriality model. The statistical methods of the error distribution, the receiver operating characteristic (ROC) curve, the Pearson's test and the Akaike's information criterion were utilized to assess the goodness-of-fit and the agreement between the pattern of the radiobiological predictions with that of the clinical records. Results: The estimated values of the radiobiological parameters of salivary glands are D50 = 25.2 Gy, γ = 0.52, s = 0.001. The statistical analysis confirmed the clinical validity of those parameters (area under the ROC curve = 0.65 and AIC = 38.3). Conclusion: The analysis proved that the treatment outcome pattern of the patient material can be reproduced by the relative seriality model and the estimated radiobiological parameters. Salivary glands were found to have strong volume dependence (low relative seriality). Diminishing the biologically effective uniform dose to salivary glands below 30 Gy may significantly reduce the risk of complications to the patients irradiated for

A mathematical model describing the glycemic response of diabetic patients to meal and i.v. infusion of insulin.

Science.gov (United States)

Fabietti, P G; Calabrese, G; Iorio, M; Bistoni, S; Brunetti, P; Sarti, E; Benedetti, M M

2001-10-01

Nine type 1 diabetic patients were studied for 24 hours. During this period they were given three calibrated meals. The glycemia was feedback-controlled by means of an artificial pancreas. The blood concentration of glucose and the infusion speed of the insulin were measured every minute. The experimental data referring to each of the three meals were used to estimate the parameters of a mathematical model suitable for describing the glycemic response of diabetic patients at meals and at the i.v. infusion of exogenous insulin. From the estimate a marked dispersion of the parameters was found, both interindividual and intraindividual. Nevertheless the models thus obtained seem to be usable for the synthesis of a feedback controller, especially in view of creating a portable artificial pancreas that now seems possible owing to the realization (so far experimental) of sufficiently reliable glucose concentration sensors.

On parameter estimation in deformable models

DEFF Research Database (Denmark)

Fisker, Rune; Carstensen, Jens Michael

1998-01-01

Deformable templates have been intensively studied in image analysis through the last decade, but despite its significance the estimation of model parameters has received little attention. We present a method for supervised and unsupervised model parameter estimation using a general Bayesian form...

THREE-PHASE TRANSFORMER PARAMETERS CALCULATION CONSIDERING THE CORE SATURATION FOR THE MATLAB-SIMULINK TRANSFORMER MODEL

Directory of Open Access Journals (Sweden)

I. V. Novash

2015-01-01

Full Text Available This article describes the parameters calculation for the three-phase two-winding power transformer model taken from the SimPowerSystems library, which is the part of the MatLab- Simulink environment. Presented methodology is based on the power transformer nameplate data usage. Particular attention is paid to the power transformer magnetization curve para- meters  calculation.  The  methodology  of  the  three-phase  two-winding  power  transformer model parameters calculation considering the magnetization curve nonlinearity isn’t presented in Russian-and English-language sources. Power transformers demo models described in the SimPowerSystems user’s guide have already calculated parameters, but without reference to the sources of their determination. A power transformer is a nonlinear element of the power system, that’s why for its performance analysis in different modes of operation is necessary to have the magnetization curve parameters.The process analysis during no-load energizing of the power transformer is of special interest. This regime is accompanied by the inrush current on the supply side of the power transformer, which is several times larger than the transformer rated current. Sharp rising of the magnetizing current is explained by the magnetic core saturation. Therefore, magnetiza- tion characteristic accounting during transformer no-load energizing modeling is a mandatory requirement. Article authors attempt to put all calculating formulas in a more convenient form and validate the power transformer nonlinear magnetization characteristics parameters calcu- lation. Inrush current oscillograms obtained during the simulation experiment confirmed the adequacy of the calculated model parameters.

Assessment of parameter regionalization methods for modeling flash floods in China

Science.gov (United States)

Ragettli, Silvan; Zhou, Jian; Wang, Haijing

2017-04-01

catchments resulted in good model performance (NSE > 0.5) in 10 and medium performance (NSE > 0.2) in 3 catchments. Optimal model parameters proofed to be relatively insensitive to different HRU configurations. This suggests that dominant controls on hydrologic parameter transfer can potentially be identified based on catchment attributes describing meteorological, geological or landscape characteristics. Parameter regionalization based on a principal component analysis (PCA) nearest neighbor search (using all available catchment attributes) resulted in a 54% success rate in transferring optimal parameter sets and still yielding acceptable model performance. Data from more catchments are required to further increase the parameter transferability success rate or to develop regionalization strategies for individual parameters.

Parameter sensitivity and identifiability for a biogeochemical model of hypoxia in the northern Gulf of Mexico

Science.gov (United States)

Local sensitivity analyses and identifiable parameter subsets were used to describe numerical constraints of a hypoxia model for bottom waters of the northern Gulf of Mexico. The sensitivity of state variables differed considerably with parameter changes, although most variables ...

Digital Materials – Evaluation of the Possibilities of using Selected Hyperelastic Models to Describe Constitutive Relations

Directory of Open Access Journals (Sweden)

Mańkowski J.

2017-08-01

Full Text Available The authors tried to identify the parameters of numerical models of digital materials, which are a kind of composite resulting from the manufacture of the product in 3D printers. With the arrangement of several heads of the printer, the new material can result from mixing of materials with radically different properties, during the process of producing single layer of the product. The new material has properties dependent on the base materials properties and their proportions. Digital materials tensile characteristics are often non-linear and qualify to be described by hyperelastic materials models. The identification was conducted based on the results of tensile tests models, its various degrees coefficients of the polynomials to various degrees coefficients of the polynomials. The Drucker’s stability criterion was also examined. Fourteen different materials were analyzed.

HERMES: A Model to Describe Deformation, Burning, Explosion, and Detonation

Energy Technology Data Exchange (ETDEWEB)

Reaugh, J E

2011-11-22

HERMES (High Explosive Response to MEchanical Stimulus) was developed to fill the need for a model to describe an explosive response of the type described as BVR (Burn to Violent Response) or HEVR (High Explosive Violent Response). Characteristically this response leaves a substantial amount of explosive unconsumed, the time to reaction is long, and the peak pressure developed is low. In contrast, detonations characteristically consume all explosive present, the time to reaction is short, and peak pressures are high. However, most of the previous models to describe explosive response were models for detonation. The earliest models to describe the response of explosives to mechanical stimulus in computer simulations were applied to intentional detonation (performance) of nearly ideal explosives. In this case, an ideal explosive is one with a vanishingly small reaction zone. A detonation is supersonic with respect to the undetonated explosive (reactant). The reactant cannot respond to the pressure of the detonation before the detonation front arrives, so the precise compressibility of the reactant does not matter. Further, the mesh sizes that were practical for the computer resources then available were large with respect to the reaction zone. As a result, methods then used to model detonations, known as {beta}-burn or program burn, were not intended to resolve the structure of the reaction zone. Instead, these methods spread the detonation front over a few finite-difference zones, in the same spirit that artificial viscosity is used to spread the shock front in inert materials over a few finite-difference zones. These methods are still widely used when the structure of the reaction zone and the build-up to detonation are unimportant. Later detonation models resolved the reaction zone. These models were applied both to performance, particularly as it is affected by the size of the charge, and to situations in which the stimulus was less than that needed for reliable

Coast-down model based on rated parameters of reactor coolant pump

International Nuclear Information System (INIS)

Jiang Maohua; Zou Zhichao; Wang Pengfei; Ruan Xiaodong

2014-01-01

For a sudden loss of power in reactor coolant pump (RCP), a calculation model of rotor speed and flow characteristics based on rated parameters was studied. The derived model was verified by comparing with the power-off experimental data of 100D RCP. The results indicate that it can be used in preliminary design calculation and verification analysis. Then a design criterion of RCP was described based on the calculation model. The moment of inertia in AP1000 RCP was verified by this criterion. (authors)

Complexity, parameter sensitivity and parameter transferability in the modelling of floodplain inundation

Science.gov (United States)

Bates, P. D.; Neal, J. C.; Fewtrell, T. J.

2012-12-01

In this we paper we consider two related questions. First, we address the issue of how much physical complexity is necessary in a model in order to simulate floodplain inundation to within validation data error. This is achieved through development of a single code/multiple physics hydraulic model (LISFLOOD-FP) where different degrees of complexity can be switched on or off. Different configurations of this code are applied to four benchmark test cases, and compared to the results of a number of industry standard models. Second we address the issue of how parameter sensitivity and transferability change with increasing complexity using numerical experiments with models of different physical and geometric intricacy. Hydraulic models are a good example system with which to address such generic modelling questions as: (1) they have a strong physical basis; (2) there is only one set of equations to solve; (3) they require only topography and boundary conditions as input data; and (4) they typically require only a single free parameter, namely boundary friction. In terms of complexity required we show that for the problem of sub-critical floodplain inundation a number of codes of different dimensionality and resolution can be found to fit uncertain model validation data equally well, and that in this situation Occam's razor emerges as a useful logic to guide model selection. We find also find that model skill usually improves more rapidly with increases in model spatial resolution than increases in physical complexity, and that standard approaches to testing hydraulic models against laboratory data or analytical solutions may fail to identify this important fact. Lastly, we find that in benchmark testing studies significant differences can exist between codes with identical numerical solution techniques as a result of auxiliary choices regarding the specifics of model implementation that are frequently unreported by code developers. As a consequence, making sound

Multi-Scale Parameter Identification of Lithium-Ion Battery Electric Models Using a PSO-LM Algorithm

Directory of Open Access Journals (Sweden)

Wen-Jing Shen

2017-03-01

Full Text Available This paper proposes a multi-scale parameter identification algorithm for the lithium-ion battery (LIB electric model by using a combination of particle swarm optimization (PSO and Levenberg-Marquardt (LM algorithms. Two-dimensional Poisson equations with unknown parameters are used to describe the potential and current density distribution (PDD of the positive and negative electrodes in the LIB electric model. The model parameters are difficult to determine in the simulation due to the nonlinear complexity of the model. In the proposed identification algorithm, PSO is used for the coarse-scale parameter identification and the LM algorithm is applied for the fine-scale parameter identification. The experiment results show that the multi-scale identification not only improves the convergence rate and effectively escapes from the stagnation of PSO, but also overcomes the local minimum entrapment drawback of the LM algorithm. The terminal voltage curves from the PDD model with the identified parameter values are in good agreement with those from the experiments at different discharge/charge rates.

Dynamic Friction Parameter Identification Method with LuGre Model for Direct-Drive Rotary Torque Motor

Directory of Open Access Journals (Sweden)

Xingjian Wang

2016-01-01

Full Text Available Attainment of high-performance motion/velocity control objectives for the Direct-Drive Rotary (DDR torque motor should fully consider practical nonlinearities in controller design, such as dynamic friction. The LuGre model has been widely utilized to describe nonlinear friction behavior; however, parameter identification for the LuGre model remains a challenge. A new dynamic friction parameter identification method for LuGre model is proposed in this study. Static parameters are identified through a series of constant velocity experiments, while dynamic parameters are obtained through a presliding process. Novel evolutionary algorithm (NEA is utilized to increase identification accuracy. Experimental results gathered from the identification experiments conducted in the study for a practical DDR torque motor control system validate the effectiveness of the proposed method.

New Values of Cross-Talk Parameters for Twisted Pair Model

OpenAIRE

Milos Kozak; Lukas Cepa; Jiri Vodrazka

2010-01-01

Near-end Crosstalk (NEXT) and Far-end Crosstalk (FEXT) of unshielded twisted pair (UTP) cable are the main factors limiting the information capacity in data transmission. Crosstalk depends mostly on the frequency. Frequency dependent transfer functions and crosstalk attenuation may be obtained by measurement, but for the analytical description of the transmission channel's parameters is useful to define functions modelling the crosstalk. The paper describes the measuri...

Inverse modeling approach for evaluation of kinetic parameters of a biofilm reactor using tabu search.

Science.gov (United States)

Kumar, B Shiva; Venkateswarlu, Ch

2014-08-01

The complex nature of biological reactions in biofilm reactors often poses difficulties in analyzing such reactors experimentally. Mathematical models could be very useful for their design and analysis. However, application of biofilm reactor models to practical problems proves somewhat ineffective due to the lack of knowledge of accurate kinetic models and uncertainty in model parameters. In this work, we propose an inverse modeling approach based on tabu search (TS) to estimate the parameters of kinetic and film thickness models. TS is used to estimate these parameters as a consequence of the validation of the mathematical models of the process with the aid of measured data obtained from an experimental fixed-bed anaerobic biofilm reactor involving the treatment of pharmaceutical industry wastewater. The results evaluated for different modeling configurations of varying degrees of complexity illustrate the effectiveness of TS for accurate estimation of kinetic and film thickness model parameters of the biofilm process. The results show that the two-dimensional mathematical model with Edward kinetics (with its optimum parameters as mu(max)rho(s)/Y = 24.57, Ks = 1.352 and Ki = 102.36) and three-parameter film thickness expression (with its estimated parameters as a = 0.289 x 10(-5), b = 1.55 x 10(-4) and c = 15.2 x 10(-6)) better describes the biofilm reactor treating the industry wastewater.

Field tests and machine learning approaches for refining algorithms and correlations of driver's model parameters.

Science.gov (United States)

Tango, Fabio; Minin, Luca; Tesauri, Francesco; Montanari, Roberto

2010-03-01

This paper describes the field tests on a driving simulator carried out to validate the algorithms and the correlations of dynamic parameters, specifically driving task demand and drivers' distraction, able to predict drivers' intentions. These parameters belong to the driver's model developed by AIDE (Adaptive Integrated Driver-vehicle InterfacE) European Integrated Project. Drivers' behavioural data have been collected from the simulator tests to model and validate these parameters using machine learning techniques, specifically the adaptive neuro fuzzy inference systems (ANFIS) and the artificial neural network (ANN). Two models of task demand and distraction have been developed, one for each adopted technique. The paper provides an overview of the driver's model, the description of the task demand and distraction modelling and the tests conducted for the validation of these parameters. A test comparing predicted and expected outcomes of the modelled parameters for each machine learning technique has been carried out: for distraction, in particular, promising results (low prediction errors) have been obtained by adopting an artificial neural network.

COMPREHENSIVE CHECK MEASUREMENT OF KEY PARAMETERS ON MODEL BELT CONVEYOR

Directory of Open Access Journals (Sweden)

Vlastimil MONI

2013-07-01

Full Text Available Complex measurements of characteristic parameters realised on a long distance model belt conveyor are described. The main objective was to complete and combine the regular measurements of electric power on drives of belt conveyors operated in Czech opencast mines with measurements of other physical quantities and to gain by this way an image of their mutual relations and relations of quantities derived from them. The paper includes a short description and results of the measurements on an experimental model conveyor with a closed material transport way.

On the role of modeling parameters in IMRT plan optimization

International Nuclear Information System (INIS)

Krause, Michael; Scherrer, Alexander; Thieke, Christian

2008-01-01

The formulation of optimization problems in intensity-modulated radiotherapy (IMRT) planning comprises the choice of various values such as function-specific parameters or constraint bounds. In current inverse planning programs that yield a single treatment plan for each optimization, it is often unclear how strongly these modeling parameters affect the resulting plan. This work investigates the mathematical concepts of elasticity and sensitivity to deal with this problem. An artificial planning case with a horse-shoe formed target with different opening angles surrounding a circular risk structure is studied. As evaluation functions the generalized equivalent uniform dose (EUD) and the average underdosage below and average overdosage beyond certain dose thresholds are used. A single IMRT plan is calculated for an exemplary parameter configuration. The elasticity and sensitivity of each parameter are then calculated without re-optimization, and the results are numerically verified. The results show the following. (1) elasticity can quantify the influence of a modeling parameter on the optimization result in terms of how strongly the objective function value varies under modifications of the parameter value. It also can describe how strongly the geometry of the involved planning structures affects the optimization result. (2) Based on the current parameter settings and corresponding treatment plan, sensitivity analysis can predict the optimization result for modified parameter values without re-optimization, and it can estimate the value intervals in which such predictions are valid. In conclusion, elasticity and sensitivity can provide helpful tools in inverse IMRT planning to identify the most critical parameters of an individual planning problem and to modify their values in an appropriate way

Comparison of two mathematical models for describing heat-induced cell killing

International Nuclear Information System (INIS)

Roti Roti, J.L.; Henle, K.J.

1980-01-01

A computer-based minimization algorithm is utilized to obtain the optimum fits of two models to hyperthermic cell killing data. The models chosen are the multitarget, single-hit equation, which is in general use, and the linear-quadratic equation, which has been applied to cell killing by ionizing irradiation but not to heat-induced cell killing. The linear-quadratic equation fits hyperthermic cell killing data as well as the multitarget, single-hit equation. Both parameters of the linear-quadratic equation obey the Arrhenius law, whereas only one of the two parameters of the multitarget, single-hit equation obeys the Arrhenius law. Thus the linear-quadratic function can completely define cell killing as a function of both time and temperature. In addition, the linear-quadratic model will provide a simplified approach to the study of the synergism between heat and X irradiation

Universally sloppy parameter sensitivities in systems biology models.

Directory of Open Access Journals (Sweden)

Ryan N Gutenkunst

2007-10-01

Full Text Available Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring in vivo biochemical parameters is difficult, and collectively fitting them to other experimental data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a "sloppy" spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our growth-factor-signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.

Universally sloppy parameter sensitivities in systems biology models.

Science.gov (United States)

Gutenkunst, Ryan N; Waterfall, Joshua J; Casey, Fergal P; Brown, Kevin S; Myers, Christopher R; Sethna, James P

2007-10-01

Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring in vivo biochemical parameters is difficult, and collectively fitting them to other experimental data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a "sloppy" spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our growth-factor-signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.

Evaluation of a cross contamination model describing transfer of Salmonella spp. and Listeria monocytogenes during grinding of pork and beef.

Science.gov (United States)

Møller, C O A; Sant'Ana, A S; Hansen, S K H; Nauta, M J; Silva, L P; Alvarenga, V O; Maffei, D; Silva, F F P; Lopes, J T; Franco, B D G M; Aabo, S; Hansen, T B

2016-06-02

In a previous study, a model was developed to describe the transfer and survival of Salmonella during grinding of pork (Møller, C.O.A., Nauta, M.J., Christensen, B.B., Dalgaard, P., Hansen, T.B., 2012. Modelling transfer of Salmonella typhimurium DT104 during simulation of grinding of pork. Journal of Applied Microbiology 112 (1), 90-98). The robustness of this model is now evaluated by studying its performance for predicting the transfer and survival of Salmonella spp. and Listeria monocytogenes during grinding of different types of meat (pork and beef), using two different grinders, different sizes and different numbers of pieces of meats to be ground. A total of 19 grinding trials were collected. Acceptable Simulation Zone (ASZ), visual inspection of the data, Quantitative Microbiological Risk Assessment (QMRA), as well as the Total Transfer Potential (TTP) were used as approaches to evaluate model performance and to access the quality of the cross contamination model predictions. Using the ASZ approach and considering that 70% of the observed counts have to be inside a defined acceptable zone of ±0.5 log10CFU per portion, it was found that the cross contamination parameters suggested by Møller et al. (2012) were not able to describe all 19 trials. However, for each of the collected grinding trials, the transfer event was well described when fitted to the model structure proposed by Møller et al. (2012). Parameter estimates obtained by fitting observed trials performed at different conditions, such as size and number of pieces of meat to be ground, may not be applied to describe cross contamination of unlike processing. Nevertheless, the risk estimates, as well as the TTP, revealed that the risk of disease may be reduced when the grinding of meat is performed in a grinder made of stainless steel (for all surfaces in contact with the meat), using a well-sharpened knife and holding at room temperatures lower than 4°C. Copyright © 2016 Elsevier B.V. All

Mathematical modelling of nonstationary processes in a regenerator with dissociating coolant at supercritical parameters

International Nuclear Information System (INIS)

Tashchilova, Eh.M.; Sharovarov, G.A.

1985-01-01

The mathematical model of nonstationary processes in heat exchangers with dissociating coolant at supercritical parameters is given. Its dimensionless criteria are deveped. The effect of NPP regenerator parameters on criteria variation is determined. The proceeding nonstationary processes are estimated qualitatively using the dimensionless parameters. Dynamics of the processes in heat exchangers is described by the energy, mass and moment-of-momentum equations for heating and heated medium taking into account heat accumulation in the heat-transfer wall and distribution of parameters along the length of a heat exchanger

A Parameter-based Model for Generating Culturally Adaptive Nonverbal Behaviors in Embodied Conversational Agents

DEFF Research Database (Denmark)

Lipi, Afia Akhter; Nakano, Yukiko; Rehm, Matthias

2009-01-01

The goal of this paper is to integrate culture as a computational term in embodied conversational agents by employing an empirical data-driven approach as well as a theoretical model-driven approach. We propose a parameter-based model that predicts nonverbal expressions appropriate for specific...... cultures. First, we introduce the Hofstede theory to describe socio-cultural characteristics of each country. Then, based on the previous studies in cultural differences of nonverbal behaviors, we propose expressive parameters to characterize nonverbal behaviors. Finally, by integrating socio-cultural...

Parameters estimation for reactive transport: A way to test the validity of a reactive model

Science.gov (United States)

Aggarwal, Mohit; Cheikh Anta Ndiaye, Mame; Carrayrou, Jérôme

The chemical parameters used in reactive transport models are not known accurately due to the complexity and the heterogeneous conditions of a real domain. We will present an efficient algorithm in order to estimate the chemical parameters using Monte-Carlo method. Monte-Carlo methods are very robust for the optimisation of the highly non-linear mathematical model describing reactive transport. Reactive transport of tributyltin (TBT) through natural quartz sand at seven different pHs is taken as the test case. Our algorithm will be used to estimate the chemical parameters of the sorption of TBT onto the natural quartz sand. By testing and comparing three models of surface complexation, we show that the proposed adsorption model cannot explain the experimental data.

Flare parameters inferred from a 3D loop model data base

Science.gov (United States)

Cuambe, Valente A.; Costa, J. E. R.; Simões, P. J. A.

2018-06-01

We developed a data base of pre-calculated flare images and spectra exploring a set of parameters which describe the physical characteristics of coronal loops and accelerated electron distribution. Due to the large number of parameters involved in describing the geometry and the flaring atmosphere in the model used, we built a large data base of models (˜250 000) to facilitate the flare analysis. The geometry and characteristics of non-thermal electrons are defined on a discrete grid with spatial resolution greater than 4 arcsec. The data base was constructed based on general properties of known solar flares and convolved with instrumental resolution to replicate the observations from the Nobeyama radio polarimeter spectra and Nobeyama radioheliograph (NoRH) brightness maps. Observed spectra and brightness distribution maps are easily compared with the modelled spectra and images in the data base, indicating a possible range of solutions. The parameter search efficiency in this finite data base is discussed. 8 out of 10 parameters analysed for 1000 simulated flare searches were recovered with a relative error of less than 20 per cent on average. In addition, from the analysis of the observed correlation between NoRH flare sizes and intensities at 17 GHz, some statistical properties were derived. From these statistics, the energy spectral index was found to be δ ˜ 3, with non-thermal electron densities showing a peak distribution ⪅107 cm-3, and Bphotosphere ⪆ 2000 G. Some bias for larger loops with heights as great as ˜2.6 × 109 cm, and looptop events were noted. An excellent match of the spectrum and the brightness distribution at 17 and 34 GHz of the 2002 May 31 flare is presented as well.

Systematic parameter inference in stochastic mesoscopic modeling

Energy Technology Data Exchange (ETDEWEB)

Lei, Huan; Yang, Xiu [Pacific Northwest National Laboratory, Richland, WA 99352 (United States); Li, Zhen [Division of Applied Mathematics, Brown University, Providence, RI 02912 (United States); Karniadakis, George Em, E-mail: george_karniadakis@brown.edu [Division of Applied Mathematics, Brown University, Providence, RI 02912 (United States)

2017-02-01

We propose a method to efficiently determine the optimal coarse-grained force field in mesoscopic stochastic simulations of Newtonian fluid and polymer melt systems modeled by dissipative particle dynamics (DPD) and energy conserving dissipative particle dynamics (eDPD). The response surfaces of various target properties (viscosity, diffusivity, pressure, etc.) with respect to model parameters are constructed based on the generalized polynomial chaos (gPC) expansion using simulation results on sampling points (e.g., individual parameter sets). To alleviate the computational cost to evaluate the target properties, we employ the compressive sensing method to compute the coefficients of the dominant gPC terms given the prior knowledge that the coefficients are “sparse”. The proposed method shows comparable accuracy with the standard probabilistic collocation method (PCM) while it imposes a much weaker restriction on the number of the simulation samples especially for systems with high dimensional parametric space. Fully access to the response surfaces within the confidence range enables us to infer the optimal force parameters given the desirable values of target properties at the macroscopic scale. Moreover, it enables us to investigate the intrinsic relationship between the model parameters, identify possible degeneracies in the parameter space, and optimize the model by eliminating model redundancies. The proposed method provides an efficient alternative approach for constructing mesoscopic models by inferring model parameters to recover target properties of the physics systems (e.g., from experimental measurements), where those force field parameters and formulation cannot be derived from the microscopic level in a straight forward way.

Prediction Model of Battery State of Charge and Control Parameter Optimization for Electric Vehicle

Directory of Open Access Journals (Sweden)

Bambang Wahono

2015-07-01

Full Text Available This paper presents the construction of a battery state of charge (SOC prediction model and the optimization method of the said model to appropriately control the number of parameters in compliance with the SOC as the battery output objectives. Research Centre for Electrical Power and Mechatronics, Indonesian Institute of Sciences has tested its electric vehicle research prototype on the road, monitoring its voltage, current, temperature, time, vehicle velocity, motor speed, and SOC during the operation. Using this experimental data, the prediction model of battery SOC was built. Stepwise method considering multicollinearity was able to efficiently develops the battery prediction model that describes the multiple control parameters in relation to the characteristic values such as SOC. It was demonstrated that particle swarm optimization (PSO succesfully and efficiently calculated optimal control parameters to optimize evaluation item such as SOC based on the model.

Mathematical modelling of simultaneous solvent evaporation and chemical curing in thermoset coatings: A parameter study

DEFF Research Database (Denmark)

Kiil, Søren

2011-01-01

A mathematical model, describing the curing behaviour of a two-component, solvent-based, thermoset coating, is used to conduct a parameter study. The model includes curing reactions, solvent intra-film diffusion and evaporation, film gelation, vitrification, and crosslinking. A case study with a ...

A Kinetic Model Describing Injury-Burden in Team Sports.

Science.gov (United States)

Fuller, Colin W

2017-12-01

Injuries in team sports are normally characterised by the incidence, severity, and location and type of injuries sustained: these measures, however, do not provide an insight into the variable injury-burden experienced during a season. Injury burden varies according to the team's match and training loads, the rate at which injuries are sustained and the time taken for these injuries to resolve. At the present time, this time-based variation of injury burden has not been modelled. To develop a kinetic model describing the time-based injury burden experienced by teams in elite team sports and to demonstrate the model's utility. Rates of injury were quantified using a large eight-season database of rugby injuries (5253) and exposure (60,085 player-match-hours) in English professional rugby. Rates of recovery from injury were quantified using time-to-recovery analysis of the injuries. The kinetic model proposed for predicting a team's time-based injury burden is based on a composite rate equation developed from the incidence of injury, a first-order rate of recovery from injury and the team's playing load. The utility of the model was demonstrated by examining common scenarios encountered in elite rugby. The kinetic model developed describes and predicts the variable injury-burden arising from match play during a season of rugby union based on the incidence of match injuries, the rate of recovery from injury and the playing load. The model is equally applicable to other team sports and other scenarios.

On the identification of multiple space dependent ionic parameters in cardiac electrophysiology modelling

Science.gov (United States)

Abidi, Yassine; Bellassoued, Mourad; Mahjoub, Moncef; Zemzemi, Nejib

2018-03-01

In this paper, we consider the inverse problem of space dependent multiple ionic parameters identification in cardiac electrophysiology modelling from a set of observations. We use the monodomain system known as a state-of-the-art model in cardiac electrophysiology and we consider a general Hodgkin-Huxley formalism to describe the ionic exchanges at the microscopic level. This formalism covers many physiological transmembrane potential models including those in cardiac electrophysiology. Our main result is the proof of the uniqueness and a Lipschitz stability estimate of ion channels conductance parameters based on some observations on an arbitrary subdomain. The key idea is a Carleman estimate for a parabolic operator with multiple coefficients and an ordinary differential equation system.

Integrated model of insulin and glucose kinetics describing both hepatic glucose and pancreatic insulin regulation

DEFF Research Database (Denmark)

Erlandsen, Mogens; Martinussen, Christoffer; Gravholt, Claus Højbjerg

2018-01-01

AbstractBackground and objectives Modeling of glucose kinetics has to a large extent been based on models with plasma insulin as a known forcing function. Furthermore, population-based statistical methods for parameter estimation in these models have mainly addressed random inter-individual varia......AbstractBackground and objectives Modeling of glucose kinetics has to a large extent been based on models with plasma insulin as a known forcing function. Furthermore, population-based statistical methods for parameter estimation in these models have mainly addressed random inter......-individual variations and not intra-individual variations in the parameters. Here we present an integrated whole-body model of glucose and insulin kinetics which extends the well-known two-compartment glucose minimal model. The population-based estimation technique allow for quantification of both random inter......- and intra-individual variation in selected parameters using simultaneous data series on glucose and insulin. Methods We extend the two-compartment glucose model into a whole-body model for both glucose and insulin using a simple model for the pancreas compartment which includes feedback of glucose on both...

Test models for improving filtering with model errors through stochastic parameter estimation

International Nuclear Information System (INIS)

Gershgorin, B.; Harlim, J.; Majda, A.J.

2010-01-01

The filtering skill for turbulent signals from nature is often limited by model errors created by utilizing an imperfect model for filtering. Updating the parameters in the imperfect model through stochastic parameter estimation is one way to increase filtering skill and model performance. Here a suite of stringent test models for filtering with stochastic parameter estimation is developed based on the Stochastic Parameterization Extended Kalman Filter (SPEKF). These new SPEKF-algorithms systematically correct both multiplicative and additive biases and involve exact formulas for propagating the mean and covariance including the parameters in the test model. A comprehensive study is presented of robust parameter regimes for increasing filtering skill through stochastic parameter estimation for turbulent signals as the observation time and observation noise are varied and even when the forcing is incorrectly specified. The results here provide useful guidelines for filtering turbulent signals in more complex systems with significant model errors.

Evaluation and Validation of a TCAT Model to Describe Non-Dilute Flow and Species Transport in Porous Media

Science.gov (United States)

Weigand, T. M.; Harrison, E.; Miller, C. T.

2017-12-01

A thermodynamically constrained averaging theory (TCAT) model has been developed to simulate non-dilute flow and species transport in porous media. This model has the advantages of a firm connection between the microscale, or pore scale, and the macroscale; a thermodynamically consistent basis; the explicit inclusion of dissipative terms that arise from spatial gradients in pressure and chemical activity; and the ability to describe both high and low concentration displacement. The TCAT model has previously been shown to provide excellent agreement for a set of laboratory data and outperformed existing macroscale models that have been used for non-dilute flow and transport. The examined experimental dataset consisted of stable brine displacements for a large range of fluid properties. This dataset however only examined one type of porous media and had a fixed flow rate for all experiments. In this work, the TCAT model is applied to a dataset that consists of two different porous media types, constant head and flow rate conditions, varying resident fluid concentrations, and internal probes that measured the pressure and salt mass fraction. Parameter estimation is performed on a subset of the experimental data for the TCAT model as well as other existing non-dilute flow and transport models. The optimized parameters are then used for forward simulations and the accuracy of the models is compared.

Site-specific parameter values for the Nuclear Regulatory Commission's food pathway dose model

International Nuclear Information System (INIS)

Hamby, D.M.

1992-01-01

Routine operations at the Savannah River Site (SRS) in Western South Carolina result in radionuclide releases to the atmosphere and to the Savannah River. The resulting radiation doses to the off-site maximum individual and the off-site population within 80 km of the SRS are estimated on a yearly basis. These estimates are currently generated using dose models prescribed for the commercial nuclear power industry by the Nuclear Regulatory Commission (NRC). The NRC provides default values for dose-model parameters for facilities without resources to develop site-specific values. A survey of land- and water-use characteristics for the Savannah River area has been conducted to determine site-specific values for water recreation, consumption, and agricultural parameters used in the NRC Regulatory Guide 1.109 (1977) dosimetric models. These site parameters include local characteristics of meat, milk, and vegetable production; recreational and commercial activities on the Savannah River; and meat, milk, vegetable, and seafood consumption rates. This paper describes how parameter data were obtained at the Savannah River Site and the impacts of such data on off-site dose. Dose estimates using site-specific parameter values are compared to estimates using the NRC default values

Study on Identification of Material Model Parameters from Compact Tension Test on Concrete Specimens

Science.gov (United States)

Hokes, Filip; Kral, Petr; Husek, Martin; Kala, Jiri

2017-10-01

Identification of a concrete material model parameters using optimization is based on a calculation of a difference between experimentally measured and numerically obtained data. Measure of the difference can be formulated via root mean squared error that is often used for determination of accuracy of a mathematical model in the field of meteorology or demography. The quality of the identified parameters is, however, determined not only by right choice of an objective function but also by the source experimental data. One of the possible way is to use load-displacement curves from three-point bending tests that were performed on concrete specimens. This option shows the significance of modulus of elasticity, tensile strength and specific fracture energy. Another possible option is to use experimental data from compact tension test. It is clear that the response in the second type of test is also dependent on the above mentioned material parameters. The question is whether the parameters identified within three-point bending test and within compact tension test will reach the same values. The presented article brings the numerical study of inverse identification of material model parameters from experimental data measured during compact tension tests. The article also presents utilization of the modified sensitivity analysis that calculates the sensitivity of the material model parameters for different parts of loading curve. The main goal of the article is to describe the process of inverse identification of parameters for plasticity-based material model of concrete and prepare data for future comparison with identified values of the material model parameters from different type of fracture tests.

Standard model parameters and the search for new physics

International Nuclear Information System (INIS)

Marciano, W.J.

1988-04-01

In these lectures, my aim is to present an up-to-date status report on the standard model and some key tests of electroweak unification. Within that context, I also discuss how and where hints of new physics may emerge. To accomplish those goals, I have organized my presentation as follows: I discuss the standard model parameters with particular emphasis on the gauge coupling constants and vector boson masses. Examples of new physics appendages are also briefly commented on. In addition, because these lectures are intended for students and thus somewhat pedagogical, I have included an appendix on dimensional regularization and a simple computational example that employs that technique. Next, I focus on weak charged current phenomenology. Precision tests of the standard model are described and up-to-date values for the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix parameters are presented. Constraints implied by those tests for a 4th generation, supersymmetry, extra Z/prime/ bosons, and compositeness are also discussed. I discuss weak neutral current phenomenology and the extraction of sin/sup 2/ /theta//sub W/ from experiment. The results presented there are based on a recently completed global analysis of all existing data. I have chosen to concentrate that discussion on radiative corrections, the effect of a heavy top quark mass, and implications for grand unified theories (GUTS). The potential for further experimental progress is also commented on. I depart from the narrowest version of the standard model and discuss effects of neutrino masses and mixings. I have chosen to concentrate on oscillations, the Mikheyev-Smirnov- Wolfenstein (MSW) effect, and electromagnetic properties of neutrinos. On the latter topic, I will describe some recent work on resonant spin-flavor precession. Finally, I conclude with a prospectus on hopes for the future. 76 refs

Exploiting intrinsic fluctuations to identify model parameters.

Science.gov (United States)

Zimmer, Christoph; Sahle, Sven; Pahle, Jürgen

2015-04-01

Parameterisation of kinetic models plays a central role in computational systems biology. Besides the lack of experimental data of high enough quality, some of the biggest challenges here are identification issues. Model parameters can be structurally non-identifiable because of functional relationships. Noise in measured data is usually considered to be a nuisance for parameter estimation. However, it turns out that intrinsic fluctuations in particle numbers can make parameters identifiable that were previously non-identifiable. The authors present a method to identify model parameters that are structurally non-identifiable in a deterministic framework. The method takes time course recordings of biochemical systems in steady state or transient state as input. Often a functional relationship between parameters presents itself by a one-dimensional manifold in parameter space containing parameter sets of optimal goodness. Although the system's behaviour cannot be distinguished on this manifold in a deterministic framework it might be distinguishable in a stochastic modelling framework. Their method exploits this by using an objective function that includes a measure for fluctuations in particle numbers. They show on three example models, immigration-death, gene expression and Epo-EpoReceptor interaction, that this resolves the non-identifiability even in the case of measurement noise with known amplitude. The method is applied to partially observed recordings of biochemical systems with measurement noise. It is simple to implement and it is usually very fast to compute. This optimisation can be realised in a classical or Bayesian fashion.

Assessing composition and structure of soft biphasic media from Kelvin-Voigt fractional derivative model parameters

Science.gov (United States)

Zhang, Hongmei; Wang, Yue; Fatemi, Mostafa; Insana, Michael F.

2017-03-01

Kelvin-Voigt fractional derivative (KVFD) model parameters have been used to describe viscoelastic properties of soft tissues. However, translating model parameters into a concise set of intrinsic mechanical properties related to tissue composition and structure remains challenging. This paper begins by exploring these relationships using a biphasic emulsion materials with known composition. Mechanical properties are measured by analyzing data from two indentation techniques—ramp-stress relaxation and load-unload hysteresis tests. Material composition is predictably correlated with viscoelastic model parameters. Model parameters estimated from the tests reveal that elastic modulus E 0 closely approximates the shear modulus for pure gelatin. Fractional-order parameter α and time constant τ vary monotonically with the volume fraction of the material’s fluid component. α characterizes medium fluidity and the rate of energy dissipation, and τ is a viscous time constant. Numerical simulations suggest that the viscous coefficient η is proportional to the energy lost during quasi-static force-displacement cycles, E A . The slope of E A versus η is determined by α and the applied indentation ramp time T r. Experimental measurements from phantom and ex vivo liver data show close agreement with theoretical predictions of the η -{{E}A} relation. The relative error is less than 20% for emulsions 22% for liver. We find that KVFD model parameters form a concise features space for biphasic medium characterization that described time-varying mechanical properties. The experimental work was carried out at the Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA. Methodological development, including numerical simulation and all data analysis, were carried out at the school of Life Science and Technology, Xi’an JiaoTong University, 710049, China.

Accounting for parameter uncertainty in the definition of parametric distributions used to describe individual patient variation in health economic models

NARCIS (Netherlands)

Degeling, Koen; Ijzerman, Maarten J.; Koopman, Miriam; Koffijberg, Hendrik

2017-01-01

Background: Parametric distributions based on individual patient data can be used to represent both stochastic and parameter uncertainty. Although general guidance is available on how parameter uncertainty should be accounted for in probabilistic sensitivity analysis, there is no comprehensive

Accounting for parameter uncertainty in the definition of parametric distributions used to describe individual patient variation in health economic models

NARCIS (Netherlands)

Degeling, Koen; IJzerman, Maarten J; Koopman, Miriam; Koffijberg, Hendrik

2017-01-01

Background Parametric distributions based on individual patient data can be used to represent both stochastic and parameter uncertainty. Although general guidance is available on how parameter uncertainty should be accounted for in probabilistic sensitivity analysis, there is no comprehensive

Parameter discovery in stochastic biological models using simulated annealing and statistical model checking.

Science.gov (United States)

Hussain, Faraz; Jha, Sumit K; Jha, Susmit; Langmead, Christopher J

2014-01-01

Stochastic models are increasingly used to study the behaviour of biochemical systems. While the structure of such models is often readily available from first principles, unknown quantitative features of the model are incorporated into the model as parameters. Algorithmic discovery of parameter values from experimentally observed facts remains a challenge for the computational systems biology community. We present a new parameter discovery algorithm that uses simulated annealing, sequential hypothesis testing, and statistical model checking to learn the parameters in a stochastic model. We apply our technique to a model of glucose and insulin metabolism used for in-silico validation of artificial pancreata and demonstrate its effectiveness by developing parallel CUDA-based implementation for parameter synthesis in this model.

Models for describing the thermal characteristics of building components

DEFF Research Database (Denmark)

Jimenez, M.J.; Madsen, Henrik

2008-01-01

, for example. For the analysis of these tests, dynamic analysis models and methods are required. However, a wide variety of models and methods exists, and the problem of choosing the most appropriate approach for each particular case is a non-trivial and interdisciplinary task. Knowledge of a large family....... The characteristics of each type of model are highlighted. Some available software tools for each of the methods described will be mentioned. A case study also demonstrating the difference between linear and nonlinear models is considered....... of these approaches may therefore be very useful for selecting a suitable approach for each particular case. This paper presents an overview of models that can be applied for modelling the thermal characteristics of buildings and building components using data from outdoor testing. The choice of approach depends...

Fractional single-phase-lagging heat conduction model for describing anomalous diffusion

Directory of Open Access Journals (Sweden)

T.N. Mishra

2016-03-01

Full Text Available The fractional single-phase-lagging (FSPL heat conduction model is obtained by combining scalar time fractional conservation equation to the single-phase-lagging (SPL heat conduction model. Based on the FSPL heat conduction model, anomalous diffusion within a finite thin film is investigated. The effect of different parameters on solution has been observed and studied the asymptotic behavior of the FSPL model. The analytical solution is obtained using Laplace transform method. The whole analysis is presented in dimensionless form. Numerical examples of particular interest have been studied and discussed in details.

Establishing statistical models of manufacturing parameters

International Nuclear Information System (INIS)

Senevat, J.; Pape, J.L.; Deshayes, J.F.

1991-01-01

This paper reports on the effect of pilgering and cold-work parameters on contractile strain ratio and mechanical properties that were investigated using a large population of Zircaloy tubes. Statistical models were established between: contractile strain ratio and tooling parameters, mechanical properties (tensile test, creep test) and cold-work parameters, and mechanical properties and stress-relieving temperature

Lumped-parameter fuel rod model for rapid thermal transients

International Nuclear Information System (INIS)

Perkins, K.R.; Ramshaw, J.D.

1975-07-01

The thermal behavior of fuel rods during simulated accident conditions is extremely sensitive to the heat transfer coefficient which is, in turn, very sensitive to the cladding surface temperature and the fluid conditions. The development of a semianalytical, lumped-parameter fuel rod model which is intended to provide accurate calculations, in a minimum amount of computer time, of the thermal response of fuel rods during a simulated loss-of-coolant accident is described. The results show good agreement with calculations from a comprehensive fuel-rod code (FRAP-T) currently in use at Aerojet Nuclear Company

The Polyakov, Nambu and Jona-Lasinio model and its applications to describe the sub-nuclear particles

International Nuclear Information System (INIS)

Blanquier, E.

2013-01-01

To study the high energy nuclear physics and the associated phenomenon, as the quark-gluon plasma / hadronic matter phase transition, the Nambu and Jona-Lasinio model (NJL) appears as an interesting alternative to the Quantum Chromodynamics, not solvable at the considered energies. Indeed, the NJL model allows the description of quarks physics, at finite temperatures and densities. Furthermore, in order to try to correct a limitation of the NJL model, i.e. the absence of confinement, it was proposed a coupling of the quarks/antiquarks to a Polyakov loop, forming the PNJL model. The objective of this thesis is to see the possibilities offered by the NJL and PNJL models, to describe relevant sub-nuclear particles (quarks, mesons, diquarks and baryons), to study their interactions, and to proceed to a dynamical study involving these particles. After a recall of the useful tools, we modeled the u, d, s effective quarks and the mesons. Then, we described the baryons as quarks-diquarks bound states. A part of the work concerned the calculations of the cross-sections associated to the possible reactions implying these particles. Then, we incorporated these results in a computer code, in order to study the cooling of a quarks/antiquarks plasma and its hadronization. In this study, each particle evolves in a system in which the temperature and the densities are local parameters. We have two types of interactions: one due to the collisions, and the other is a remote interaction, notably between quarks. Finally, we studied the properties of our approach: qualities, limitations, and possible evolutions. (author)

Some tests for parameter constancy in cointegrated VAR-models

DEFF Research Database (Denmark)

Hansen, Henrik; Johansen, Søren

1999-01-01

Some methods for the evaluation of parameter constancy in vector autoregressive (VAR) models are discussed. Two different ways of re-estimating the VAR model are proposed; one in which all parameters are estimated recursively based upon the likelihood function for the first observations, and anot...... be applied to test the constancy of the long-run parameters in the cointegrated VAR-model. All results are illustrated using a model for the term structure of interest rates on US Treasury securities. ......Some methods for the evaluation of parameter constancy in vector autoregressive (VAR) models are discussed. Two different ways of re-estimating the VAR model are proposed; one in which all parameters are estimated recursively based upon the likelihood function for the first observations......, and another in which the cointegrating relations are estimated recursively from a likelihood function, where the short-run parameters have been concentrated out. We suggest graphical procedures based on recursively estimated eigenvalues to evaluate the constancy of the long-run parameters in the model...

Edge Modeling by Two Blur Parameters in Varying Contrasts.

Science.gov (United States)

Seo, Suyoung

2018-06-01

This paper presents a method of modeling edge profiles with two blur parameters, and estimating and predicting those edge parameters with varying brightness combinations and camera-to-object distances (COD). First, the validity of the edge model is proven mathematically. Then, it is proven experimentally with edges from a set of images captured for specifically designed target sheets and with edges from natural images. Estimation of the two blur parameters for each observed edge profile is performed with a brute-force method to find parameters that produce global minimum errors. Then, using the estimated blur parameters, actual blur parameters of edges with arbitrary brightness combinations are predicted using a surface interpolation method (i.e., kriging). The predicted surfaces show that the two blur parameters of the proposed edge model depend on both dark-side edge brightness and light-side edge brightness following a certain global trend. This is similar across varying CODs. The proposed edge model is compared with a one-blur parameter edge model using experiments of the root mean squared error for fitting the edge models to each observed edge profile. The comparison results suggest that the proposed edge model has superiority over the one-blur parameter edge model in most cases where edges have varying brightness combinations.

Approaches to highly parameterized inversion: A guide to using PEST for model-parameter and predictive-uncertainty analysis

Science.gov (United States)

Doherty, John E.; Hunt, Randall J.; Tonkin, Matthew J.

2010-01-01

Analysis of the uncertainty associated with parameters used by a numerical model, and with predictions that depend on those parameters, is fundamental to the use of modeling in support of decisionmaking. Unfortunately, predictive uncertainty analysis with regard to models can be very computationally demanding, due in part to complex constraints on parameters that arise from expert knowledge of system properties on the one hand (knowledge constraints) and from the necessity for the model parameters to assume values that allow the model to reproduce historical system behavior on the other hand (calibration constraints). Enforcement of knowledge and calibration constraints on parameters used by a model does not eliminate the uncertainty in those parameters. In fact, in many cases, enforcement of calibration constraints simply reduces the uncertainties associated with a number of broad-scale combinations of model parameters that collectively describe spatially averaged system properties. The uncertainties associated with other combinations of parameters, especially those that pertain to small-scale parameter heterogeneity, may not be reduced through the calibration process. To the extent that a prediction depends on system-property detail, its postcalibration variability may be reduced very little, if at all, by applying calibration constraints; knowledge constraints remain the only limits on the variability of predictions that depend on such detail. Regrettably, in many common modeling applications, these constraints are weak. Though the PEST software suite was initially developed as a tool for model calibration, recent developments have focused on the evaluation of model-parameter and predictive uncertainty. As a complement to functionality that it provides for highly parameterized inversion (calibration) by means of formal mathematical regularization techniques, the PEST suite provides utilities for linear and nonlinear error-variance and uncertainty analysis in

The fitting parameters extraction of conversion model of the low dose rate effect in bipolar devices

International Nuclear Information System (INIS)

Bakerenkov, Alexander

2011-01-01

The Enhanced Low Dose Rate Sensitivity (ELDRS) in bipolar devices consists of in base current degradation of NPN and PNP transistors increase as the dose rate is decreased. As a result of almost 20-year studying, the some physical models of effect are developed, being described in detail. Accelerated test methods, based on these models use in standards. The conversion model of the effect, that allows to describe the inverse S-shaped excess base current dependence versus dose rate, was proposed. This paper presents the problem of conversion model fitting parameters extraction.

Identifyability measures to select the parameters to be estimated in a solid-state fermentation distributed parameter model.

Science.gov (United States)

da Silveira, Christian L; Mazutti, Marcio A; Salau, Nina P G

2016-07-08

Process modeling can lead to of advantages such as helping in process control, reducing process costs and product quality improvement. This work proposes a solid-state fermentation distributed parameter model composed by seven differential equations with seventeen parameters to represent the process. Also, parameters estimation with a parameters identifyability analysis (PIA) is performed to build an accurate model with optimum parameters. Statistical tests were made to verify the model accuracy with the estimated parameters considering different assumptions. The results have shown that the model assuming substrate inhibition better represents the process. It was also shown that eight from the seventeen original model parameters were nonidentifiable and better results were obtained with the removal of these parameters from the estimation procedure. Therefore, PIA can be useful to estimation procedure, since it may reduce the number of parameters that can be evaluated. Further, PIA improved the model results, showing to be an important procedure to be taken. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:905-917, 2016. © 2016 American Institute of Chemical Engineers.

Sensitivity analysis of a radionuclide transfer model describing contaminated vegetation in Fukushima prefecture, using Morris and Sobol' - Application of sensitivity analysis on a radionuclides transfer model in the environment describing weeds contamination in Fukushima Prefecture, using Morris method and Sobol' indices indices

Energy Technology Data Exchange (ETDEWEB)

Nicoulaud-Gouin, V.; Metivier, J.M.; Gonze, M.A. [Institut de Radioprotection et de Surete Nucleaire-PRP-ENV/SERIS/LM2E (France); Garcia-Sanchez, L. [Institut de Radioprotection et de Surete Nucleaire-PRPENV/SERIS/L2BT (France)

2014-07-01

The increasing spatial and temporal complexity of models demands methods capable of ranking the influence of their large numbers of parameters. This question specifically arises in assessment studies on the consequences of the Fukushima accident. Sensitivity analysis aims at measuring the influence of input variability on the output response. Generally, two main approaches are distinguished (Saltelli, 2001, Iooss, 2011): - Screening approach, less expensive in computation time and allowing to identify non influential parameters; - Measures of importance, introducing finer quantitative indices. In this category, there are regression-based methods, assuming a linear or monotonic response (Pearson coefficient, Spearman coefficient), and variance-based methods, without assumptions on the model but requiring an increasingly prohibitive number of evaluations when the number of parameters increases. These approaches are available in various statistical programs (notably R) but are still poorly integrated in modelling platforms of radioecological risk assessment. This work aimed at illustrating the benefits of sensitivity analysis in the course of radioecological risk assessments This study used two complementary state-of-art global sensitivity analysis methods: - The screening method of Morris (Morris, 1991; Campolongo et al., 2007) based on limited model evaluations with a one-at-a-time (OAT) design; - The variance-based Sobol' sensitivity analysis (Saltelli, 2002) based a large number of model evaluations in the parameter space with a quasi-random sampling (Owen, 2003). Sensitivity analyses were applied on a dynamic Soil-Plant Deposition Model (Gonze et al., submitted to this conference) predicting foliar concentration in weeds after atmospheric radionuclide fallout. The Soil-Plant Deposition Model considers two foliage pools and a root pool, and describes foliar biomass growth with a Verhulst model. The developed semi-analytic formulation of foliar concentration

Model structural uncertainty quantification and hydrologic parameter and prediction error analysis using airborne electromagnetic data

DEFF Research Database (Denmark)

Minsley, B. J.; Christensen, Nikolaj Kruse; Christensen, Steen

Model structure, or the spatial arrangement of subsurface lithological units, is fundamental to the hydrological behavior of Earth systems. Knowledge of geological model structure is critically important in order to make informed hydrological predictions and management decisions. Model structure...... is never perfectly known, however, and incorrect assumptions can be a significant source of error when making model predictions. We describe a systematic approach for quantifying model structural uncertainty that is based on the integration of sparse borehole observations and large-scale airborne...... electromagnetic (AEM) data. Our estimates of model structural uncertainty follow a Bayesian framework that accounts for both the uncertainties in geophysical parameter estimates given AEM data, and the uncertainties in the relationship between lithology and geophysical parameters. Using geostatistical sequential...

Parameter estimation in stochastic rainfall-runoff models

DEFF Research Database (Denmark)

Jonsdottir, Harpa; Madsen, Henrik; Palsson, Olafur Petur

2006-01-01

A parameter estimation method for stochastic rainfall-runoff models is presented. The model considered in the paper is a conceptual stochastic model, formulated in continuous-discrete state space form. The model is small and a fully automatic optimization is, therefore, possible for estimating all...... the parameter values are optimal for simulation or prediction. The data originates from Iceland and the model is designed for Icelandic conditions, including a snow routine for mountainous areas. The model demands only two input data series, precipitation and temperature and one output data series...

What Population Reveals about Individual Cell Identity: Single-Cell Parameter Estimation of Models of Gene Expression in Yeast.

Directory of Open Access Journals (Sweden)

Artémis Llamosi

2016-02-01

Full Text Available Significant cell-to-cell heterogeneity is ubiquitously observed in isogenic cell populations. Consequently, parameters of models of intracellular processes, usually fitted to population-averaged data, should rather be fitted to individual cells to obtain a population of models of similar but non-identical individuals. Here, we propose a quantitative modeling framework that attributes specific parameter values to single cells for a standard model of gene expression. We combine high quality single-cell measurements of the response of yeast cells to repeated hyperosmotic shocks and state-of-the-art statistical inference approaches for mixed-effects models to infer multidimensional parameter distributions describing the population, and then derive specific parameters for individual cells. The analysis of single-cell parameters shows that single-cell identity (e.g. gene expression dynamics, cell size, growth rate, mother-daughter relationships is, at least partially, captured by the parameter values of gene expression models (e.g. rates of transcription, translation and degradation. Our approach shows how to use the rich information contained into longitudinal single-cell data to infer parameters that can faithfully represent single-cell identity.

How Mathematics Describes Life

Science.gov (United States)

Teklu, Abraham

2017-01-01

The circle of life is something we have all heard of from somewhere, but we don't usually try to calculate it. For some time we have been working on analyzing a predator-prey model to better understand how mathematics can describe life, in particular the interaction between two different species. The model we are analyzing is called the Holling-Tanner model, and it cannot be solved analytically. The Holling-Tanner model is a very common model in population dynamics because it is a simple descriptor of how predators and prey interact. The model is a system of two differential equations. The model is not specific to any particular set of species and so it can describe predator-prey species ranging from lions and zebras to white blood cells and infections. One thing all these systems have in common are critical points. A critical point is a value for both populations that keeps both populations constant. It is important because at this point the differential equations are equal to zero. For this model there are two critical points, a predator free critical point and a coexistence critical point. Most of the analysis we did is on the coexistence critical point because the predator free critical point is always unstable and frankly less interesting than the coexistence critical point. What we did is consider two regimes for the differential equations, large B and small B. B, A, and C are parameters in the differential equations that control the system where B measures how responsive the predators are to change in the population, A represents predation of the prey, and C represents the satiation point of the prey population. For the large B case we were able to approximate the system of differential equations by a single scalar equation. For the small B case we were able to predict the limit cycle. The limit cycle is a process of the predator and prey populations growing and shrinking periodically. This model has a limit cycle in the regime of small B, that we solved for

Application of lumped-parameter models

DEFF Research Database (Denmark)

Ibsen, Lars Bo; Liingaard, Morten

This technical report concerns the lumped-parameter models for a suction caisson with a ratio between skirt length and foundation diameter equal to 1/2, embedded into an viscoelastic soil. The models are presented for three different values of the shear modulus of the subsoil (section 1.1). Subse...

A simulation of water pollution model parameter estimation

Science.gov (United States)

Kibler, J. F.

1976-01-01

A parameter estimation procedure for a water pollution transport model is elaborated. A two-dimensional instantaneous-release shear-diffusion model serves as representative of a simple transport process. Pollution concentration levels are arrived at via modeling of a remote-sensing system. The remote-sensed data are simulated by adding Gaussian noise to the concentration level values generated via the transport model. Model parameters are estimated from the simulated data using a least-squares batch processor. Resolution, sensor array size, and number and location of sensor readings can be found from the accuracies of the parameter estimates.

Identification of ecosystem parameters by SDE-modelling

DEFF Research Database (Denmark)

Stochastic differential equations (SDEs) for ecosystem modelling have attracted increasing attention during recent years. The modelling has mostly been through simulation experiments in order to analyse how system noise propagates through the ordinary differential equation formulation of ecosystem...... models. Estimation of parameters in SDEs is, however, possible by combining Kalman filter techniques and likelihood estimation. By modelling parameters as random walks it is possible to identify linear as well as non-linear interactions between ecosystem components. By formulating a simple linear SDE...

Han's model parameters for microalgae grown under intermittent illumination: Determined using particle swarm optimization.

Science.gov (United States)

Pozzobon, Victor; Perre, Patrick

2018-01-21

This work provides a model and the associated set of parameters allowing for microalgae population growth computation under intermittent lightning. Han's model is coupled with a simple microalgae growth model to yield a relationship between illumination and population growth. The model parameters were obtained by fitting a dataset available in literature using Particle Swarm Optimization method. In their work, authors grew microalgae in excess of nutrients under flashing conditions. Light/dark cycles used for these experimentations are quite close to those found in photobioreactor, i.e. ranging from several seconds to one minute. In this work, in addition to producing the set of parameters, Particle Swarm Optimization robustness was assessed. To do so, two different swarm initialization techniques were used, i.e. uniform and random distribution throughout the search-space. Both yielded the same results. In addition, swarm distribution analysis reveals that the swarm converges to a unique minimum. Thus, the produced set of parameters can be trustfully used to link light intensity to population growth rate. Furthermore, the set is capable to describe photodamages effects on population growth. Hence, accounting for light overexposure effect on algal growth. Copyright © 2017 Elsevier Ltd. All rights reserved.

Interpretation of hydraulic conductivity data and parameter evaluation for groundwater flow models

International Nuclear Information System (INIS)

Niemi, A.

1991-01-01

The report reviews recent developments in evaluating effective permeabilities for groundwater flow models, starting from methods of well test interpretation for and proceeding to the principles of parameter estimation. Basic concepts of parameter evaluation as well as expressions derived for effective permeabilities in traditional porous medium are described. Due to the assumptions made, these do often not apply for fractured media. Specific features of fractured medium are discussed, including approaches used determining the size of a possible REV and questions related to the application of stochastic theories. Due to the difficulties encountered when applying traditional deterministic models for fractured media, stochastic and fracture network approaches have been developed. The application of these techniques is still under development, the main questions to be resolved being related to the scarcity of data

Spatio-temporal modeling of nonlinear distributed parameter systems

CERN Document Server

Li, Han-Xiong

2011-01-01

The purpose of this volume is to provide a brief review of the previous work on model reduction and identifi cation of distributed parameter systems (DPS), and develop new spatio-temporal models and their relevant identifi cation approaches. In this book, a systematic overview and classifi cation on the modeling of DPS is presented fi rst, which includes model reduction, parameter estimation and system identifi cation. Next, a class of block-oriented nonlinear systems in traditional lumped parameter systems (LPS) is extended to DPS, which results in the spatio-temporal Wiener and Hammerstein s

Models for estimating photosynthesis parameters from in situ production profiles

Science.gov (United States)

Kovač, Žarko; Platt, Trevor; Sathyendranath, Shubha; Antunović, Suzana

2017-12-01

The rate of carbon assimilation in phytoplankton primary production models is mathematically prescribed with photosynthesis irradiance functions, which convert a light flux (energy) into a material flux (carbon). Information on this rate is contained in photosynthesis parameters: the initial slope and the assimilation number. The exactness of parameter values is crucial for precise calculation of primary production. Here we use a model of the daily production profile based on a suite of photosynthesis irradiance functions and extract photosynthesis parameters from in situ measured daily production profiles at the Hawaii Ocean Time-series station Aloha. For each function we recover parameter values, establish parameter distributions and quantify model skill. We observe that the choice of the photosynthesis irradiance function to estimate the photosynthesis parameters affects the magnitudes of parameter values as recovered from in situ profiles. We also tackle the problem of parameter exchange amongst the models and the effect it has on model performance. All models displayed little or no bias prior to parameter exchange, but significant bias following parameter exchange. The best model performance resulted from using optimal parameter values. Model formulation was extended further by accounting for spectral effects and deriving a spectral analytical solution for the daily production profile. The daily production profile was also formulated with time dependent growing biomass governed by a growth equation. The work on parameter recovery was further extended by exploring how to extract photosynthesis parameters from information on watercolumn production. It was demonstrated how to estimate parameter values based on a linearization of the full analytical solution for normalized watercolumn production and from the solution itself, without linearization. The paper complements previous works on photosynthesis irradiance models by analysing the skill and consistency of

A new settling velocity model to describe secondary sedimentation.

Science.gov (United States)

Ramin, Elham; Wágner, Dorottya S; Yde, Lars; Binning, Philip J; Rasmussen, Michael R; Mikkelsen, Peter Steen; Plósz, Benedek Gy

2014-12-01

Secondary settling tanks (SSTs) are the most hydraulically sensitive unit operations in biological wastewater treatment plants. The maximum permissible inflow to the plant depends on the efficiency of SSTs in separating and thickening the activated sludge. The flow conditions and solids distribution in SSTs can be predicted using computational fluid dynamics (CFD) tools. Despite extensive studies on the compression settling behaviour of activated sludge and the development of advanced settling velocity models for use in SST simulations, these models are not often used, due to the challenges associated with their calibration. In this study, we developed a new settling velocity model, including hindered, transient and compression settling, and showed that it can be calibrated to data from a simple, novel settling column experimental set-up using the Bayesian optimization method DREAM(ZS). In addition, correlations between the Herschel-Bulkley rheological model parameters and sludge concentration were identified with data from batch rheological experiments. A 2-D axisymmetric CFD model of a circular SST containing the new settling velocity and rheological model was validated with full-scale measurements. Finally, it was shown that the representation of compression settling in the CFD model can significantly influence the prediction of sludge distribution in the SSTs under dry- and wet-weather flow conditions. Copyright © 2014 Elsevier Ltd. All rights reserved.

An analytical mechanical model to describe the response of NiTi rotary endodontic files in a curved root canal.

Science.gov (United States)

Leroy, Agnès Marie Françoise; Bahia, Maria Guiomar de Azevedo; Ehrlacher, Alain; Buono, Vicente Tadeu Lopes

2012-08-01

To build a mathematical model describing the mechanical behavior of NiTi rotary files while they are rotating in a root canal. The file was seen as a beam undergoing large transformations. The instrument was assumed to be rotating steadily in the root canal, and the geometry of the canal was considered as a known parameter of the problem. The formulae of large transformations mechanics then allowed the calculation of the Green-Lagrange strain field in the file. The non-linear mechanical behavior of NiTi was modeled as a continuous piecewise linear function, assuming that the material did not reach plastic deformation. Criteria locating the changes of behavior of NiTi were established and the tension field in the file, and the external efforts applied on it were calculated. The unknown variable of torsion was deduced from the equilibrium equation system using a Coulomb contact law which solved the problem on a cycle of rotation. In order to verify that the model described well reality, three-point bending experiments were managed on superelastic NiTi wires, whose results were compared to the theoretical ones. It appeared that the model gave a good mentoring of the empirical results in the range of bending angles that interested us. Knowing the geometry of the root canal, one is now able to write the equations of the strain and stress fields in the endodontic instrument, and to quantify the impact of each macroscopic parameter of the problem on its response. This should be useful to predict failure of the files under rotating bending fatigue, and to optimize the geometry of the files. Copyright © 2012 Elsevier B.V. All rights reserved.

GEMSFITS: Code package for optimization of geochemical model parameters and inverse modeling

International Nuclear Information System (INIS)

Miron, George D.; Kulik, Dmitrii A.; Dmytrieva, Svitlana V.; Wagner, Thomas

2015-01-01

Highlights: • Tool for generating consistent parameters against various types of experiments. • Handles a large number of experimental data and parameters (is parallelized). • Has a graphical interface and can perform statistical analysis on the parameters. • Tested on fitting the standard state Gibbs free energies of aqueous Al species. • Example on fitting interaction parameters of mixing models and thermobarometry. - Abstract: GEMSFITS is a new code package for fitting internally consistent input parameters of GEM (Gibbs Energy Minimization) geochemical–thermodynamic models against various types of experimental or geochemical data, and for performing inverse modeling tasks. It consists of the gemsfit2 (parameter optimizer) and gfshell2 (graphical user interface) programs both accessing a NoSQL database, all developed with flexibility, generality, efficiency, and user friendliness in mind. The parameter optimizer gemsfit2 includes the GEMS3K chemical speciation solver ( (http://gems.web.psi.ch/GEMS3K)), which features a comprehensive suite of non-ideal activity- and equation-of-state models of solution phases (aqueous electrolyte, gas and fluid mixtures, solid solutions, (ad)sorption. The gemsfit2 code uses the robust open-source NLopt library for parameter fitting, which provides a selection between several nonlinear optimization algorithms (global, local, gradient-based), and supports large-scale parallelization. The gemsfit2 code can also perform comprehensive statistical analysis of the fitted parameters (basic statistics, sensitivity, Monte Carlo confidence intervals), thus supporting the user with powerful tools for evaluating the quality of the fits and the physical significance of the model parameters. The gfshell2 code provides menu-driven setup of optimization options (data selection, properties to fit and their constraints, measured properties to compare with computed counterparts, and statistics). The practical utility, efficiency, and

Identification of parameters of discrete-continuous models

International Nuclear Information System (INIS)

Cekus, Dawid; Warys, Pawel

2015-01-01

In the paper, the parameters of a discrete-continuous model have been identified on the basis of experimental investigations and formulation of optimization problem. The discrete-continuous model represents a cantilever stepped Timoshenko beam. The mathematical model has been formulated and solved according to the Lagrange multiplier formalism. Optimization has been based on the genetic algorithm. The presented proceeding’s stages make the identification of any parameters of discrete-continuous systems possible

Identification of parameters of discrete-continuous models

Energy Technology Data Exchange (ETDEWEB)

Cekus, Dawid, E-mail: cekus@imipkm.pcz.pl; Warys, Pawel, E-mail: warys@imipkm.pcz.pl [Institute of Mechanics and Machine Design Foundations, Czestochowa University of Technology, Dabrowskiego 73, 42-201 Czestochowa (Poland)

2015-03-10

In the paper, the parameters of a discrete-continuous model have been identified on the basis of experimental investigations and formulation of optimization problem. The discrete-continuous model represents a cantilever stepped Timoshenko beam. The mathematical model has been formulated and solved according to the Lagrange multiplier formalism. Optimization has been based on the genetic algorithm. The presented proceeding’s stages make the identification of any parameters of discrete-continuous systems possible.

Identifying the connective strength between model parameters and performance criteria

Directory of Open Access Journals (Sweden)

B. Guse

2017-11-01

Full Text Available In hydrological models, parameters are used to represent the time-invariant characteristics of catchments and to capture different aspects of hydrological response. Hence, model parameters need to be identified based on their role in controlling the hydrological behaviour. For the identification of meaningful parameter values, multiple and complementary performance criteria are used that compare modelled and measured discharge time series. The reliability of the identification of hydrologically meaningful model parameter values depends on how distinctly a model parameter can be assigned to one of the performance criteria. To investigate this, we introduce the new concept of connective strength between model parameters and performance criteria. The connective strength assesses the intensity in the interrelationship between model parameters and performance criteria in a bijective way. In our analysis of connective strength, model simulations are carried out based on a latin hypercube sampling. Ten performance criteria including Nash–Sutcliffe efficiency (NSE, Kling–Gupta efficiency (KGE and its three components (alpha, beta and r as well as RSR (the ratio of the root mean square error to the standard deviation for different segments of the flow duration curve (FDC are calculated. With a joint analysis of two regression tree (RT approaches, we derive how a model parameter is connected to different performance criteria. At first, RTs are constructed using each performance criterion as the target variable to detect the most relevant model parameters for each performance criterion. Secondly, RTs are constructed using each parameter as the target variable to detect which performance criteria are impacted by changes in the values of one distinct model parameter. Based on this, appropriate performance criteria are identified for each model parameter. In this study, a high bijective connective strength between model parameters and performance criteria

Full-range stress–strain behaviour of contemporary pipeline steels: Part II. Estimation of model parameters

International Nuclear Information System (INIS)

Hertelé, Stijn; De Waele, Wim; Denys, Rudi; Verstraete, Matthias

2012-01-01

Contemporary pipeline steels with a yield-to-tensile ratio above 0.80 often show two-stages of strain hardening, which cannot be simultaneously described by the standardized Ramberg–Osgood model. A companion paper (Part I) showed that the recently developed UGent model provides more accurate descriptions than the Ramberg–Osgood model, as it succeeds in describing both strain hardening stages. However, it may be challenging to obtain an optimal model fit in absence of full stress–strain data. This paper discusses on how to find suited parameter values for the UGent model, given a set of measurable tensile test characteristics. The proposed methodology shows good results for an extensive set of investigated experimental stress–strain curves. Next to some common tensile test characteristics, the 1.0% proof stress is needed. The authors therefore encourage the acquisition of this stress during tensile tests. - Highlights: ► An analytical procedure estimates UGent model parameters. ► The procedure requires a set of tensile test characteristics. ► The UGent model performs better than the Ramberg–Osgood model. ► Apart from common characteristics, the 1.0% proof stress is required. ► The authors encourage the acquisition of this 1.0% proof stress.

Brownian motion model with stochastic parameters for asset prices

Science.gov (United States)

Ching, Soo Huei; Hin, Pooi Ah

2013-09-01

The Brownian motion model may not be a completely realistic model for asset prices because in real asset prices the drift μ and volatility σ may change over time. Presently we consider a model in which the parameter x = (μ,σ) is such that its value x (t + Δt) at a short time Δt ahead of the present time t depends on the value of the asset price at time t + Δt as well as the present parameter value x(t) and m-1 other parameter values before time t via a conditional distribution. The Malaysian stock prices are used to compare the performance of the Brownian motion model with fixed parameter with that of the model with stochastic parameter.

Study on Parameters Modeling of Wind Turbines Using SCADA Data

Directory of Open Access Journals (Sweden)

Yonglong YAN

2014-08-01

Full Text Available Taking the advantage of the current massive monitoring data from Supervisory Control and Data Acquisition (SCADA system of wind farm, it is of important significance for anomaly detection, early warning and fault diagnosis to build the data model of state parameters of wind turbines (WTs. The operational conditions and the relationships between the state parameters of wind turbines are complex. It is difficult to establish the model of state parameter accurately, and the modeling method of state parameters of wind turbines considering parameter selection is proposed. Firstly, by analyzing the characteristic of SCADA data, a reasonable range of data and monitoring parameters are chosen. Secondly, neural network algorithm is adapted, and the selection method of input parameters in the model is presented. Generator bearing temperature and cooling air temperature are regarded as target parameters, and the two models are built and input parameters of the models are selected, respectively. Finally, the parameter selection method in this paper and the method using genetic algorithm-partial least square (GA-PLS are analyzed comparatively, and the results show that the proposed methods are correct and effective. Furthermore, the modeling of two parameters illustrate that the method in this paper can applied to other state parameters of wind turbines.

Computational Models Describing Possible Mechanisms for Generation of Excessive Beta Oscillations in Parkinson's Disease.

Directory of Open Access Journals (Sweden)

Alex Pavlides

2015-12-01

Full Text Available In Parkinson's disease, an increase in beta oscillations within the basal ganglia nuclei has been shown to be associated with difficulty in movement initiation. An important role in the generation of these oscillations is thought to be played by the motor cortex and by a network composed of the subthalamic nucleus (STN and the external segment of globus pallidus (GPe. Several alternative models have been proposed to describe the mechanisms for generation of the Parkinsonian beta oscillations. However, a recent experimental study of Tachibana and colleagues yielded results which are challenging for all published computational models of beta generation. That study investigated how the presence of beta oscillations in a primate model of Parkinson's disease is affected by blocking different connections of the STN-GPe circuit. Due to a large number of experimental conditions, the study provides strong constraints that any mechanistic model of beta generation should satisfy. In this paper we present two models consistent with the data of Tachibana et al. The first model assumes that Parkinsonian beta oscillation are generated in the cortex and the STN-GPe circuits resonates at this frequency. The second model additionally assumes that the feedback from STN-GPe circuit to cortex is important for maintaining the oscillations in the network. Predictions are made about experimental evidence that is required to differentiate between the two models, both of which are able to reproduce firing rates, oscillation frequency and effects of lesions carried out by Tachibana and colleagues. Furthermore, an analysis of the models reveals how the amplitude and frequency of the generated oscillations depend on parameters.

Genetic Algorithms for Estimating Effective Parameters in a Lumped Reactor Model for Reactivity Predictions

International Nuclear Information System (INIS)

Marseguerra, Marzio; Zio, Enrico

2001-01-01

The control system of a reactor should be able to predict, in real time, the amount of reactivity to be inserted (e.g., by control rod movements and boron injection and dilution) to respond to a given electrical load demand or to undesired, accidental transients. The real-time constraint renders impractical the use of a large, detailed dynamic reactor code. One has, then, to resort to simplified analytical models with lumped effective parameters suitably estimated from the reactor data.The simple and well-known Chernick model for describing the reactor power evolution in the presence of xenon is considered and the feasibility of using genetic algorithms for estimating the effective nuclear parameters involved and the initial nonmeasurable xenon and iodine conditions is investigated. This approach has the advantage of counterbalancing the inherent model simplicity with the periodic reestimation of the effective parameter values pertaining to each reactor on the basis of its recent history. By so doing, other effects, such as burnup, are automatically taken into account

Lumped-parameter Model of a Bucket Foundation

DEFF Research Database (Denmark)

Andersen, Lars; Ibsen, Lars Bo; Liingaard, Morten

2009-01-01

efficient model that can be applied in aero-elastic codes for fast evaluation of the dynamic structural response of wind turbines. The target solutions, utilised for calibration of the lumped-parameter models, are obtained by a coupled finite-element/boundaryelement scheme in the frequency domain......, and the quality of the models are tested in the time and frequency domains. It is found that precise results are achieved by lumped-parameter models with two to four internal degrees of freedom per displacement or rotation of the foundation. Further, coupling between the horizontal sliding and rocking cannot...

Source term modelling parameters for Project-90

International Nuclear Information System (INIS)

Shaw, W.; Smith, G.; Worgan, K.; Hodgkinson, D.; Andersson, K.

1992-04-01

This document summarises the input parameters for the source term modelling within Project-90. In the first place, the parameters relate to the CALIBRE near-field code which was developed for the Swedish Nuclear Power Inspectorate's (SKI) Project-90 reference repository safety assessment exercise. An attempt has been made to give best estimate values and, where appropriate, a range which is related to variations around base cases. It should be noted that the data sets contain amendments to those considered by KBS-3. In particular, a completely new set of inventory data has been incorporated. The information given here does not constitute a complete set of parameter values for all parts of the CALIBRE code. Rather, it gives the key parameter values which are used in the constituent models within CALIBRE and the associated studies. For example, the inventory data acts as an input to the calculation of the oxidant production rates, which influence the generation of a redox front. The same data is also an initial value data set for the radionuclide migration component of CALIBRE. Similarly, the geometrical parameters of the near-field are common to both sub-models. The principal common parameters are gathered here for ease of reference and avoidance of unnecessary duplication and transcription errors. (au)

Incorporation of FcRn-mediated disposition model to describe the population pharmacokinetics of therapeutic monoclonal IgG antibody in clinical patients.

Science.gov (United States)

Ng, Chee M

2016-03-01

The two-compartment linear model used to describe the population pharmacokinetics (PK) of many therapeutic monoclonal antibodies (TMAbs) offered little biological insight to antibody disposition in humans. The purpose of this study is to develop a semi-mechanistic FcRn-mediated IgG disposition model to describe the population PK of TMAbs in clinical patients. A standard two-compartment linear PK model from a previously published population PK model of pertuzumab was used to simulate intensive PK data of 100 subjects for model development. Two different semi-mechanistic FcRn-mediated IgG disposition models were developed and First Order Conditional Estimation (FOCE) with the interaction method in NONMEM was used to obtain the final model estimates. The performances of these models were then compared with the two-compartment linear PK model used to simulate the data for model development. A semi-mechanistic FcRn-mediated IgG disposition model consisting of a peripheral tissue compartment and FcRn-containing endosomes in the central compartment best describes the simulated pertuzumab population PK data. This developed semi-mechanistic population PK model had the same number of model parameters, produced very similar concentration-time profiles but provided additional biological insight to the FcRn-mediated IgG disposition in human subjects compared with the standard linear two-compartment linear PK model. This first reported semi-mechanistic model may serve as an important model framework for developing future population PK models of TMAbs in clinical patients. Copyright © 2015 John Wiley & Sons, Ltd.

A new model describing the curves for repair of both DNA double-strand breaks and chromosome damage

International Nuclear Information System (INIS)

Foray, N.; Badie, C.; Alsbeih, G.; Malaise, E.P.; Fertil, B.

1996-01-01

A review of reports dealing with fittings of the data for repair of DNA double-strand breaks (DSBs) and excess chromosome fragments (ECFs) shows that several models are used to fit the repair curves. Since DSBs and ECFs are correleated, it is worth developing a model describing both phenomena. The curve-fitting models used most extensively, the two repair half-times model for DSBs and the monoexponential plus residual model for ECFs, appear to be too inflexible to describe the repair curves for both DSBs and ECFs. We have therefore developed a new concept based on a variable repair half-time. According to this concept, the repair curve is continuously bending and dependent on time and probably reflects a continuous spectrum of damage repairability. The fits of the curves for DSB repair to the variable repair half-time and the variable repair half-time plus residual models were compared to those obtained with the two half-times plus residual and two half-times models. Similarly, the fits of the curves for ECF repair to the variable repair half-time and variable half-time plus residual models were compared to that obtained with the monoexponential plus residual model. The quality of fit and the dependence of adjustable parameters on the portion of the curve fitted were used as comparison criteria. We found that: (a) It is useful to postulate the existence of a residual term for unrepairable lesions, regardless of the model adopted. (b) With the two cell lines tested (a normal and a hypersensitive one), data for both DSBs and ECTs are best fitted to the variable repair half-time plus residual model, whatever the repair time range. 47 refs., 3 figs., 3 tabs

An analytical mechanical model to describe the response of NiTi rotary endodontic files in a curved root canal

Energy Technology Data Exchange (ETDEWEB)

Leroy, Agnes Marie Francoise [Department of Metallurgical and Materials Engineering, School of Engineering, Universidade Federal de Minas Gerais (UFMG), Belo Horizonte, MG (Brazil); Department of Mechanical and Materials Engineering, Ecole des Ponts Paristech (ENPC), Champs-sur-Marne (France); Bahia, Maria Guiomar de Azevedo [Department of Restorative Dentistry, Faculty of Dentistry, Universidade Federal de Minas Gerais (UFMG), Belo Horizonte, MG (Brazil); Ehrlacher, Alain [Department of Mechanical and Materials Engineering, Ecole des Ponts Paristech (ENPC), Champs-sur-Marne (France); Buono, Vicente Tadeu Lopes, E-mail: vbuono@demet.ufmg.br [Department of Metallurgical and Materials Engineering, School of Engineering, Universidade Federal de Minas Gerais (UFMG), Belo Horizonte, MG (Brazil)

2012-08-01

Aim: To build a mathematical model describing the mechanical behavior of NiTi rotary files while they are rotating in a root canal. Methodology: The file was seen as a beam undergoing large transformations. The instrument was assumed to be rotating steadily in the root canal, and the geometry of the canal was considered as a known parameter of the problem. The formulae of large transformations mechanics then allowed the calculation of the Green-Lagrange strain field in the file. The non-linear mechanical behavior of NiTi was modeled as a continuous piecewise linear function, assuming that the material did not reach plastic deformation. Criteria locating the changes of behavior of NiTi were established and the tension field in the file, and the external efforts applied on it were calculated. The unknown variable of torsion was deduced from the equilibrium equation system using a Coulomb contact law which solved the problem on a cycle of rotation. Results: In order to verify that the model described well reality, three-point bending experiments were managed on superelastic NiTi wires, whose results were compared to the theoretical ones. It appeared that the model gave a good mentoring of the empirical results in the range of bending angles that interested us. Conclusions: Knowing the geometry of the root canal, one is now able to write the equations of the strain and stress fields in the endodontic instrument, and to quantify the impact of each macroscopic parameter of the problem on its response. This should be useful to predict failure of the files under rotating bending fatigue, and to optimize the geometry of the files. - Highlights: Black-Right-Pointing-Pointer A mechanical model of the behavior of a NiTi endodontic instrument was developed. Black-Right-Pointing-Pointer The model was validated with results of three-point bending tests on NiTi wires. Black-Right-Pointing-Pointer The model is appropriate for the optimization of instruments' geometry.

An analytical mechanical model to describe the response of NiTi rotary endodontic files in a curved root canal

International Nuclear Information System (INIS)

Leroy, Agnès Marie Françoise; Bahia, Maria Guiomar de Azevedo; Ehrlacher, Alain; Buono, Vicente Tadeu Lopes

2012-01-01

Aim: To build a mathematical model describing the mechanical behavior of NiTi rotary files while they are rotating in a root canal. Methodology: The file was seen as a beam undergoing large transformations. The instrument was assumed to be rotating steadily in the root canal, and the geometry of the canal was considered as a known parameter of the problem. The formulae of large transformations mechanics then allowed the calculation of the Green–Lagrange strain field in the file. The non-linear mechanical behavior of NiTi was modeled as a continuous piecewise linear function, assuming that the material did not reach plastic deformation. Criteria locating the changes of behavior of NiTi were established and the tension field in the file, and the external efforts applied on it were calculated. The unknown variable of torsion was deduced from the equilibrium equation system using a Coulomb contact law which solved the problem on a cycle of rotation. Results: In order to verify that the model described well reality, three-point bending experiments were managed on superelastic NiTi wires, whose results were compared to the theoretical ones. It appeared that the model gave a good mentoring of the empirical results in the range of bending angles that interested us. Conclusions: Knowing the geometry of the root canal, one is now able to write the equations of the strain and stress fields in the endodontic instrument, and to quantify the impact of each macroscopic parameter of the problem on its response. This should be useful to predict failure of the files under rotating bending fatigue, and to optimize the geometry of the files. - Highlights: ► A mechanical model of the behavior of a NiTi endodontic instrument was developed. ► The model was validated with results of three-point bending tests on NiTi wires. ► The model is appropriate for the optimization of instruments' geometry.

Improving fault image by determination of optimum seismic survey parameters using ray-based modeling

Science.gov (United States)

Saffarzadeh, Sadegh; Javaherian, Abdolrahim; Hasani, Hossein; Talebi, Mohammad Ali

2018-06-01

In complex structures such as faults, salt domes and reefs, specifying the survey parameters is more challenging and critical owing to the complicated wave field behavior involved in such structures. In the petroleum industry, detecting faults has become crucial for reservoir potential where faults can act as traps for hydrocarbon. In this regard, seismic survey modeling is employed to construct a model close to the real structure, and obtain very realistic synthetic seismic data. Seismic modeling software, the velocity model and parameters pre-determined by conventional methods enable a seismic survey designer to run a shot-by-shot virtual survey operation. A reliable velocity model of structures can be constructed by integrating the 2D seismic data, geological reports and the well information. The effects of various survey designs can be investigated by the analysis of illumination maps and flower plots. Also, seismic processing of the synthetic data output can describe the target image using different survey parameters. Therefore, seismic modeling is one of the most economical ways to establish and test the optimum acquisition parameters to obtain the best image when dealing with complex geological structures. The primary objective of this study is to design a proper 3D seismic survey orientation to achieve fault zone structures through ray-tracing seismic modeling. The results prove that a seismic survey designer can enhance the image of fault planes in a seismic section by utilizing the proposed modeling and processing approach.

Steady And Unsteady Lumped-Parameter Models For Determination of Groundwater Residence Time Distribution

International Nuclear Information System (INIS)

Oezyurt, N.N.

2002-01-01

Groundwater's residence time distribution is an important information to identify the transport mechanism in aquifer systems. In the absence or scarcity of geometric, hydraulic and geohydrologic data needed to describe a flow system, lumped parameter models, that handle the flow system as a whole, exist as an alternative to determine the residence time distribution. Lumped parametre models comprise of idealized models of piston and well-mixed flow and their combinations. Aquifer properties such as, dead volume and by-pass flow can also be included in these models. With the aid of these models, conceptual aquifer models can be tested and residence time distribution of groundwater can be determined

Global sensitivity analysis of a model related to memory formation in synapses: Model reduction based on epistemic parameter uncertainties and related issues.

Science.gov (United States)

Kulasiri, Don; Liang, Jingyi; He, Yao; Samarasinghe, Sandhya

2017-04-21

We investigate the epistemic uncertainties of parameters of a mathematical model that describes the dynamics of CaMKII-NMDAR complex related to memory formation in synapses using global sensitivity analysis (GSA). The model, which was published in this journal, is nonlinear and complex with Ca 2+ patterns with different level of frequencies as inputs. We explore the effects of parameter on the key outputs of the model to discover the most sensitive ones using GSA and partial ranking correlation coefficient (PRCC) and to understand why they are sensitive and others are not based on the biology of the problem. We also extend the model to add presynaptic neurotransmitter vesicles release to have action potentials as inputs of different frequencies. We perform GSA on this extended model to show that the parameter sensitivities are different for the extended model as shown by PRCC landscapes. Based on the results of GSA and PRCC, we reduce the original model to a less complex model taking the most important biological processes into account. We validate the reduced model against the outputs of the original model. We show that the parameter sensitivities are dependent on the inputs and GSA would make us understand the sensitivities and the importance of the parameters. A thorough phenomenological understanding of the relationships involved is essential to interpret the results of GSA and hence for the possible model reduction. Copyright © 2017 Elsevier Ltd. All rights reserved.

The mobilisation model and parameter sensitivity

International Nuclear Information System (INIS)

Blok, B.M.

1993-12-01

In the PRObabillistic Safety Assessment (PROSA) of radioactive waste in a salt repository one of the nuclide release scenario's is the subrosion scenario. A new subrosion model SUBRECN has been developed. In this model the combined effect of a depth-dependent subrosion, glass dissolution, and salt rise has been taken into account. The subrosion model SUBRECN and the implementation of this model in the German computer program EMOS4 is presented. A new computer program PANTER is derived from EMOS4. PANTER models releases of radionuclides via subrosion from a disposal site in a salt pillar into the biosphere. For uncertainty and sensitivity analyses the new subrosion model Latin Hypercube Sampling has been used for determine the different values for the uncertain parameters. The influence of the uncertainty in the parameters on the dose calculations has been investigated by the following sensitivity techniques: Spearman Rank Correlation Coefficients, Partial Rank Correlation Coefficients, Standardised Rank Regression Coefficients, and the Smirnov Test. (orig./HP)

Parameter estimations in predictive microbiology: Statistically sound modelling of the microbial growth rate.

Science.gov (United States)

Akkermans, Simen; Logist, Filip; Van Impe, Jan F

2018-04-01

When building models to describe the effect of environmental conditions on the microbial growth rate, parameter estimations can be performed either with a one-step method, i.e., directly on the cell density measurements, or in a two-step method, i.e., via the estimated growth rates. The two-step method is often preferred due to its simplicity. The current research demonstrates that the two-step method is, however, only valid if the correct data transformation is applied and a strict experimental protocol is followed for all experiments. Based on a simulation study and a mathematical derivation, it was demonstrated that the logarithm of the growth rate should be used as a variance stabilizing transformation. Moreover, the one-step method leads to a more accurate estimation of the model parameters and a better approximation of the confidence intervals on the estimated parameters. Therefore, the one-step method is preferred and the two-step method should be avoided. Copyright © 2017. Published by Elsevier Ltd.

Models and parameters for environmental radiological assessments

International Nuclear Information System (INIS)

Miller, C.W.

1983-01-01

This article reviews the forthcoming book Models and Parameters for Environmental Radiological Assessments, which presents a unified compilation of models and parameters for assessing the impact on man of radioactive discharges, both routine and accidental, into the environment. Models presented in this book include those developed for the prediction of atmospheric and hydrologic transport and deposition, for terrestrial and aquatic food-chain bioaccumulation, and for internal and external dosimetry. Summaries are presented for each of the transport and dosimetry areas previously for each of the transport and dosimetry areas previously mentioned, and details are available in the literature cited. A chapter of example problems illustrates many of the methodologies presented throughout the text. Models and parameters presented are based on the results of extensive literature reviews and evaluations performed primarily by the staff of the Health and Safety Research Division of Oak Ridge National Laboratory

Icosahedral symmetry described by an incommensurately modulated crystal structure model

DEFF Research Database (Denmark)

Wolny, Janusz; Lebech, Bente

1986-01-01

A crystal structure model of an incommensurately modulated structure is presented. Although six different reciprocal vectors are used to describe the model, all calculations are done in three dimensions making calculation of the real-space structure trivial. Using this model, it is shown that both...... the positions of the bragg reflections and information about the relative intensities of these reflections are in full accordance with the diffraction patterns reported for microcrystals of the rapidly quenched Al86Mn14 alloy. It is also shown that at least the local structure possesses full icosahedral...

A Model Describing Stable Coherent Synchrotron Radiation in Storage Rings

International Nuclear Information System (INIS)

Sannibale, F.

2004-01-01

We present a model describing high power stable broadband coherent synchrotron radiation (CSR) in the terahertz frequency region in an electron storage ring. The model includes distortion of bunch shape from the synchrotron radiation (SR), which enhances higher frequency coherent emission, and limits to stable emission due to an instability excited by the SR wakefield. It gives a quantitative explanation of several features of the recent observations of CSR at the BESSY II storage ring. We also use this model to optimize the performance of a source for stable CSR emission

A model describing stable coherent synchrotron radiation in storage rings

International Nuclear Information System (INIS)

Sannibale, F.; Byrd, J.M.; Loftsdottir, A.; Venturini, M.; Abo-Bakr, M.; Feikes, J.; Holldack, K.; Kuske, P.; Wuestefeld, G.; Huebers, H.-W.; Warnock, R.

2004-01-01

We present a model describing high power stable broadband coherent synchrotron radiation (CSR) in the terahertz frequency region in an electron storage ring. The model includes distortion of bunch shape from the synchrotron radiation (SR), which enhances higher frequency coherent emission, and limits to stable emission due to an instability excited by the SR wakefield. It gives a quantitative explanation of several features of the recent observations of CSR at the BESSY II storage ring. We also use this model to optimize the performance of a source for stable CSR emission

Multi-parameter variational calculations for the (2+1)-dimensional U(1) lattice gauge theory and the XY model

International Nuclear Information System (INIS)

Heys, D.W.; Stump, D.R.

1987-01-01

Variational calculations are described that use multi-parameter trial wave functions for the U(1) lattice gauge theory in two space dimensions, and for the XY model. The trial functions are constructed as the exponential of a linear combination of states from the strong-coupling basis of the model, with the coefficients treated as variational parameters. The expectation of the hamiltonian is computed by the Monte Carlo method, using a reweighting technique to evaluate expectation values in finite patches of the parameter space. The trial function for the U(1) gauge theory involves six variational parameters, and its weak-coupling behaviour is in reasonable agreement with theoretical expectations. (orig.)

Parameter Estimation of Nonlinear Models in Forestry.

OpenAIRE

Fekedulegn, Desta; Mac Siúrtáin, Máirtín Pádraig; Colbert, Jim J.

1999-01-01

Partial derivatives of the negative exponential, monomolecular, Mitcherlich, Gompertz, logistic, Chapman-Richards, von Bertalanffy, Weibull and the Richard’s nonlinear growth models are presented. The application of these partial derivatives in estimating the model parameters is illustrated. The parameters are estimated using the Marquardt iterative method of nonlinear regression relating top height to age of Norway spruce (Picea abies L.) from the Bowmont Norway Spruce Thinnin...

Learning about physical parameters: the importance of model discrepancy

International Nuclear Information System (INIS)

Brynjarsdóttir, Jenný; O'Hagan, Anthony

2014-01-01

Science-based simulation models are widely used to predict the behavior of complex physical systems. It is also common to use observations of the physical system to solve the inverse problem, that is, to learn about the values of parameters within the model, a process which is often called calibration. The main goal of calibration is usually to improve the predictive performance of the simulator but the values of the parameters in the model may also be of intrinsic scientific interest in their own right. In order to make appropriate use of observations of the physical system it is important to recognize model discrepancy, the difference between reality and the simulator output. We illustrate through a simple example that an analysis that does not account for model discrepancy may lead to biased and over-confident parameter estimates and predictions. The challenge with incorporating model discrepancy in statistical inverse problems is being confounded with calibration parameters, which will only be resolved with meaningful priors. For our simple example, we model the model-discrepancy via a Gaussian process and demonstrate that through accounting for model discrepancy our prediction within the range of data is correct. However, only with realistic priors on the model discrepancy do we uncover the true parameter values. Through theoretical arguments we show that these findings are typical of the general problem of learning about physical parameters and the underlying physical system using science-based mechanistic models. (paper)

Wind Farm Decentralized Dynamic Modeling With Parameters

DEFF Research Database (Denmark)

Soltani, Mohsen; Shakeri, Sayyed Mojtaba; Grunnet, Jacob Deleuran

2010-01-01

Development of dynamic wind flow models for wind farms is part of the research in European research FP7 project AEOLUS. The objective of this report is to provide decentralized dynamic wind flow models with parameters. The report presents a structure for decentralized flow models with inputs from...... local models. The results of this report are especially useful, but not limited, to design a decentralized wind farm controller, since in centralized controller design one can also use the model and update it in a central computing node.......Development of dynamic wind flow models for wind farms is part of the research in European research FP7 project AEOLUS. The objective of this report is to provide decentralized dynamic wind flow models with parameters. The report presents a structure for decentralized flow models with inputs from...

Optimizing incomplete sample designs for item response model parameters

NARCIS (Netherlands)

van der Linden, Willem J.

Several models for optimizing incomplete sample designs with respect to information on the item parameters are presented. The following cases are considered: (1) known ability parameters; (2) unknown ability parameters; (3) item sets with multiple ability scales; and (4) response models with

A model describing water and salt migration in concrete during wetting/drying cycles

NARCIS (Netherlands)

Arends, T.; Taher, A.; van der Zanden, A.J.J.; Brouwers, H.J.H.; Bilek, V.; Kersner, Z.

2014-01-01

In order to predict the life span of concrete structures, models describing the migration of chloride are needed. In this paper, a start is made with a simple, theoretical model describing water and chloride transport in a concrete sample. First, transport of water in concrete is considered with

Recommended Parameter Values for GENII Modeling of Radionuclides in Routine Air and Water Releases

Energy Technology Data Exchange (ETDEWEB)

Snyder, Sandra F. [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Arimescu, Carmen [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Napier, Bruce A. [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Hay, Tristan R. [Pacific Northwest National Lab. (PNNL), Richland, WA (United States)

2012-11-01

The GENII v2 code is used to estimate dose to individuals or populations from the release of radioactive materials into air or water. Numerous parameter values are required for input into this code. User-defined parameters cover the spectrum from chemical data, meteorological data, agricultural data, and behavioral data. This document is a summary of parameter values that reflect conditions in the United States. Reasonable regional and age-dependent data is summarized. Data availability and quality varies. The set of parameters described address scenarios for chronic air emissions or chronic releases to public waterways. Considerations for the special tritium and carbon-14 models are briefly addressed. GENIIv2.10.0 is the current software version that this document supports.

The methodology of choice Cam-Clay model parameters for loess subsoil

Science.gov (United States)

Nepelski, Krzysztof; Błazik-Borowa, Ewa

2018-01-01

The paper deals with the calibration method of FEM subsoil model described by the constitutive Cam-Clay model. The four-storey residential building and solid substrate are modelled. Identification of the substrate is made using research drilling, CPT static tests, DMT Marchetti dilatometer, and laboratory tests. Latter are performed on the intact soil specimens which are taken from the wide planning trench at the depth of foundation. The real building settlements was measured as the vertical displacement of benchmarks. These measurements were carried out periodically during the erection of the building and its operation. Initially, the Cam Clay model parameters were determined on the basis of the laboratory tests, and later, they were corrected by taking into consideration numerical analyses results (whole building and its parts) and real building settlements.

Ordering dynamics of microscopic models with nonconserved order parameter of continuous symmetry

DEFF Research Database (Denmark)

Zhang, Z.; Mouritsen, Ole G.; Zuckermann, Martin J.

1993-01-01

crystals. For both models, which have a nonconserved order parameter, it is found that the linear scale, R(t), of the evolving order, following quenches to below the transition temperature, grows at late times in an effectively algebraic fashion, R(t)∼tn, with exponent values which are strongly temperature......Numerical Monte Carlo temperature-quenching experiments have been performed on two three-dimensional classical lattice models with continuous ordering symmetry: the Lebwohl-Lasher model [Phys. Rev. A 6, 426 (1972)] and the ferromagnetic isotropic Heisenberg model. Both models describe a transition...... from a disordered phase to an orientationally ordered phase of continuous symmetry. The Lebwohl-Lasher model accounts for the orientational ordering properties of the nematic-isotropic transition in liquid crystals and the Heisenberg model for the ferromagnetic-paramagnetic transition in magnetic...

Parameters and error of a theoretical model

International Nuclear Information System (INIS)

Moeller, P.; Nix, J.R.; Swiatecki, W.

1986-09-01

We propose a definition for the error of a theoretical model of the type whose parameters are determined from adjustment to experimental data. By applying a standard statistical method, the maximum-likelihoodlmethod, we derive expressions for both the parameters of the theoretical model and its error. We investigate the derived equations by solving them for simulated experimental and theoretical quantities generated by use of random number generators. 2 refs., 4 tabs

Application of multi-parameter chorus and plasmaspheric hiss wave models in radiation belt modeling

Science.gov (United States)

Aryan, H.; Kang, S. B.; Balikhin, M. A.; Fok, M. C. H.; Agapitov, O. V.; Komar, C. M.; Kanekal, S. G.; Nagai, T.; Sibeck, D. G.

2017-12-01

Numerical simulation studies of the Earth's radiation belts are important to understand the acceleration and loss of energetic electrons. The Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model along with many other radiation belt models require inputs for pitch angle, energy, and cross diffusion of electrons, due to chorus and plasmaspheric hiss waves. These parameters are calculated using statistical wave distribution models of chorus and plasmaspheric hiss amplitudes. In this study we incorporate recently developed multi-parameter chorus and plasmaspheric hiss wave models based on geomagnetic index and solar wind parameters. We perform CIMI simulations for two geomagnetic storms and compare the flux enhancement of MeV electrons with data from the Van Allen Probes and Akebono satellites. We show that the relativistic electron fluxes calculated with multi-parameter wave models resembles the observations more accurately than the relativistic electron fluxes calculated with single-parameter wave models. This indicates that wave models based on a combination of geomagnetic index and solar wind parameters are more effective as inputs to radiation belt models.

A Markov chain model for N-linked protein glycosylation – towards a low-parameter tool for model-driven glycoengineering

DEFF Research Database (Denmark)

Spahn, Philipp N.; Hansen, Anders Holmgaard; Hansen, Henning Gram

2016-01-01

Glycosylation is a critical quality attribute of most recombinant biotherapeutics. Consequently, drug development requires careful control of glycoforms to meet bioactivity and biosafety requirements. However, glycoengineering can be extraordinarily difficult given the complex reaction networks...... present a novel low-parameter approach to describe glycosylation using flux-balance and Markov chain modeling. The model recapitulates the biological complexity of glycosylation, but does not require user-provided kinetic information. We use this method to predict and experimentally validate glycoprofiles...

WINKLER'S SINGLE-PARAMETER SUBGRADE MODEL FROM ...

African Journals Online (AJOL)

Preferred Customer

Page 1 ... corresponding single-parameter Winkler model presented in this work. Keywords: Heterogeneous subgrade, Reissner's simplified continuum, Shear interaction, Simplified continuum, Winkler ... model in practical applications and its long time familiarity among practical engineers, its usage has endured to this date ...

The method for determination of parameters of the phenomenological continual model of soil organic matter transformation

Directory of Open Access Journals (Sweden)

S. I. Bartsev

2015-06-01

Full Text Available A possible method for experimental determination of parameters of the previously proposed continual mathematical model of soil organic matter transformation is theoretically considered in this paper. The previously proposed by the authors continual model of soil organic matter transformation, based on using the rate of matter transformation as a continual scale of its recalcitrance, describes the transformation process phenomenologically without going into detail of microbiological mechanisms of transformation. Thereby simplicity of the model is achieved. The model is represented in form of one differential equation in first­order partial derivatives, which has an analytical solution in elementary functions. The model equation contains a small number of empirical parameters which generally characterize environmental conditions where the matter transformation process occurs and initial properties of the plant litter. Given the values of these parameters, it is possible to calculate dynamics of soil organic matter stocks and its distribution over transformation rate. In the present study, possible approaches for determination of the model parameters are considered and a simple method of their experimental measurement is proposed. An experiment of an incubation of chemically homogeneous samples in soil and multiple sequential measurement of the sample mass loss with time is proposed. An equation of time dynamics of mass loss of incubated homogeneous sample is derived from the basic assumption of the presented soil organic matter transformation model. Thus, fitting by the least squares method the parameters of sample mass loss curve calculated according the proposed mass loss dynamics equation allows to determine the parameters of the general equation of soil organic transformation model.

Insight into model mechanisms through automatic parameter fitting: a new methodological framework for model development.

Science.gov (United States)

Tøndel, Kristin; Niederer, Steven A; Land, Sander; Smith, Nicolas P

2014-05-20

Striking a balance between the degree of model complexity and parameter identifiability, while still producing biologically feasible simulations using modelling is a major challenge in computational biology. While these two elements of model development are closely coupled, parameter fitting from measured data and analysis of model mechanisms have traditionally been performed separately and sequentially. This process produces potential mismatches between model and data complexities that can compromise the ability of computational frameworks to reveal mechanistic insights or predict new behaviour. In this study we address this issue by presenting a generic framework for combined model parameterisation, comparison of model alternatives and analysis of model mechanisms. The presented methodology is based on a combination of multivariate metamodelling (statistical approximation of the input-output relationships of deterministic models) and a systematic zooming into biologically feasible regions of the parameter space by iterative generation of new experimental designs and look-up of simulations in the proximity of the measured data. The parameter fitting pipeline includes an implicit sensitivity analysis and analysis of parameter identifiability, making it suitable for testing hypotheses for model reduction. Using this approach, under-constrained model parameters, as well as the coupling between parameters within the model are identified. The methodology is demonstrated by refitting the parameters of a published model of cardiac cellular mechanics using a combination of measured data and synthetic data from an alternative model of the same system. Using this approach, reduced models with simplified expressions for the tropomyosin/crossbridge kinetics were found by identification of model components that can be omitted without affecting the fit to the parameterising data. Our analysis revealed that model parameters could be constrained to a standard deviation of on

Modelling of intermittent microwave convective drying: parameter sensitivity

Directory of Open Access Journals (Sweden)

Zhang Zhijun

2017-06-01

Full Text Available The reliability of the predictions of a mathematical model is a prerequisite to its utilization. A multiphase porous media model of intermittent microwave convective drying is developed based on the literature. The model considers the liquid water, gas and solid matrix inside of food. The model is simulated by COMSOL software. Its sensitivity parameter is analysed by changing the parameter values by ±20%, with the exception of several parameters. The sensitivity analysis of the process of the microwave power level shows that each parameter: ambient temperature, effective gas diffusivity, and evaporation rate constant, has significant effects on the process. However, the surface mass, heat transfer coefficient, relative and intrinsic permeability of the gas, and capillary diffusivity of water do not have a considerable effect. The evaporation rate constant has minimal parameter sensitivity with a ±20% value change, until it is changed 10-fold. In all results, the temperature and vapour pressure curves show the same trends as the moisture content curve. However, the water saturation at the medium surface and in the centre show different results. Vapour transfer is the major mass transfer phenomenon that affects the drying process.

Endorsement of Models Describing Sexual Response of Men and Women with a Sexual Partner

DEFF Research Database (Denmark)

Giraldi, Annamaria; Kristensen, Ellids; Sand, Michael

2015-01-01

, erectile dysfunction and dissatisfaction with sexual life were significantly related to endorsement of the Basson model or none of the models (P = 0.01). CONCLUSIONS: No single model of sexual response could describe men's and women's sexual responses. The majority of men and women with no sexual......INTRODUCTION: Several models have been used to describe men's and women's sexual responses. These models have been conceptualized as linear or circular models. The circular models were proposed to describe women's sexual function best. AIM: This study aims to determine whether men and women thought...... that current theoretical models of sexual responses accurately reflected their own sexual experience and to what extent this was influenced by sexual dysfunction. METHODS: A cross-sectional study of a large, broadly sampled, nonclinical population, cohort of Danish men and women. The Female Sexual Function...

Parameter identification in ODE models with oscillatory dynamics: a Fourier regularization approach

Science.gov (United States)

Chiara D'Autilia, Maria; Sgura, Ivonne; Bozzini, Benedetto

2017-12-01

In this paper we consider a parameter identification problem (PIP) for data oscillating in time, that can be described in terms of the dynamics of some ordinary differential equation (ODE) model, resulting in an optimization problem constrained by the ODEs. In problems with this type of data structure, simple application of the direct method of control theory (discretize-then-optimize) yields a least-squares cost function exhibiting multiple ‘low’ minima. Since in this situation any optimization algorithm is liable to fail in the approximation of a good solution, here we propose a Fourier regularization approach that is able to identify an iso-frequency manifold {{ S}} of codimension-one in the parameter space \

Parameter identification in multinomial processing tree models

NARCIS (Netherlands)

Schmittmann, V.D.; Dolan, C.V.; Raijmakers, M.E.J.; Batchelder, W.H.

2010-01-01

Multinomial processing tree models form a popular class of statistical models for categorical data that have applications in various areas of psychological research. As in all statistical models, establishing which parameters are identified is necessary for model inference and selection on the basis

Pharmacodynamic Model To Describe the Concentration-Dependent Selection of Cefotaxime-Resistant Escherichia coli

Science.gov (United States)

Olofsson, Sara K.; Geli, Patricia; Andersson, Dan I.; Cars, Otto

2005-01-01

Antibiotic dosing regimens may vary in their capacity to select mutants. Our hypothesis was that selection of a more resistant bacterial subpopulation would increase with the time within a selective window (SW), i.e., when drug concentrations fall between the MICs of two strains. An in vitro kinetic model was used to study the selection of two Escherichia coli strains with different susceptibilities to cefotaxime. The bacterial mixtures were exposed to cefotaxime for 24 h and SWs of 1, 2, 4, 8, and 12 h. A mathematical model was developed that described the selection of preexisting and newborn mutants and the post-MIC effect (PME) as functions of pharmacokinetic parameters. Our main conclusions were as follows: (i) the selection between preexisting mutants increased with the time within the SW; (ii) the emergence and selection of newborn mutants increased with the time within the SW (with a short time, only 4% of the preexisting mutants were replaced by newborn mutants, compared to the longest times, where 100% were replaced); and (iii) PME increased with the area under the concentration-time curve (AUC) and was slightly more pronounced with a long elimination half-life (T1/2) than with a short T1/2 situation, when AUC is fixed. We showed that, in a dynamic competition between strains with different levels of resistance, the appearance of newborn high-level resistant mutants from the parental strains and the PME can strongly affect the outcome of the selection and that pharmacodynamic models can be used to predict the outcome of resistance development. PMID:16304176

An inverse modeling approach to estimate groundwater flow and transport model parameters at a research site at Vandenberg AFB, CA

Science.gov (United States)

Rasa, E.; Foglia, L.; Mackay, D. M.; Ginn, T. R.; Scow, K. M.

2009-12-01

A numerical groundwater fate and transport model was developed for analyses of data from field experiments evaluating the impacts of ethanol on the natural attenuation of benzene, toluene, ethylbenzene, and xylenes (BTEX) and methyl tert-butyl ether (MTBE) at Vandenberg Air Force Base, Site 60. We used the U.S. Geological Survey (USGS) groundwater flow (MODFLOW2000) and transport (MT3DMS) models in conjunction with the USGS universal inverse modeling code (UCODE) to jointly determine flow and transport parameters using bromide tracer data from multiple experiments in the same location. The key flow and transport parameters include hydraulic conductivity of aquifer and aquitard layers, porosity, and transverse and longitudinal dispersivity. Aquifer and aquitard layers were assumed homogenous in this study. Therefore, the calibration parameters were not spatially variable within each layer. A total of 162 monitoring wells in seven transects perpendicular to the mean flow direction were monitored over the course of ten months, resulting in 1,766 bromide concentration data points and 149 head values used as observations for the inverse modeling. The results showed the significance of the concentration observation data in predicting the flow model parameters and indicated the sensitivity of the hydraulic conductivity of different zones in the aquifer including the excavated former contaminant zone. The model has already been used to evaluate alternative designs for further experiments on in situ bioremediation of the tert-butyl alcohol (TBA) plume remaining at the site. We describe the recent applications of the model and future work, including adding reaction submodels to the calibrated flow model.

Methodology for identifying parameters for the TRNSYS model Type 210 - wood pellet stoves and boilers

Energy Technology Data Exchange (ETDEWEB)

Persson, Tomas; Fiedler, Frank; Nordlander, Svante

2006-05-15

This report describes a method how to perform measurements on boilers and stoves and how to identify parameters from the measurements for the boiler/stove-model TRNSYS Type 210. The model can be used for detailed annual system simulations using TRNSYS. Experience from measurements on three different pellet stoves and four boilers were used to develop this methodology. Recommendations for the set up of measurements are given and the required combustion theory for the data evaluation and data preparation are given. The data evaluation showed that the uncertainties are quite large for the measured flue gas flow rate and for boilers and stoves with high fraction of energy going to the water jacket also the calculated heat rate to the room may have large uncertainties. A methodology for the parameter identification process and identified parameters for two different stoves and three boilers are given. Finally the identified models are compared with measured data showing that the model generally agreed well with measured data during both stationary and dynamic conditions.

Reduction of the number of parameters needed for a polynomial random regression test-day model

NARCIS (Netherlands)

Pool, M.H.; Meuwissen, T.H.E.

2000-01-01

Legendre polynomials were used to describe the (co)variance matrix within a random regression test day model. The goodness of fit depended on the polynomial order of fit, i.e., number of parameters to be estimated per animal but is limited by computing capacity. Two aspects: incomplete lactation

Optimal parameters for the FFA-Beddoes dynamic stall model

Energy Technology Data Exchange (ETDEWEB)

Bjoerck, A; Mert, M [FFA, The Aeronautical Research Institute of Sweden, Bromma (Sweden); Madsen, H A [Risoe National Lab., Roskilde (Denmark)

1999-03-01

Unsteady aerodynamic effects, like dynamic stall, must be considered in calculation of dynamic forces for wind turbines. Models incorporated in aero-elastic programs are of semi-empirical nature. Resulting aerodynamic forces therefore depend on values used for the semi-empiricial parameters. In this paper a study of finding appropriate parameters to use with the Beddoes-Leishman model is discussed. Minimisation of the `tracking error` between results from 2D wind tunnel tests and simulation with the model is used to find optimum values for the parameters. The resulting optimum parameters show a large variation from case to case. Using these different sets of optimum parameters in the calculation of blade vibrations, give rise to quite different predictions of aerodynamic damping which is discussed. (au)

A parameter identification problem arising from a two-dimensional airfoil section model

International Nuclear Information System (INIS)

Cerezo, G.M.

1994-01-01

The development of state space models for aeroelastic systems, including unsteady aerodynamics, is particularly important for the design of highly maneuverable aircraft. In this work we present a state space formulation for a special class of singular neutral functional differential equations (SNFDE) with initial data in C(-1, 0). This work is motivated by the two-dimensional airfoil model presented by Burns, Cliff and Herdman in. In the same authors discuss the validity of the assumptions under which the model was formulated. They pay special attention to the derivation of the evolution equation for the circulation on the airfoil. This equation was coupled to the rigid-body dynamics of the airfoil in order to obtain a complete set of functional differential equations that describes the composite system. The resulting mathematical model for the aeroelastic system has a weakly singular component. In this work we consider a finite delay approximation to the model presented in. We work with a scalar model in which we consider the weak singularity appearing in the original problem. The main goal of this work is to develop numerical techniques for the identification of the parameters appearing in the kernel of the associated scalar integral equation. Clearly this is the first step in the study of parameter identification for the original model and the corresponding validation of this model for the aeroelastic system

Modelling of nonhomogeneous atmosphere in NPP containment using lumped-parameter model based on CFD calculations

International Nuclear Information System (INIS)

Kljenak, I.; Mavko, B.; Babic, M.

2005-01-01

steam, air and helium were injected during different phases of the experiment. The thermal-hydraulic behaviour was determined by the dominant physical phenomena: gas injection, steam condensation, heat transfer and buoyant flow. During certain phases of the experiment, intermediate steady states were reached when the steam condensation rate became equal to the steam injection rate, with boundary conditions remaining constant. In the proposed work, steady states were simulated independently by applying the approach described above. The code CFX4.4 was used as a CFD code, whereas the code CONTAIN was used as a lumped-parameter code. The simulation with the CFX code was performed with a two-dimensional axisymmetric model of the TOSQAN vessel. The proposed work may be considered as a first step towards development of a suitable methodology for modelling nonhomogeneous atmosphere in NPP containments with lumped-parameter codes. (authors)

A model to describe the performance of the UASB reactor.

Science.gov (United States)

Rodríguez-Gómez, Raúl; Renman, Gunno; Moreno, Luis; Liu, Longcheng

2014-04-01

A dynamic model to describe the performance of the Upflow Anaerobic Sludge Blanket (UASB) reactor was developed. It includes dispersion, advection, and reaction terms, as well as the resistances through which the substrate passes before its biotransformation. The UASB reactor is viewed as several continuous stirred tank reactors connected in series. The good agreement between experimental and simulated results shows that the model is able to predict the performance of the UASB reactor (i.e. substrate concentration, biomass concentration, granule size, and height of the sludge bed).

Lumped-parameters equivalent circuit for condenser microphones modeling.

Science.gov (United States)

Esteves, Josué; Rufer, Libor; Ekeom, Didace; Basrour, Skandar

2017-10-01

This work presents a lumped parameters equivalent model of condenser microphone based on analogies between acoustic, mechanical, fluidic, and electrical domains. Parameters of the model were determined mainly through analytical relations and/or finite element method (FEM) simulations. Special attention was paid to the air gap modeling and to the use of proper boundary condition. Corresponding lumped-parameters were obtained as results of FEM simulations. Because of its simplicity, the model allows a fast simulation and is readily usable for microphone design. This work shows the validation of the equivalent circuit on three real cases of capacitive microphones, including both traditional and Micro-Electro-Mechanical Systems structures. In all cases, it has been demonstrated that the sensitivity and other related data obtained from the equivalent circuit are in very good agreement with available measurement data.

Analysis of Modeling Parameters on Threaded Screws.

Energy Technology Data Exchange (ETDEWEB)

Vigil, Miquela S. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Brake, Matthew Robert [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Vangoethem, Douglas [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2015-06-01

Assembled mechanical systems often contain a large number of bolted connections. These bolted connections (joints) are integral aspects of the load path for structural dynamics, and, consequently, are paramount for calculating a structure's stiffness and energy dissipation prop- erties. However, analysts have not found the optimal method to model appropriately these bolted joints. The complexity of the screw geometry cause issues when generating a mesh of the model. This paper will explore different approaches to model a screw-substrate connec- tion. Model parameters such as mesh continuity, node alignment, wedge angles, and thread to body element size ratios are examined. The results of this study will give analysts a better understanding of the influences of these parameters and will aide in finding the optimal method to model bolted connections.

Parameter Estimates in Differential Equation Models for Chemical Kinetics

Science.gov (United States)

Winkel, Brian

2011-01-01

We discuss the need for devoting time in differential equations courses to modelling and the completion of the modelling process with efforts to estimate the parameters in the models using data. We estimate the parameters present in several differential equation models of chemical reactions of order n, where n = 0, 1, 2, and apply more general…

Simultaneous inference for model averaging of derived parameters

DEFF Research Database (Denmark)

Jensen, Signe Marie; Ritz, Christian

2015-01-01

Model averaging is a useful approach for capturing uncertainty due to model selection. Currently, this uncertainty is often quantified by means of approximations that do not easily extend to simultaneous inference. Moreover, in practice there is a need for both model averaging and simultaneous...... inference for derived parameters calculated in an after-fitting step. We propose a method for obtaining asymptotically correct standard errors for one or several model-averaged estimates of derived parameters and for obtaining simultaneous confidence intervals that asymptotically control the family...

Parameter identification in the logistic STAR model

DEFF Research Database (Denmark)

Ekner, Line Elvstrøm; Nejstgaard, Emil

We propose a new and simple parametrization of the so-called speed of transition parameter of the logistic smooth transition autoregressive (LSTAR) model. The new parametrization highlights that a consequence of the well-known identification problem of the speed of transition parameter is that th...

Computer modeling describes gravity-related adaptation in cell cultures.

Science.gov (United States)

Alexandrov, Ludmil B; Alexandrova, Stoyana; Usheva, Anny

2009-12-16

Questions about the changes of biological systems in response to hostile environmental factors are important but not easy to answer. Often, the traditional description with differential equations is difficult due to the overwhelming complexity of the living systems. Another way to describe complex systems is by simulating them with phenomenological models such as the well-known evolutionary agent-based model (EABM). Here we developed an EABM to simulate cell colonies as a multi-agent system that adapts to hyper-gravity in starvation conditions. In the model, the cell's heritable characteristics are generated and transferred randomly to offspring cells. After a qualitative validation of the model at normal gravity, we simulate cellular growth in hyper-gravity conditions. The obtained data are consistent with previously confirmed theoretical and experimental findings for bacterial behavior in environmental changes, including the experimental data from the microgravity Atlantis and the Hypergravity 3000 experiments. Our results demonstrate that it is possible to utilize an EABM with realistic qualitative description to examine the effects of hypergravity and starvation on complex cellular entities.

Regionalization of SWAT Model Parameters for Use in Ungauged Watersheds

Directory of Open Access Journals (Sweden)

Indrajeet Chaubey

2010-11-01

Full Text Available There has been a steady shift towards modeling and model-based approaches as primary methods of assessing watershed response to hydrologic inputs and land management, and of quantifying watershed-wide best management practice (BMP effectiveness. Watershed models often require some degree of calibration and validation to achieve adequate watershed and therefore BMP representation. This is, however, only possible for gauged watersheds. There are many watersheds for which there are very little or no monitoring data available, thus the question as to whether it would be possible to extend and/or generalize model parameters obtained through calibration of gauged watersheds to ungauged watersheds within the same region. This study explored the possibility of developing regionalized model parameter sets for use in ungauged watersheds. The study evaluated two regionalization methods: global averaging, and regression-based parameters, on the SWAT model using data from priority watersheds in Arkansas. Resulting parameters were tested and model performance determined on three gauged watersheds. Nash-Sutcliffe efficiencies (NS for stream flow obtained using regression-based parameters (0.53–0.83 compared well with corresponding values obtained through model calibration (0.45–0.90. Model performance obtained using global averaged parameter values was also generally acceptable (0.4 ≤ NS ≤ 0.75. Results from this study indicate that regionalized parameter sets for the SWAT model can be obtained and used for making satisfactory hydrologic response predictions in ungauged watersheds.

Bayesian uncertainty analysis for complex systems biology models: emulation, global parameter searches and evaluation of gene functions.

Science.gov (United States)

Vernon, Ian; Liu, Junli; Goldstein, Michael; Rowe, James; Topping, Jen; Lindsey, Keith

2018-01-02

Many mathematical models have now been employed across every area of systems biology. These models increasingly involve large numbers of unknown parameters, have complex structure which can result in substantial evaluation time relative to the needs of the analysis, and need to be compared to observed data of various forms. The correct analysis of such models usually requires a global parameter search, over a high dimensional parameter space, that incorporates and respects the most important sources of uncertainty. This can be an extremely difficult task, but it is essential for any meaningful inference or prediction to be made about any biological system. It hence represents a fundamental challenge for the whole of systems biology. Bayesian statistical methodology for the uncertainty analysis of complex models is introduced, which is designed to address the high dimensional global parameter search problem. Bayesian emulators that mimic the systems biology model but which are extremely fast to evaluate are embeded within an iterative history match: an efficient method to search high dimensional spaces within a more formal statistical setting, while incorporating major sources of uncertainty. The approach is demonstrated via application to a model of hormonal crosstalk in Arabidopsis root development, which has 32 rate parameters, for which we identify the sets of rate parameter values that lead to acceptable matches between model output and observed trend data. The multiple insights into the model's structure that this analysis provides are discussed. The methodology is applied to a second related model, and the biological consequences of the resulting comparison, including the evaluation of gene functions, are described. Bayesian uncertainty analysis for complex models using both emulators and history matching is shown to be a powerful technique that can greatly aid the study of a large class of systems biology models. It both provides insight into model behaviour

Optimization of Saturn paraboloid magnetospheric field model parameters using Cassini equatorial magnetic field data

Directory of Open Access Journals (Sweden)

E. S. Belenkaya

2016-07-01

Full Text Available The paraboloid model of Saturn's magnetosphere describes the magnetic field as being due to the sum of contributions from the internal field of the planet, the ring current, and the tail current, all contained by surface currents inside a magnetopause boundary which is taken to be a paraboloid of revolution about the planet-Sun line. The parameters of the model have previously been determined by comparison with data from a few passes through Saturn's magnetosphere in compressed and expanded states, depending on the prevailing dynamic pressure of the solar wind. Here we significantly expand such comparisons through examination of Cassini magnetic field data from 18 near-equatorial passes that span wide ranges of local time, focusing on modelling the co-latitudinal field component that defines the magnetic flux passing through the equatorial plane. For 12 of these passes, spanning pre-dawn, via noon, to post-midnight, the spacecraft crossed the magnetopause during the pass, thus allowing an estimate of the concurrent subsolar radial distance of the magnetopause R1 to be made, considered to be the primary parameter defining the scale size of the system. The best-fit model parameters from these passes are then employed to determine how the parameters vary with R1, using least-squares linear fits, thus providing predictive model parameters for any value of R1 within the range. We show that the fits obtained using the linear approximation parameters are of the same order as those for the individually selected parameters. We also show that the magnetic flux mapping to the tail lobes in these models is generally in good accord with observations of the location of the open-closed field line boundary in Saturn's ionosphere, and the related position of the auroral oval. We then investigate the field data on six passes through the nightside magnetosphere, for which the spacecraft did not cross the magnetopause, such that in this case we compare the

Modeling and Parameter Estimation of a Small Wind Generation System

Directory of Open Access Journals (Sweden)

Carlos A. Ramírez Gómez

2013-11-01

Full Text Available The modeling and parameter estimation of a small wind generation system is presented in this paper. The system consists of a wind turbine, a permanent magnet synchronous generator, a three phase rectifier, and a direct current load. In order to estimate the parameters wind speed data are registered in a weather station located in the Fraternidad Campus at ITM. Wind speed data were applied to a reference model programed with PSIM software. From that simulation, variables were registered to estimate the parameters. The wind generation system model together with the estimated parameters is an excellent representation of the detailed model, but the estimated model offers a higher flexibility than the programed model in PSIM software.

Improving weather predictability by including land-surface model parameter uncertainty

Science.gov (United States)

Orth, Rene; Dutra, Emanuel; Pappenberger, Florian

2016-04-01

The land surface forms an important component of Earth system models and interacts nonlinearly with other parts such as ocean and atmosphere. To capture the complex and heterogenous hydrology of the land surface, land surface models include a large number of parameters impacting the coupling to other components of the Earth system model. Focusing on ECMWF's land-surface model HTESSEL we present in this study a comprehensive parameter sensitivity evaluation using multiple observational datasets in Europe. We select 6 poorly constrained effective parameters (surface runoff effective depth, skin conductivity, minimum stomatal resistance, maximum interception, soil moisture stress function shape, total soil depth) and explore their sensitivity to model outputs such as soil moisture, evapotranspiration and runoff using uncoupled simulations and coupled seasonal forecasts. Additionally we investigate the possibility to construct ensembles from the multiple land surface parameters. In the uncoupled runs we find that minimum stomatal resistance and total soil depth have the most influence on model performance. Forecast skill scores are moreover sensitive to the same parameters as HTESSEL performance in the uncoupled analysis. We demonstrate the robustness of our findings by comparing multiple best performing parameter sets and multiple randomly chosen parameter sets. We find better temperature and precipitation forecast skill with the best-performing parameter perturbations demonstrating representativeness of model performance across uncoupled (and hence less computationally demanding) and coupled settings. Finally, we construct ensemble forecasts from ensemble members derived with different best-performing parameterizations of HTESSEL. This incorporation of parameter uncertainty in the ensemble generation yields an increase in forecast skill, even beyond the skill of the default system. Orth, R., E. Dutra, and F. Pappenberger, 2016: Improving weather predictability by

Parameters-related uncertainty in modeling sugar cane yield with an agro-Land Surface Model

Science.gov (United States)

Valade, A.; Ciais, P.; Vuichard, N.; Viovy, N.; Ruget, F.; Gabrielle, B.

2012-12-01

Agro-Land Surface Models (agro-LSM) have been developed from the coupling of specific crop models and large-scale generic vegetation models. They aim at accounting for the spatial distribution and variability of energy, water and carbon fluxes within soil-vegetation-atmosphere continuum with a particular emphasis on how crop phenology and agricultural management practice influence the turbulent fluxes exchanged with the atmosphere, and the underlying water and carbon pools. A part of the uncertainty in these models is related to the many parameters included in the models' equations. In this study, we quantify the parameter-based uncertainty in the simulation of sugar cane biomass production with the agro-LSM ORCHIDEE-STICS on a multi-regional approach with data from sites in Australia, La Reunion and Brazil. First, the main source of uncertainty for the output variables NPP, GPP, and sensible heat flux (SH) is determined through a screening of the main parameters of the model on a multi-site basis leading to the selection of a subset of most sensitive parameters causing most of the uncertainty. In a second step, a sensitivity analysis is carried out on the parameters selected from the screening analysis at a regional scale. For this, a Monte-Carlo sampling method associated with the calculation of Partial Ranked Correlation Coefficients is used. First, we quantify the sensitivity of the output variables to individual input parameters on a regional scale for two regions of intensive sugar cane cultivation in Australia and Brazil. Then, we quantify the overall uncertainty in the simulation's outputs propagated from the uncertainty in the input parameters. Seven parameters are identified by the screening procedure as driving most of the uncertainty in the agro-LSM ORCHIDEE-STICS model output at all sites. These parameters control photosynthesis (optimal temperature of photosynthesis, optimal carboxylation rate), radiation interception (extinction coefficient), root

Hybrid artificial bee colony algorithm for parameter optimization of five-parameter bidirectional reflectance distribution function model.

Science.gov (United States)

Wang, Qianqian; Zhao, Jing; Gong, Yong; Hao, Qun; Peng, Zhong

2017-11-20

A hybrid artificial bee colony (ABC) algorithm inspired by the best-so-far solution and bacterial chemotaxis was introduced to optimize the parameters of the five-parameter bidirectional reflectance distribution function (BRDF) model. To verify the performance of the hybrid ABC algorithm, we measured BRDF of three kinds of samples and simulated the undetermined parameters of the five-parameter BRDF model using the hybrid ABC algorithm and the genetic algorithm, respectively. The experimental results demonstrate that the hybrid ABC algorithm outperforms the genetic algorithm in convergence speed, accuracy, and time efficiency under the same conditions.

SPOTting Model Parameters Using a Ready-Made Python Package.

Directory of Open Access Journals (Sweden)

Tobias Houska

Full Text Available The choice for specific parameter estimation methods is often more dependent on its availability than its performance. We developed SPOTPY (Statistical Parameter Optimization Tool, an open source python package containing a comprehensive set of methods typically used to calibrate, analyze and optimize parameters for a wide range of ecological models. SPOTPY currently contains eight widely used algorithms, 11 objective functions, and can sample from eight parameter distributions. SPOTPY has a model-independent structure and can be run in parallel from the workstation to large computation clusters using the Message Passing Interface (MPI. We tested SPOTPY in five different case studies to parameterize the Rosenbrock, Griewank and Ackley functions, a one-dimensional physically based soil moisture routine, where we searched for parameters of the van Genuchten-Mualem function and a calibration of a biogeochemistry model with different objective functions. The case studies reveal that the implemented SPOTPY methods can be used for any model with just a minimal amount of code for maximal power of parameter optimization. They further show the benefit of having one package at hand that includes number of well performing parameter search methods, since not every case study can be solved sufficiently with every algorithm or every objective function.

Hydrological model parameter dimensionality is a weak measure of prediction uncertainty

Science.gov (United States)

Pande, S.; Arkesteijn, L.; Savenije, H.; Bastidas, L. A.

2015-04-01

This paper shows that instability of hydrological system representation in response to different pieces of information and associated prediction uncertainty is a function of model complexity. After demonstrating the connection between unstable model representation and model complexity, complexity is analyzed in a step by step manner. This is done measuring differences between simulations of a model under different realizations of input forcings. Algorithms are then suggested to estimate model complexity. Model complexities of the two model structures, SAC-SMA (Sacramento Soil Moisture Accounting) and its simplified version SIXPAR (Six Parameter Model), are computed on resampled input data sets from basins that span across the continental US. The model complexities for SIXPAR are estimated for various parameter ranges. It is shown that complexity of SIXPAR increases with lower storage capacity and/or higher recession coefficients. Thus it is argued that a conceptually simple model structure, such as SIXPAR, can be more complex than an intuitively more complex model structure, such as SAC-SMA for certain parameter ranges. We therefore contend that magnitudes of feasible model parameters influence the complexity of the model selection problem just as parameter dimensionality (number of parameters) does and that parameter dimensionality is an incomplete indicator of stability of hydrological model selection and prediction problems.

Parameter identification and global sensitivity analysis of Xin'anjiang model using meta-modeling approach

Directory of Open Access Journals (Sweden)

Xiao-meng Song

2013-01-01

Full Text Available 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.

A two component model describing nucleon structure functions in the low-x region

Energy Technology Data Exchange (ETDEWEB)

Bugaev, E.V. [Institute for Nuclear Research of the Russian Academy of Sciences, 7a, 60th October Anniversary prospect, Moscow 117312 (Russian Federation); Mangazeev, B.V. [Irkutsk State University, 1, Karl Marx Street, Irkutsk 664003 (Russian Federation)

2009-12-15

A two component model describing the electromagnetic nucleon structure functions in the low-x region, based on generalized vector dominance and color dipole approaches is briefly described. The model operates with the mesons of rho-family having the mass spectrum of the form m{sub n}{sup 2}=m{sub r}ho{sup 2}(1+2n) and takes into account the nondiagonal transitions in meson-nucleon scattering. The special cut-off factors are introduced in the model, to exclude the gamma-qq-bar-V transitions in the case of narrow qq-bar-pairs. For the color dipole part of the model the well known FKS-parameterization is used.

The strong prognostic value of KELIM, a model-based parameter from CA 125 kinetics in ovarian cancer

DEFF Research Database (Denmark)

You, Benoit; Colomban, Olivier; Heywood, Mark

2013-01-01

Unexpected results were recently reported about the poor surrogacy of Gynecologic Cancer Intergroup (GCIG) defined CA-125 response in recurrent ovarian cancer (ROC) patients. Mathematical modeling may help describe CA-125 decline dynamically and discriminate prognostic kinetic parameters....

Temporal variation and scaling of parameters for a monthly hydrologic model

Science.gov (United States)

Deng, Chao; Liu, Pan; Wang, Dingbao; Wang, Weiguang

2018-03-01

The temporal variation of model parameters is affected by the catchment conditions and has a significant impact on hydrological simulation. This study aims to evaluate the seasonality and downscaling of model parameter across time scales based on monthly and mean annual water balance models with a common model framework. Two parameters of the monthly model, i.e., k and m, are assumed to be time-variant at different months. Based on the hydrological data set from 121 MOPEX catchments in the United States, we firstly analyzed the correlation between parameters (k and m) and catchment properties (NDVI and frequency of rainfall events, α). The results show that parameter k is positively correlated with NDVI or α, while the correlation is opposite for parameter m, indicating that precipitation and vegetation affect monthly water balance by controlling temporal variation of parameters k and m. The multiple linear regression is then used to fit the relationship between ε and the means and coefficient of variations of parameters k and m. Based on the empirical equation and the correlations between the time-variant parameters and NDVI, the mean annual parameter ε is downscaled to monthly k and m. The results show that it has lower NSEs than these from model with time-variant k and m being calibrated through SCE-UA, while for several study catchments, it has higher NSEs than that of the model with constant parameters. The proposed method is feasible and provides a useful tool for temporal scaling of model parameter.

Evaluation of physiological parameters and their influence on doses calculated from two alternative dosimetric models for the gastrointestinal tract

International Nuclear Information System (INIS)

Lessard, E.T.; Skrable, K.W.

1981-01-01

Two dosimetric models, the catenary compartmental model and the slug flow model are examined using three sets of physiological parameters. The impact of physiological parameters on the dosimetry of the tract is illustrated by comparing calculated maximum permissible daily activity ingestion rates for single, unabsorbed, particle emitting radionuclides with an effective energy term of unity. The conclusions drawn from this intercomparison of six different cases are: (1) Current dosimetric models which use physiological parameters described in this article do not significantly disagree, and (2) For the determination of average dose equivalent rates to segments of the tract due to chronic, long term ingestion of any radionuclide, the catenary compartmental model is a mathematically simpler approach. The catenary model in addition has certain advantages for the calculation of the photon dose contribution to one segment from cumulated activity (disintegrations) in another segment

An improved state-parameter analysis of ecosystem models using data assimilation

Science.gov (United States)

Chen, M.; Liu, S.; Tieszen, L.L.; Hollinger, D.Y.

2008-01-01

Much of the effort spent in developing data assimilation methods for carbon dynamics analysis has focused on estimating optimal values for either model parameters or state variables. The main weakness of estimating parameter values alone (i.e., without considering state variables) is that all errors from input, output, and model structure are attributed to model parameter uncertainties. On the other hand, the accuracy of estimating state variables may be lowered if the temporal evolution of parameter values is not incorporated. This research develops a smoothed ensemble Kalman filter (SEnKF) by combining ensemble Kalman filter with kernel smoothing technique. SEnKF has following characteristics: (1) to estimate simultaneously the model states and parameters through concatenating unknown parameters and state variables into a joint state vector; (2) to mitigate dramatic, sudden changes of parameter values in parameter sampling and parameter evolution process, and control narrowing of parameter variance which results in filter divergence through adjusting smoothing factor in kernel smoothing algorithm; (3) to assimilate recursively data into the model and thus detect possible time variation of parameters; and (4) to address properly various sources of uncertainties stemming from input, output and parameter uncertainties. The SEnKF is tested by assimilating observed fluxes of carbon dioxide and environmental driving factor data from an AmeriFlux forest station located near Howland, Maine, USA, into a partition eddy flux model. Our analysis demonstrates that model parameters, such as light use efficiency, respiration coefficients, minimum and optimum temperatures for photosynthetic activity, and others, are highly constrained by eddy flux data at daily-to-seasonal time scales. The SEnKF stabilizes parameter values quickly regardless of the initial values of the parameters. Potential ecosystem light use efficiency demonstrates a strong seasonality. Results show that the

A lumped-parameter model for cryo-adsorber hydrogen storage tank

Energy Technology Data Exchange (ETDEWEB)

Senthil Kumar, V.; Raghunathan, K. [India Science Lab, General Motors R and D, Creator Building, International Technology Park, Bangalore 560066 (India); Kumar, Sudarshan [Chemical and Environmental Sciences Lab, General Motors R and D, 30500 Mound Road, Warren, MI 48090 (United States)

2009-07-15

One of the primary requirements for commercialization of hydrogen fuel-cell vehicles is the on-board storage of hydrogen in sufficient quantities. On-board storage of hydrogen by adsorption on nano-porous adsorbents at around liquid nitrogen temperatures and moderate pressures is considered viable and competitive with other storage technologies: liquid hydrogen, compressed gas, and metallic or complex hydrides. The four cryo-adsorber fuel tank processes occur over different time scales: refueling over a few minutes, discharge over a few hours, dormancy over a few days, and venting over a few weeks. The slower processes i.e. discharge, dormancy and venting are expected to have negligible temperature gradients within the bed, and hence are amenable to a lumped-parameter analysis. Here we report a quasi-static lumped-parameter model for the cryo-adsorber fuel tank, and discuss the results for these slower processes. We also describe an alternative solution method for dormancy and venting based on the thermodynamic state description. (author)

HF Parameters of Induction Motor

Directory of Open Access Journals (Sweden)

M. N. Benallal

2017-09-01

Full Text Available This article describes the results of experimental studies of HF input and primary parameters. A simulation model in Matlab SimulinkTM of multiphase windings as ladder circuit of coils is developed. A method for determining the primary parameters of ladder equivalent circuits is presented.

Experimental investigation of statistical models describing distribution of counts

International Nuclear Information System (INIS)

Salma, I.; Zemplen-Papp, E.

1992-01-01

The binomial, Poisson and modified Poisson models which are used for describing the statistical nature of the distribution of counts are compared theoretically, and conclusions for application are considered. The validity of the Poisson and the modified Poisson statistical distribution for observing k events in a short time interval is investigated experimentally for various measuring times. The experiments to measure the influence of the significant radioactive decay were performed with 89 Y m (T 1/2 =16.06 s), using a multichannel analyser (4096 channels) in the multiscaling mode. According to the results, Poisson statistics describe the counting experiment for short measuring times (up to T=0.5T 1/2 ) and its application is recommended. However, analysis of the data demonstrated, with confidence, that for long measurements (T≥T 1/2 ) Poisson distribution is not valid and the modified Poisson function is preferable. The practical implications in calculating uncertainties and in optimizing the measuring time are discussed. Differences between the standard deviations evaluated on the basis of the Poisson and binomial models are especially significant for experiments with long measuring time (T/T 1/2 ≥2) and/or large detection efficiency (ε>0.30). Optimization of the measuring time for paired observations yields the same solution for either the binomial or the Poisson distribution. (orig.)

Use of multilevel modeling for determining optimal parameters of heat supply systems

Science.gov (United States)

Stennikov, V. A.; Barakhtenko, E. A.; Sokolov, D. V.

2017-07-01

The problem of finding optimal parameters of a heat-supply system (HSS) is in ensuring the required throughput capacity of a heat network by determining pipeline diameters and characteristics and location of pumping stations. Effective methods for solving this problem, i.e., the method of stepwise optimization based on the concept of dynamic programming and the method of multicircuit optimization, were proposed in the context of the hydraulic circuit theory developed at Melentiev Energy Systems Institute (Siberian Branch, Russian Academy of Sciences). These methods enable us to determine optimal parameters of various types of piping systems due to flexible adaptability of the calculation procedure to intricate nonlinear mathematical models describing features of used equipment items and methods of their construction and operation. The new and most significant results achieved in developing methodological support and software for finding optimal parameters of complex heat supply systems are presented: a new procedure for solving the problem based on multilevel decomposition of a heat network model that makes it possible to proceed from the initial problem to a set of interrelated, less cumbersome subproblems with reduced dimensionality; a new algorithm implementing the method of multicircuit optimization and focused on the calculation of a hierarchical model of a heat supply system; the SOSNA software system for determining optimum parameters of intricate heat-supply systems and implementing the developed methodological foundation. The proposed procedure and algorithm enable us to solve engineering problems of finding the optimal parameters of multicircuit heat supply systems having large (real) dimensionality, and are applied in solving urgent problems related to the optimal development and reconstruction of these systems. The developed methodological foundation and software can be used for designing heat supply systems in the Central and the Admiralty regions in

Conceptual hierarchical modeling to describe wetland plant community organization

Science.gov (United States)

Little, A.M.; Guntenspergen, G.R.; Allen, T.F.H.

2010-01-01

Using multivariate analysis, we created a hierarchical modeling process that describes how differently-scaled environmental factors interact to affect wetland-scale plant community organization in a system of small, isolated wetlands on Mount Desert Island, Maine. We followed the procedure: 1) delineate wetland groups using cluster analysis, 2) identify differently scaled environmental gradients using non-metric multidimensional scaling, 3) order gradient hierarchical levels according to spatiotem-poral scale of fluctuation, and 4) assemble hierarchical model using group relationships with ordination axes and post-hoc tests of environmental differences. Using this process, we determined 1) large wetland size and poor surface water chemistry led to the development of shrub fen wetland vegetation, 2) Sphagnum and water chemistry differences affected fen vs. marsh / sedge meadows status within small wetlands, and 3) small-scale hydrologic differences explained transitions between forested vs. non-forested and marsh vs. sedge meadow vegetation. This hierarchical modeling process can help explain how upper level contextual processes constrain biotic community response to lower-level environmental changes. It creates models with more nuanced spatiotemporal complexity than classification and regression tree procedures. Using this process, wetland scientists will be able to generate more generalizable theories of plant community organization, and useful management models. ?? Society of Wetland Scientists 2009.

Ground level enhancement (GLE) energy spectrum parameters model

Science.gov (United States)

Qin, G.; Wu, S.

2017-12-01

We study the ground level enhancement (GLE) events in solar cycle 23 with the four energy spectra parameters, the normalization parameter C, low-energy power-law slope γ 1, high-energy power-law slope γ 2, and break energy E0, obtained by Mewaldt et al. 2012 who fit the observations to the double power-law equation. we divide the GLEs into two groups, one with strong acceleration by interplanetary (IP) shocks and another one without strong acceleration according to the condition of solar eruptions. We next fit the four parameters with solar event conditions to get models of the parameters for the two groups of GLEs separately. So that we would establish a model of energy spectrum for GLEs for the future space weather prediction.

A New Model for Describing the Rheological Behavior of Heavy and Extra Heavy Crude Oils in the Presence of Nanoparticles

Directory of Open Access Journals (Sweden)

Esteban A. Taborda

2017-12-01

Full Text Available The present work proposes for the first time a mathematical model for describing the rheological behavior of heavy and extra-heavy crude oils in the presence of nanoparticles. This model results from the combination of two existing mathematical models. The first one applies to the rheology of pseudoplastic substances, i.e., the Herschel-Bulkley model. The second one was previously developed by our research group to model the rheology of suspensions, namely the modified Pal and Rhodes model. The proposed model is applied to heavy and extra heavy crude oils in the presence of nanoparticles, considering the effects of nanoparticles concentration and surface chemical nature, temperature, and crude oil type. All the experimental data evaluated exhibited compelling goodness of fitting, and the physical parameters in the model follow correlate well with variations in viscosity. The new model is dependent of share rate and opens new possibilities for phenomenologically understanding viscosity reduction in heavy crude by adding solid nanoparticles and favoring the scale-up in enhanced oil recovery (EOR and/or improved oil recovery (IOR process.

Parameter Identification of Piecewise Linear Plasticity Metal Models Used in Numerical Modeling of Structures Under Plastic Deformation and Failure

Directory of Open Access Journals (Sweden)

A. V. Shmeliov

2016-01-01

Full Text Available The article describes the models of metallic materials used in the calculation of deformation and destruction of engineering structures. The reliability of material models can adequately assess the strength characteristics of the designs of new technology in its designing and certification.The article deals with contingencies and true mechanical properties of materials and presents equations of their relationship. It notes that in the software systems mechanical characteristics of materials are given in the true sense.The paper considers the linear and exponential models of materials, their characteristics, and methods to implement them. It considers the models of Johnson-Cook Steinberg-Guinan, Zerilli-Armstrong, Cowper-Symonds, Gurson-Tvergaard that take into account the strain rate and temperature of the material. Describes their applications, advantages and disadvantages. Considers single- and multi-parameter criteria of materials fracture, the prospects for their use. Gives a rational justification for using a piecewise linear plasticity material model *MAT_PIECEWISE_LINEAR_PLASTICITY (024, LS-DYNA software package for the engineering industry, and presents its main parameters.A technique to identify parameters of piecewise linear plasticity metal material models has been developed. The technique consists of the stages, based on the equations of transition from the conventional stress and strain values to the true ones. Taking into consideration the stressstrain state in the neck of the sample is a distinctive feature of the technique.Tensile tests of the round material samples have been conducted. To test the developed technique in the software package ANSYS LS-DYNA PC have been made tensile sample modeling and results comparison to show high convergence.Further improvement of the technique can be achieved through the development of a statistical approach to the analysis of the results of a series of tests. This will allow a kind of

Global sensitivity analysis for identifying important parameters of nitrogen nitrification and denitrification under model uncertainty and scenario uncertainty

Science.gov (United States)

Chen, Zhuowei; Shi, Liangsheng; Ye, Ming; Zhu, Yan; Yang, Jinzhong

2018-06-01

Nitrogen reactive transport modeling is subject to uncertainty in model parameters, structures, and scenarios. By using a new variance-based global sensitivity analysis method, this paper identifies important parameters for nitrogen reactive transport with simultaneous consideration of these three uncertainties. A combination of three scenarios of soil temperature and two scenarios of soil moisture creates a total of six scenarios. Four alternative models describing the effect of soil temperature and moisture content are used to evaluate the reduction functions used for calculating actual reaction rates. The results show that for nitrogen reactive transport problem, parameter importance varies substantially among different models and scenarios. Denitrification and nitrification process is sensitive to soil moisture content status rather than to the moisture function parameter. Nitrification process becomes more important at low moisture content and low temperature. However, the changing importance of nitrification activity with respect to temperature change highly relies on the selected model. Model-averaging is suggested to assess the nitrification (or denitrification) contribution by reducing the possible model error. Despite the introduction of biochemical heterogeneity or not, fairly consistent parameter importance rank is obtained in this study: optimal denitrification rate (Kden) is the most important parameter; reference temperature (Tr) is more important than temperature coefficient (Q10); empirical constant in moisture response function (m) is the least important one. Vertical distribution of soil moisture but not temperature plays predominant role controlling nitrogen reaction. This study provides insight into the nitrogen reactive transport modeling and demonstrates an effective strategy of selecting the important parameters when future temperature and soil moisture carry uncertainties or when modelers face with multiple ways of establishing nitrogen

Comparison of Parameter Identification Techniques

Directory of Open Access Journals (Sweden)

Eder Rafael

2016-01-01

Full Text Available Model-based control of mechatronic systems requires excellent knowledge about the physical behavior of each component. For several types of components of a system, e.g. mechanical or electrical ones, the dynamic behavior can be described by means of a mathematic model consisting of a set of differential equations, difference equations and/or algebraic constraint equations. The knowledge of a realistic mathematic model and its parameter values is essential to represent the behaviour of a mechatronic system. Frequently it is hard or impossible to obtain all required values of the model parameters from the producer, so an appropriate parameter estimation technique is required to compute missing parameters. A manifold of parameter identification techniques can be found in the literature, but their suitability depends on the mathematic model. Previous work dealt with the automatic assembly of mathematical models of serial and parallel robots with drives and controllers within the dynamic multibody simulation code HOTINT as fully-fledged mechatronic simulation. Several parameters of such robot models were identified successfully by our embedded algorithm. The present work proposes an improved version of the identification algorithm with higher performance. The quality of the identified parameter values and the computation effort are compared with another standard technique.

Users' guide to system dynamics model describing Coho salmon survival in Olema Creek, Point Reyes National Seashore, Marin County, California

Science.gov (United States)

Woodward, Andrea; Torregrosa, Alicia; Madej, Mary Ann; Reichmuth, Michael; Fong, Darren

2014-01-01

The system dynamics model described in this report is the result of a collaboration between U.S. Geological Survey (USGS) scientists and National Park Service (NPS) San Francisco Bay Area Network (SFAN) staff, whose goal was to develop a methodology to integrate inventory and monitoring data to better understand ecosystem dynamics and trends using salmon in Olema Creek, Marin County, California, as an example case. The SFAN began monitoring multiple life stages of coho salmon (Oncorhynchus kisutch) in Olema Creek during 2003 (Carlisle and others, 2013), building on previous monitoring of spawning fish and redds. They initiated water-quality and habitat monitoring, and had access to flow and weather data from other sources. This system dynamics model of the freshwater portion of the coho salmon life cycle in Olema Creek integrated 8 years of existing monitoring data, literature values, and expert opinion to investigate potential factors limiting survival and production, identify data gaps, and improve monitoring and restoration prescriptions. A system dynamics model is particularly effective when (1) data are insufficient in time series length and/or measured parameters for a statistical or mechanistic model, and (2) the model must be easily accessible by users who are not modelers. These characteristics helped us meet the following overarching goals for this model: Summarize and synthesize NPS monitoring data with data and information from other sources to describe factors and processes affecting freshwater survival of coho salmon in Olema Creek. Provide a model that can be easily manipulated to experiment with alternative values of model parameters and novel scenarios of environmental drivers. Although the model describes the ecological dynamics of Olema Creek, these dynamics are structurally similar to numerous other coastal streams along the California coast that also contain anadromous fish populations. The model developed for Olema can be used, at least as a

A new mathematical model of gastrointestinal transit incorporating age- and gender-dependent physiological parameters

International Nuclear Information System (INIS)

Stubbs, J.B.

1992-01-01

As part of the revision by the International Commission on Radiological Protection (ICRP) of its report on Reference Man, an extensive review of the literature regarding anatomy and morphology of the gastrointestinal (GI) tract has been completed. Data on age- and gender-dependent GI physiology and motility may be included in the proposed ICRP report. A new mathematical model describing the transit of substances through the GI tract as well as the absorption and secretion of material in the GI tract has been developed. This mathematical description of GI tract kinetics utilizes more physiologically accurate transit processes than the mathematically simple, but nonphysiological, GI tract model that was used in ICRP Report 30. The proposed model uses a combination of zero- and first-order kinetics to describe motility. Some of the physiological parameters that the new model accounts for include sex, age, pathophysiological condition and meal phase (solid versus liquid). A computer algorithm, written in BASIC, based on this new model has been derived and results are compared to those of the ICRP-30 model

Rapid-relocation model for describing high-fluence retention of rare gases implanted in solids

Science.gov (United States)

Wittmaack, K.

2009-09-01

to be due to bombardment induced relocation and reemission, only the remaining 10% (or less) can be attributed to sputter erosion. The relocation efficiency is interpreted as the 'speed' of radiation enhanced diffusion towards the surface. The directionality of diffusion is attributed to the gradient of the defect density on the large-depth side of the damage distribution where most of the implanted rare gas atoms come to rest. Based on SRIM calculations, two representative parameters are defined, the peak number of lattice displacements, Nd,m, and the spacing, △ zr,d, between the peaks of the range and the damage distributions. Support in favour of rapid rare gas relocation by radiation enhanced diffusion is provided by the finding that the relocation efficiencies for Ar and Xe, which vary by up to one order of magnitude, scale as Ψ=kN/Δz, independent to the implantation energy (10-80 keV Ar, 10-500 keV Xe), within an error margin of only ± 15%. The parameter k contains the properties of the implanted rare gas atoms. A recently described computer simulation model, which assumed that the pressure established by the implanted gas drives reemission, is shown to reproduce measured Xe profiles quite well, but only at that energy at which the fitting parameter of the model was determined (140 keV). Using the same parameter at other energies, deviations by up to a factor of four are observed.

Voltage stability, bifurcation parameters and continuation methods

Energy Technology Data Exchange (ETDEWEB)

Alvarado, F L [Wisconsin Univ., Madison, WI (United States)

1994-12-31

This paper considers the importance of the choice of bifurcation parameter in the determination of the voltage stability limit and the maximum power load ability of a system. When the bifurcation parameter is power demand, the two limits are equivalent. However, when other types of load models and bifurcation parameters are considered, the two concepts differ. The continuation method is considered as a method for determination of voltage stability margins. Three variants of the continuation method are described: the continuation parameter is the bifurcation parameter the continuation parameter is initially the bifurcation parameter, but is free to change, and the continuation parameter is a new `arc length` parameter. Implementations of voltage stability software using continuation methods are described. (author) 23 refs., 9 figs.

Describing the processes of propagation and eliminating wildfires with the use of agent models

Directory of Open Access Journals (Sweden)

G. A. Dorrer

2017-10-01

Full Text Available A new method of describing the processes of propagation and elimination of wildfires on the basis of agent-based modeling is proposed. The main structural units of the creation of such models are the classes of active objects (agents. Agent approach, combined with Geographic Information Systems (GIS can effectively describe the interaction of a large number of participants in the process to combat wildfires: fire spreading, fire crews, mechanization, aerial means and other. In this paper we propose a multi-agent model to predict the spread of wildfire edge and simulate the direct method of extinguishing a ground fire with non-mechanized crews. The model consist with two classes of agents, designated A and B. The burning fire edge is represented as a chain of A-agents, each of which simulates the burning of an elementary portion of vegetation fuel. Fire front movement (moving the A-agent described by the Hamilton-Jacobi equation with using the indicatrises of normal front rate of spread (figurotris. The configuration of the front calculated on basis the algorithm of mobile grids. Agents other type, B-agents, described extinguishing process; they move to the agents of A type and act on them, reducing the combustion intensity to zero. Modeling system presented as two-level coloured nested Petri Net, which describes the agents’ interaction semantics. This model is implemented as a GIS-oriented software system that can be useful both in the fire fighting management as well as in staff training tactics to fighting wildfires. Some examples of modeling decision making on а ground fire extinguishing are presented.

Parameter estimation in nonlinear models for pesticide degradation

International Nuclear Information System (INIS)

Richter, O.; Pestemer, W.; Bunte, D.; Diekkrueger, B.

1991-01-01

A wide class of environmental transfer models is formulated as ordinary or partial differential equations. With the availability of fast computers, the numerical solution of large systems became feasible. The main difficulty in performing a realistic and convincing simulation of the fate of a substance in the biosphere is not the implementation of numerical techniques but rather the incomplete data basis for parameter estimation. Parameter estimation is a synonym for statistical and numerical procedures to derive reasonable numerical values for model parameters from data. The classical method is the familiar linear regression technique which dates back to the 18th century. Because it is easy to handle, linear regression has long been established as a convenient tool for analysing relationships. However, the wide use of linear regression has led to an overemphasis of linear relationships. In nature, most relationships are nonlinear and linearization often gives a poor approximation of reality. Furthermore, pure regression models are not capable to map the dynamics of a process. Therefore, realistic models involve the evolution in time (and space). This leads in a natural way to the formulation of differential equations. To establish the link between data and dynamical models, numerical advanced parameter identification methods have been developed in recent years. This paper demonstrates the application of these techniques to estimation problems in the field of pesticide dynamics. (7 refs., 5 figs., 2 tabs.)

Biological parameters for lung cancer in mathematical models of carcinogenesis

International Nuclear Information System (INIS)

Jacob, P.; Jacob, V.

2003-01-01

Applications of the two-step model of carcinogenesis with clonal expansion (TSCE) to lung cancer data are reviewed, including those on atomic bomb survivors from Hiroshima and Nagasaki, British doctors, Colorado Plateau miners, and Chinese tin miners. Different sets of identifiable model parameters are used in the literature. The parameter set which could be determined with the lowest uncertainty consists of the net proliferation rate gamma of intermediate cells, the hazard h 55 at an intermediate age, and the hazard H? at an asymptotically large age. Also, the values of these three parameters obtained in the various studies are more consistent than other identifiable combinations of the biological parameters. Based on representative results for these three parameters, implications for the biological parameters in the TSCE model are derived. (author)

Parameter sensitivity and uncertainty analysis for a storm surge and wave model

Directory of Open Access Journals (Sweden)

L. A. Bastidas

2016-09-01

Full Text Available Development and simulation of synthetic hurricane tracks is a common methodology used to estimate hurricane hazards in the absence of empirical coastal surge and wave observations. Such methods typically rely on numerical models to translate stochastically generated hurricane wind and pressure forcing into coastal surge and wave estimates. The model output uncertainty associated with selection of appropriate model parameters must therefore be addressed. The computational overburden of probabilistic surge hazard estimates is exacerbated by the high dimensionality of numerical surge and wave models. We present a model parameter sensitivity analysis of the Delft3D model for the simulation of hazards posed by Hurricane Bob (1991 utilizing three theoretical wind distributions (NWS23, modified Rankine, and Holland. The sensitive model parameters (of 11 total considered include wind drag, the depth-induced breaking γB, and the bottom roughness. Several parameters show no sensitivity (threshold depth, eddy viscosity, wave triad parameters, and depth-induced breaking αB and can therefore be excluded to reduce the computational overburden of probabilistic surge hazard estimates. The sensitive model parameters also demonstrate a large number of interactions between parameters and a nonlinear model response. While model outputs showed sensitivity to several parameters, the ability of these parameters to act as tuning parameters for calibration is somewhat limited as proper model calibration is strongly reliant on accurate wind and pressure forcing data. A comparison of the model performance with forcings from the different wind models is also presented.

Calibration of discrete element model parameters: soybeans

Science.gov (United States)

Ghodki, Bhupendra M.; Patel, Manish; Namdeo, Rohit; Carpenter, Gopal

2018-05-01

Discrete element method (DEM) simulations are broadly used to get an insight of flow characteristics of granular materials in complex particulate systems. DEM input parameters for a model are the critical prerequisite for an efficient simulation. Thus, the present investigation aims to determine DEM input parameters for Hertz-Mindlin model using soybeans as a granular material. To achieve this aim, widely acceptable calibration approach was used having standard box-type apparatus. Further, qualitative and quantitative findings such as particle profile, height of kernels retaining the acrylic wall, and angle of repose of experiments and numerical simulations were compared to get the parameters. The calibrated set of DEM input parameters includes the following (a) material properties: particle geometric mean diameter (6.24 mm); spherical shape; particle density (1220 kg m^{-3} ), and (b) interaction parameters such as particle-particle: coefficient of restitution (0.17); coefficient of static friction (0.26); coefficient of rolling friction (0.08), and particle-wall: coefficient of restitution (0.35); coefficient of static friction (0.30); coefficient of rolling friction (0.08). The results may adequately be used to simulate particle scale mechanics (grain commingling, flow/motion, forces, etc) of soybeans in post-harvest machinery and devices.

A consilience model to describe N2O production during biological N removal

DEFF Research Database (Denmark)

Domingo Felez, Carlos; Smets, Barth F.

2016-01-01

Nitrous oxide (N2O), a potent greenhouse gas, is produced during biological nitrogen conversion in wastewater treatment operations. Complex mechanisms underlie N2O production by autotrophic and heterotrophic organisms, which continue to be unravelled. Mathematical models that describe nitric oxide...... (NO) and N2O dynamics have been proposed. Here, a first comprehensive model that considers all relevant NO and N2O production and consumption mechanisms is proposed. The model describes autotrophic NO production by ammonia oxidizing bacteria associated with ammonia oxidation and with nitrite reduction......, followed by NO reduction to N2O. It also considers NO and N2O as intermediates in heterotrophic denitrification in a 4-step model. Three biological NO and N2O production pathways are accounted for, improving the capabilities of existing models while not increasing their complexity. Abiotic contributions...

A three-dimensional cohesive sediment transport model with data assimilation: Model development, sensitivity analysis and parameter estimation

Science.gov (United States)

Wang, Daosheng; Cao, Anzhou; Zhang, Jicai; Fan, Daidu; Liu, Yongzhi; Zhang, Yue

2018-06-01

Based on the theory of inverse problems, a three-dimensional sigma-coordinate cohesive sediment transport model with the adjoint data assimilation is developed. In this model, the physical processes of cohesive sediment transport, including deposition, erosion and advection-diffusion, are parameterized by corresponding model parameters. These parameters are usually poorly known and have traditionally been assigned empirically. By assimilating observations into the model, the model parameters can be estimated using the adjoint method; meanwhile, the data misfit between model results and observations can be decreased. The model developed in this work contains numerous parameters; therefore, it is necessary to investigate the parameter sensitivity of the model, which is assessed by calculating a relative sensitivity function and the gradient of the cost function with respect to each parameter. The results of parameter sensitivity analysis indicate that the model is sensitive to the initial conditions, inflow open boundary conditions, suspended sediment settling velocity and resuspension rate, while the model is insensitive to horizontal and vertical diffusivity coefficients. A detailed explanation of the pattern of sensitivity analysis is also given. In ideal twin experiments, constant parameters are estimated by assimilating 'pseudo' observations. The results show that the sensitive parameters are estimated more easily than the insensitive parameters. The conclusions of this work can provide guidance for the practical applications of this model to simulate sediment transport in the study area.

Uncertainty of Modal Parameters Estimated by ARMA Models

DEFF Research Database (Denmark)

Jensen, Jacob Laigaard; Brincker, Rune; Rytter, Anders

1990-01-01

In this paper the uncertainties of identified modal parameters such as eidenfrequencies and damping ratios are assed. From the measured response of dynamic excited structures the modal parameters may be identified and provide important structural knowledge. However the uncertainty of the parameters...... by simulation study of a lightly damped single degree of freedom system. Identification by ARMA models has been choosen as system identification method. It is concluded that both the sampling interval and number of sampled points may play a significant role with respect to the statistical errors. Furthermore......, it is shown that the model errors may also contribute significantly to the uncertainty....

Consistent Stochastic Modelling of Meteocean Design Parameters

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Sterndorff, M. J.

2000-01-01

Consistent stochastic models of metocean design parameters and their directional dependencies are essential for reliability assessment of offshore structures. In this paper a stochastic model for the annual maximum values of the significant wave height, and the associated wind velocity, current...

Can the "standard" unitarized Regge models describe the TOTEM data?

CERN Document Server

Alkin, A; Martynov, E

2013-01-01

The standard Regge poles are considered as inputs for two unitarization methods: eikonal and U-matrix. It is shown that only models with three input pomerons and two input odderons can describe the high energy data on $pp$ and $\\bar pp$ elastic scattering including the new data from Tevatron and LHC. However, it seems that the both considered models require a further modification (e.g. nonlinear reggeon trajectories and/or nonexponential vertex functions) for a more satisfactory description of the data at 19.0 GeV$\\leq \\sqrt{s}\\leq$ 7 TeV and 0.01 $\\leq |t|\\leq $14.2 GeV$^{2}$.

LINEAR MIXED MODEL TO DESCRIBE THE BASAL AREA INCREMENT FOR INDIVUDUAL CEDRO (Cedrela odorata L.TREES IN OCCIDENTAL AMAZON, BRAZIL

Directory of Open Access Journals (Sweden)

Thiago Augusto da Cunha

2013-01-01

Full Text Available Reliable growth data from trees are important to establish a rational forest management. Characteristics from trees, like the size, crown architecture and competition indices have been used to mathematically describe the increment efficiently when associated with them. However, the precise role of these effects in the growth-modeling destined to tropical trees needs to be further studied. Here it is reconstructed the basal area increment (BAI of individual Cedrela odorata trees, sampled at Amazon forest, to develop a growth- model using potential-predictors like: (1 classical tree size; (2 morphometric data; (3 competition and (4 social position including liana loads. Despite the large variation in tree size and growth, we observed that these kinds of predictor variables described well the BAI in level of individual tree. The fitted mixed model achieve a high efficiency (R2=92.7 % and predicted 3-years BAI over bark for trees of Cedrela odorata ranging from 10 to 110 cm at diameter at breast height. Tree height, steam slenderness and crown formal demonstrated high influence in the BAI growth model and explaining most of the growth variance (Partial R2=87.2%. Competition variables had negative influence on the BAI, however, explained about 7% of the total variation. The introduction of a random parameter on the regressions model (mixed modelprocedure has demonstrated a better significance approach to the data observed and showed more realistic predictions than the fixed model.

A statistical model describing combined irreversible electroporation and electroporation-induced blood-brain barrier disruption.

Science.gov (United States)

Sharabi, Shirley; Kos, Bor; Last, David; Guez, David; Daniels, Dianne; Harnof, Sagi; Mardor, Yael; Miklavcic, Damijan

2016-03-01

Electroporation-based therapies such as electrochemotherapy (ECT) and irreversible electroporation (IRE) are emerging as promising tools for treatment of tumors. When applied to the brain, electroporation can also induce transient blood-brain-barrier (BBB) disruption in volumes extending beyond IRE, thus enabling efficient drug penetration. The main objective of this study was to develop a statistical model predicting cell death and BBB disruption induced by electroporation. This model can be used for individual treatment planning. Cell death and BBB disruption models were developed based on the Peleg-Fermi model in combination with numerical models of the electric field. The model calculates the electric field thresholds for cell kill and BBB disruption and describes the dependence on the number of treatment pulses. The model was validated using in vivo experimental data consisting of rats brains MRIs post electroporation treatments. Linear regression analysis confirmed that the model described the IRE and BBB disruption volumes as a function of treatment pulses number (r(2) = 0.79; p disruption, the ratio increased with the number of pulses. BBB disruption radii were on average 67% ± 11% larger than IRE volumes. The statistical model can be used to describe the dependence of treatment-effects on the number of pulses independent of the experimental setup.

Modeling bioaccumulation in humans using poly-parameter linear free energy relationships (PPLFERS)

Energy Technology Data Exchange (ETDEWEB)

Undeman, Emma, E-mail: emma.undeman@itm.su.se; Czub, Gertje; McLachlan, Michael S.

2011-04-01

Chemical partition coefficients between environmental media and biological tissues are a key component of bioaccumulation models. The single-parameter linear free energy relationships (spLFERs) commonly used for predicting partitioning are often derived using apolar chemicals and may not accurately capture polar chemicals. In this study, a poly-parameter LFER (ppLFER) based model of organic chemical bioaccumulation in humans is presented. Chemical partitioning was described by an air-body partition coefficient that was a volume weighted average of ppLFER based partition coefficients for the major organs and tissues constituting the human body. This model was compared to a spLFER model treating the body as a mixture of lipid ({approx} octanol) and water. Although model agreement was good for hydrophobic chemicals (average difference 15% for log K{sub OW} > 4 and log K{sub OA} > 8), the ppLFER model predicted {approx} 90% lower body burdens for hydrophilic chemicals (log K{sub OW} < 0). This was mainly due to lower predictions of muscle and adipose tissue sorption capacity for these chemicals. A comparison of the predicted muscle and adipose tissue sorption capacities of hydrophilic chemicals with measurements indicated that the ppLFER and spLFER models' uncertainties were similar. Consequently, little benefit from the implementation of ppLFERs in this model was identified. - Research Highlights: {yields}Implementation of ppLFERs resulted in on average 90% lower predicted body burdens. {yields}Uncertainties in spLFER and ppLFER predictions were similar. {yields}The benefit from implementation of ppLFERs in bioaccumulation models was limited.

Development of computation model on the GoldSim platform for the radionuclide transport in the geosphere with the time-dependent parameters

International Nuclear Information System (INIS)

Koo, Shigeru; Inagaki, Manabu

2010-06-01

In the high-level radioactive waste (HLW) disposal system, numerical evaluation for radionuclide transport with the time-dependent parameters is necessary to evaluate various scenarios. In H12 report, numerical calculation code MESHNOTE and TIGER were used for the evaluation of some natural phenomena scenarios that had to handle the time-dependent parameters. In the future, the necessity of handling the time-dependent parameters will be expected to increase, and more efficient calculation and improvement of quality control of input/output parameters will be required. Therefore, for the purpose of corresponding this requirement, a radionuclide transport model has been developed on the GoldSim platform. The GoldSim is a general simulation software, that was used for the computation modeling of Yucca Mountain Project. The conceptual model, the mathematical model and the verification of the GoldSim model are described in this report. In the future, application resources on this report will be able to upgrade for perturbation scenarios analysis model and other conceptual models. (author)

Global parameter estimation for thermodynamic models of transcriptional regulation.

Science.gov (United States)

Suleimenov, Yerzhan; Ay, Ahmet; Samee, Md Abul Hassan; Dresch, Jacqueline M; Sinha, Saurabh; Arnosti, David N

2013-07-15

Deciphering the mechanisms involved in gene regulation holds the key to understanding the control of central biological processes, including human disease, population variation, and the evolution of morphological innovations. New experimental techniques including whole genome sequencing and transcriptome analysis have enabled comprehensive modeling approaches to study gene regulation. In many cases, it is useful to be able to assign biological significance to the inferred model parameters, but such interpretation should take into account features that affect these parameters, including model construction and sensitivity, the type of fitness calculation, and the effectiveness of parameter estimation. This last point is often neglected, as estimation methods are often selected for historical reasons or for computational ease. Here, we compare the performance of two parameter estimation techniques broadly representative of local and global approaches, namely, a quasi-Newton/Nelder-Mead simplex (QN/NMS) method and a covariance matrix adaptation-evolutionary strategy (CMA-ES) method. The estimation methods were applied to a set of thermodynamic models of gene transcription applied to regulatory elements active in the Drosophila embryo. Measuring overall fit, the global CMA-ES method performed significantly better than the local QN/NMS method on high quality data sets, but this difference was negligible on lower quality data sets with increased noise or on data sets simplified by stringent thresholding. Our results suggest that the choice of parameter estimation technique for evaluation of gene expression models depends both on quality of data, the nature of the models [again, remains to be established] and the aims of the modeling effort. Copyright © 2013 Elsevier Inc. All rights reserved.

A distributed approach for parameters estimation in System Biology models

International Nuclear Information System (INIS)

Mosca, E.; Merelli, I.; Alfieri, R.; Milanesi, L.

2009-01-01

Due to the lack of experimental measurements, biological variability and experimental errors, the value of many parameters of the systems biology mathematical models is yet unknown or uncertain. A possible computational solution is the parameter estimation, that is the identification of the parameter values that determine the best model fitting respect to experimental data. We have developed an environment to distribute each run of the parameter estimation algorithm on a different computational resource. The key feature of the implementation is a relational database that allows the user to swap the candidate solutions among the working nodes during the computations. The comparison of the distributed implementation with the parallel one showed that the presented approach enables a faster and better parameter estimation of systems biology models.

Assessment of Food Chain Pathway Parameters in Biosphere Models: Annual Progress Report for Fiscal Year 2004

Energy Technology Data Exchange (ETDEWEB)

Napier, Bruce A.; Krupka, Kenneth M.; Fellows, Robert J.; Cataldo, Dominic A.; Valenta, Michelle M.; Gilmore, Tyler J.

2004-12-02

This Annual Progress Report describes the work performed and summarizes some of the key observations to date on the U.S. Nuclear Regulatory Commission’s project Assessment of Food Chain Pathway Parameters in Biosphere Models, which was established to assess and evaluate a number of key parameters used in the food-chain models used in performance assessments of radioactive waste disposal facilities. Section 2 of this report describes activities undertaken to collect samples of soils from three regions of the United States, the Southeast, Northwest, and Southwest, and perform analyses to characterize their physical and chemical properties. Section 3 summarizes information gathered regarding agricultural practices and common and unusual crops grown in each of these three areas. Section 4 describes progress in studying radionuclide uptake in several representative crops from the three soil types in controlled laboratory conditions. Section 5 describes a range of international coordination activities undertaken by Project staff in order to support the underlying data needs of the Project. Section 6 provides a very brief summary of the status of the GENII Version 2 computer program, which is a “client” of the types of data being generated by the Project, and for which the Project will be providing training to the US NRC staff in the coming Fiscal Year. Several appendices provide additional supporting information.

Generating and Describing Affective Eye Behaviors

Science.gov (United States)

Mao, Xia; Li, Zheng

The manner of a person's eye movement conveys much about nonverbal information and emotional intent beyond speech. This paper describes work on expressing emotion through eye behaviors in virtual agents based on the parameters selected from the AU-Coded facial expression database and real-time eye movement data (pupil size, blink rate and saccade). A rule-based approach to generate primary (joyful, sad, angry, afraid, disgusted and surprise) and intermediate emotions (emotions that can be represented as the mixture of two primary emotions) utilized the MPEG4 FAPs (facial animation parameters) is introduced. Meanwhile, based on our research, a scripting tool, named EEMML (Emotional Eye Movement Markup Language) that enables authors to describe and generate emotional eye movement of virtual agents, is proposed.

Parameters Estimation of Geographically Weighted Ordinal Logistic Regression (GWOLR) Model

Science.gov (United States)

Zuhdi, Shaifudin; Retno Sari Saputro, Dewi; Widyaningsih, Purnami

2017-06-01

A regression model is the representation of relationship between independent variable and dependent variable. The dependent variable has categories used in the logistic regression model to calculate odds on. The logistic regression model for dependent variable has levels in the logistics regression model is ordinal. GWOLR model is an ordinal logistic regression model influenced the geographical location of the observation site. Parameters estimation in the model needed to determine the value of a population based on sample. The purpose of this research is to parameters estimation of GWOLR model using R software. Parameter estimation uses the data amount of dengue fever patients in Semarang City. Observation units used are 144 villages in Semarang City. The results of research get GWOLR model locally for each village and to know probability of number dengue fever patient categories.

Reflector modelization for neutronic diffusion and parameters identification

International Nuclear Information System (INIS)

Argaud, J.P.

1993-04-01

Physical parameters of neutronic diffusion equations can be adjusted to decrease calculations-measurements errors. The reflector being always difficult to modelize, we choose to elaborate a new reflector model and to use the parameters of this model as adjustment coefficients in the identification procedure. Using theoretical results, and also the physical behaviour of neutronic flux solutions, the reflector model consists then in its replacement by boundary conditions for the diffusion equations on the core only. This theoretical result of non-local operator relations leads then to some discrete approximations by taking into account the multiscaled behaviour, on the core-reflector interface, of neutronic diffusion solutions. The resulting model of this approach is then compared with previous reflector modelizations, and first results indicate that this new model gives the same representation of reflector for the core than previous. (author). 12 refs

Transient dynamic and modeling parameter sensitivity analysis of 1D solid oxide fuel cell model

International Nuclear Information System (INIS)

Huangfu, Yigeng; Gao, Fei; Abbas-Turki, Abdeljalil; Bouquain, David; Miraoui, Abdellatif

2013-01-01

Highlights: • A multiphysics, 1D, dynamic SOFC model is developed. • The presented model is validated experimentally in eight different operating conditions. • Electrochemical and thermal dynamic transient time expressions are given in explicit forms. • Parameter sensitivity is discussed for different semi-empirical parameters in the model. - Abstract: In this paper, a multiphysics solid oxide fuel cell (SOFC) dynamic model is developed by using a one dimensional (1D) modeling approach. The dynamic effects of double layer capacitance on the electrochemical domain and the dynamic effect of thermal capacity on thermal domain are thoroughly considered. The 1D approach allows the model to predict the non-uniform distributions of current density, gas pressure and temperature in SOFC during its operation. The developed model has been experimentally validated, under different conditions of temperature and gas pressure. Based on the proposed model, the explicit time constant expressions for different dynamic phenomena in SOFC have been given and discussed in detail. A parameters sensitivity study has also been performed and discussed by using statistical Multi Parameter Sensitivity Analysis (MPSA) method, in order to investigate the impact of parameters on the modeling accuracy

Coupled 1D-2D hydrodynamic inundation model for sewer overflow: Influence of modeling parameters

Directory of Open Access Journals (Sweden)

Adeniyi Ganiyu Adeogun

2015-10-01

Full Text Available This paper presents outcome of our investigation on the influence of modeling parameters on 1D-2D hydrodynamic inundation model for sewer overflow, developed through coupling of an existing 1D sewer network model (SWMM and 2D inundation model (BREZO. The 1D-2D hydrodynamic model was developed for the purpose of examining flood incidence due to surcharged water on overland surface. The investigation was carried out by performing sensitivity analysis on the developed model. For the sensitivity analysis, modeling parameters, such as mesh resolution Digital Elevation Model (DEM resolution and roughness were considered. The outcome of the study shows the model is sensitive to changes in these parameters. The performance of the model is significantly influenced, by the Manning's friction value, the DEM resolution and the area of the triangular mesh. Also, changes in the aforementioned modeling parameters influence the Flood characteristics, such as the inundation extent, the flow depth and the velocity across the model domain. Keywords: Inundation, DEM, Sensitivity analysis, Model coupling, Flooding

On the validity of evolutionary models with site-specific parameters.

Directory of Open Access Journals (Sweden)

Konrad Scheffler

Full Text Available Evolutionary models that make use of site-specific parameters have recently been criticized on the grounds that parameter estimates obtained under such models can be unreliable and lack theoretical guarantees of convergence. We present a simulation study providing empirical evidence that a simple version of the models in question does exhibit sensible convergence behavior and that additional taxa, despite not being independent of each other, lead to improved parameter estimates. Although it would be desirable to have theoretical guarantees of this, we argue that such guarantees would not be sufficient to justify the use of these models in practice. Instead, we emphasize the importance of taking the variance of parameter estimates into account rather than blindly trusting point estimates - this is standardly done by using the models to construct statistical hypothesis tests, which are then validated empirically via simulation studies.

Climate change decision-making: Model & parameter uncertainties explored

Energy Technology Data Exchange (ETDEWEB)

Dowlatabadi, H.; Kandlikar, M.; Linville, C.

1995-12-31

A critical aspect of climate change decision-making is uncertainties in current understanding of the socioeconomic, climatic and biogeochemical processes involved. Decision-making processes are much better informed if these uncertainties are characterized and their implications understood. Quantitative analysis of these uncertainties serve to inform decision makers about the likely outcome of policy initiatives, and help set priorities for research so that outcome ambiguities faced by the decision-makers are reduced. A family of integrated assessment models of climate change have been developed at Carnegie Mellon. These models are distinguished from other integrated assessment efforts in that they were designed from the outset to characterize and propagate parameter, model, value, and decision-rule uncertainties. The most recent of these models is ICAM 2.1. This model includes representation of the processes of demographics, economic activity, emissions, atmospheric chemistry, climate and sea level change and impacts from these changes and policies for emissions mitigation, and adaptation to change. The model has over 800 objects of which about one half are used to represent uncertainty. In this paper we show, that when considering parameter uncertainties, the relative contribution of climatic uncertainties are most important, followed by uncertainties in damage calculations, economic uncertainties and direct aerosol forcing uncertainties. When considering model structure uncertainties we find that the choice of policy is often dominated by model structure choice, rather than parameter uncertainties.

Error propagation of partial least squares for parameters optimization in NIR modeling

Science.gov (United States)

Du, Chenzhao; Dai, Shengyun; Qiao, Yanjiang; Wu, Zhisheng

2018-03-01

A novel methodology is proposed to determine the error propagation of partial least-square (PLS) for parameters optimization in near-infrared (NIR) modeling. The parameters include spectral pretreatment, latent variables and variable selection. In this paper, an open source dataset (corn) and a complicated dataset (Gardenia) were used to establish PLS models under different modeling parameters. And error propagation of modeling parameters for water quantity in corn and geniposide quantity in Gardenia were presented by both type І and type II error. For example, when variable importance in the projection (VIP), interval partial least square (iPLS) and backward interval partial least square (BiPLS) variable selection algorithms were used for geniposide in Gardenia, compared with synergy interval partial least squares (SiPLS), the error weight varied from 5% to 65%, 55% and 15%. The results demonstrated how and what extent the different modeling parameters affect error propagation of PLS for parameters optimization in NIR modeling. The larger the error weight, the worse the model. Finally, our trials finished a powerful process in developing robust PLS models for corn and Gardenia under the optimal modeling parameters. Furthermore, it could provide a significant guidance for the selection of modeling parameters of other multivariate calibration models.

Error propagation of partial least squares for parameters optimization in NIR modeling.

Science.gov (United States)

Du, Chenzhao; Dai, Shengyun; Qiao, Yanjiang; Wu, Zhisheng

2018-03-05

A novel methodology is proposed to determine the error propagation of partial least-square (PLS) for parameters optimization in near-infrared (NIR) modeling. The parameters include spectral pretreatment, latent variables and variable selection. In this paper, an open source dataset (corn) and a complicated dataset (Gardenia) were used to establish PLS models under different modeling parameters. And error propagation of modeling parameters for water quantity in corn and geniposide quantity in Gardenia were presented by both type І and type II error. For example, when variable importance in the projection (VIP), interval partial least square (iPLS) and backward interval partial least square (BiPLS) variable selection algorithms were used for geniposide in Gardenia, compared with synergy interval partial least squares (SiPLS), the error weight varied from 5% to 65%, 55% and 15%. The results demonstrated how and what extent the different modeling parameters affect error propagation of PLS for parameters optimization in NIR modeling. The larger the error weight, the worse the model. Finally, our trials finished a powerful process in developing robust PLS models for corn and Gardenia under the optimal modeling parameters. Furthermore, it could provide a significant guidance for the selection of modeling parameters of other multivariate calibration models. Copyright © 2017. Published by Elsevier B.V.

Evaluation of physiological parameters and their influence on doses calculated from two alternative dosimetric models for the gastrointestinal tract

International Nuclear Information System (INIS)

Lessard, E.T.; Skrable, K.W.

1981-01-01

Two dosimetric models, the catenary compartmental model (Be70) and the slug flow model (Sk75), are examined using three sets of physiological parameters: those proposed by Eve, those proposed by ICRP, and those obtained from the Textbook of Physiology and Biochemistry by Bell et al. The impact of physiological parameters on the dosimetry of the tract is illustrated by comparing calculated maximum permissible daily activity ingestion rates for single, unabsorbed, particle emitting radionuclides with an effective energy term of unity. The conclusions drawn from this intercomparison of six different cases are: Current dosimetric models which use physiological parameters described in this article do not significantly disagree, and for the determination of average dose equivalent rates to segments of the tract due to chronic, long term ingestion of any radionuclide, the catenary compartmental model is a mathematically simpler approach. The catenary model in addition has certain advantages for the calculation of the photon dose contribution to one segment from cumulated activity (disintegrations) in another segment

Parameter Optimisation for the Behaviour of Elastic Models over Time

DEFF Research Database (Denmark)

Mosegaard, Jesper

2004-01-01

Optimisation of parameters for elastic models is essential for comparison or finding equivalent behaviour of elastic models when parameters cannot simply be transferred or converted. This is the case with a large range of commonly used elastic models. In this paper we present a general method tha...

Determination of probability density functions for parameters in the Munson-Dawson model for creep behavior of salt

International Nuclear Information System (INIS)

Pfeifle, T.W.; Mellegard, K.D.; Munson, D.E.

1992-10-01

The modified Munson-Dawson (M-D) constitutive model that describes the creep behavior of salt will be used in performance assessment calculations to assess compliance of the Waste Isolation Pilot Plant (WIPP) facility with requirements governing the disposal of nuclear waste. One of these standards requires that the uncertainty of future states of the system, material model parameters, and data be addressed in the performance assessment models. This paper presents a method in which measurement uncertainty and the inherent variability of the material are characterized by treating the M-D model parameters as random variables. The random variables can be described by appropriate probability distribution functions which then can be used in Monte Carlo or structural reliability analyses. Estimates of three random variables in the M-D model were obtained by fitting a scalar form of the model to triaxial compression creep data generated from tests of WIPP salt. Candidate probability distribution functions for each of the variables were then fitted to the estimates and their relative goodness-of-fit tested using the Kolmogorov-Smirnov statistic. A sophisticated statistical software package obtained from BMDP Statistical Software, Inc. was used in the M-D model fitting. A separate software package, STATGRAPHICS, was used in fitting the candidate probability distribution functions to estimates of the variables. Skewed distributions, i.e., lognormal and Weibull, were found to be appropriate for the random variables analyzed

ADJUSTMENT OF MORPHOMETRIC PARAMETERS OF WATER BASINS BASED ON DIGITAL TERRAIN MODELS

Directory of Open Access Journals (Sweden)

Krasil'nikov Vitaliy Mikhaylovich

2012-10-01

Full Text Available The authors argue that effective use of water resources requires accurate morphometric characteristics of water basins. Accurate parameters are needed to analyze their condition, and to assure their appropriate control and operation. Today multiple water basins need their morphometric characteristics to be adjusted and properly stored. The procedure employed so far is based on plane geometric horizontals depicted onto topographic maps. It is described in the procedural guidelines issued in respect of the «Application of water resource regulations governing the operation of waterworks facilities of power plants». The technology described there is obsolete due to the availability of specialized software. The computer technique is based on a digital terrain model. The authors provide an overview of the technique implemented at Rybinsk and Gorkiy water basins in this article. Thus, the digital terrain model generated on the basis of the field data is used at Gorkiy water basin, while the model based on maps and charts is applied at Rybinsk water basin. The authors believe that the software technique can be applied to any other water basin on the basis of the analysis and comparison of morphometric characteristics of the two water basins.

Recommended direct simulation Monte Carlo collision model parameters for modeling ionized air transport processes

Energy Technology Data Exchange (ETDEWEB)

Swaminathan-Gopalan, Krishnan; Stephani, Kelly A., E-mail: ksteph@illinois.edu [Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801 (United States)

2016-02-15

A systematic approach for calibrating the direct simulation Monte Carlo (DSMC) collision model parameters to achieve consistency in the transport processes is presented. The DSMC collision cross section model parameters are calibrated for high temperature atmospheric conditions by matching the collision integrals from DSMC against ab initio based collision integrals that are currently employed in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and Data Parallel Line Relaxation (DPLR) high temperature computational fluid dynamics solvers. The DSMC parameter values are computed for the widely used Variable Hard Sphere (VHS) and the Variable Soft Sphere (VSS) models using the collision-specific pairing approach. The recommended best-fit VHS/VSS parameter values are provided over a temperature range of 1000-20 000 K for a thirteen-species ionized air mixture. Use of the VSS model is necessary to achieve consistency in transport processes of ionized gases. The agreement of the VSS model transport properties with the transport properties as determined by the ab initio collision integral fits was found to be within 6% in the entire temperature range, regardless of the composition of the mixture. The recommended model parameter values can be readily applied to any gas mixture involving binary collisional interactions between the chemical species presented for the specified temperature range.

On-line scheme for parameter estimation of nonlinear lithium ion battery equivalent circuit models using the simplified refined instrumental variable method for a modified Wiener continuous-time model

International Nuclear Information System (INIS)

Allafi, Walid; Uddin, Kotub; Zhang, Cheng; Mazuir Raja Ahsan Sha, Raja; Marco, James

2017-01-01

Highlights: •Off-line estimation approach for continuous-time domain for non-invertible function. •Model reformulated to multi-input-single-output; nonlinearity described by sigmoid. •Method directly estimates parameters of nonlinear ECM from the measured-data. •Iterative on-line technique leads to smoother convergence. •The model is validated off-line and on-line using NCA battery. -- Abstract: The accuracy of identifying the parameters of models describing lithium ion batteries (LIBs) in typical battery management system (BMS) applications is critical to the estimation of key states such as the state of charge (SoC) and state of health (SoH). In applications such as electric vehicles (EVs) where LIBs are subjected to highly demanding cycles of operation and varying environmental conditions leading to non-trivial interactions of ageing stress factors, this identification is more challenging. This paper proposes an algorithm that directly estimates the parameters of a nonlinear battery model from measured input and output data in the continuous time-domain. The simplified refined instrumental variable method is extended to estimate the parameters of a Wiener model where there is no requirement for the nonlinear function to be invertible. To account for nonlinear battery dynamics, in this paper, the typical linear equivalent circuit model (ECM) is enhanced by a block-oriented Wiener configuration where the nonlinear memoryless block following the typical ECM is defined to be a sigmoid static nonlinearity. The nonlinear Weiner model is reformulated in the form of a multi-input, single-output linear model. This linear form allows the parameters of the nonlinear model to be estimated using any linear estimator such as the well-established least squares (LS) algorithm. In this paper, the recursive least square (RLS) method is adopted for online parameter estimation. The approach was validated on experimental data measured from an 18650-type Graphite

Assessment of structural model and parameter uncertainty with a multi-model system for soil water balance models

Science.gov (United States)

Michalik, Thomas; Multsch, Sebastian; Frede, Hans-Georg; Breuer, Lutz

2016-04-01

Water for agriculture is strongly limited in arid and semi-arid regions and often of low quality in terms of salinity. The application of saline waters for irrigation increases the salt load in the rooting zone and has to be managed by leaching to maintain a healthy soil, i.e. to wash out salts by additional irrigation. Dynamic simulation models are helpful tools to calculate the root zone water fluxes and soil salinity content in order to investigate best management practices. However, there is little information on structural and parameter uncertainty for simulations regarding the water and salt balance of saline irrigation. Hence, we established a multi-model system with four different models (AquaCrop, RZWQM, SWAP, Hydrus1D/UNSATCHEM) to analyze the structural and parameter uncertainty by using the Global Likelihood and Uncertainty Estimation (GLUE) method. Hydrus1D/UNSATCHEM and SWAP were set up with multiple sets of different implemented functions (e.g. matric and osmotic stress for root water uptake) which results in a broad range of different model structures. The simulations were evaluated against soil water and salinity content observations. The posterior distribution of the GLUE analysis gives behavioral parameters sets and reveals uncertainty intervals for parameter uncertainty. Throughout all of the model sets, most parameters accounting for the soil water balance show a low uncertainty, only one or two out of five to six parameters in each model set displays a high uncertainty (e.g. pore-size distribution index in SWAP and Hydrus1D/UNSATCHEM). The differences between the models and model setups reveal the structural uncertainty. The highest structural uncertainty is observed for deep percolation fluxes between the model sets of Hydrus1D/UNSATCHEM (~200 mm) and RZWQM (~500 mm) that are more than twice as high for the latter. The model sets show a high variation in uncertainty intervals for deep percolation as well, with an interquartile range (IQR) of

Comparison of three nonlinear models to describe long-term tag shedding by lake trout

Science.gov (United States)

Fabrizio, Mary C.; Swanson, Bruce L.; Schram, Stephen T.; Hoff, Michael H.

1996-01-01

We estimated long-term tag-shedding rates for lake trout Salvelinus namaycush using two existing models and a model we developed to account for the observed permanence of some tags. Because tag design changed over the course of the study, we examined tag-shedding rates for three types of numbered anchor tags (Floy tags FD-67, FD-67C, and FD-68BC) and an unprinted anchor tag (FD-67F). Lake trout from the Gull Island Shoal region, Lake Superior, were double-tagged, and subsequent recaptures were monitored in annual surveys conducted from 1974 to 1992. We modeled tag-shedding rates, using time at liberty and probabilities of tag shedding estimated from fish released in 1974 and 1978–1983 and later recaptured. Long-term shedding of numbered anchor tags in lake trout was best described by a nonlinear model with two parameters: an instantaneous tag-shedding rate and a constant representing the proportion of tags that were never shed. Although our estimates of annual shedding rates varied with tag type (0.300 for FD-67, 0.441 for FD-67C, and 0.656 for FD-68BC), differences were not significant. About 36% of tags remained permanently affixed to the fish. Of the numbered tags that were shed (about 64%), two mechanisms contributed to tag loss: disintegration and dislodgment. Tags from about 11% of recaptured fish had disintegrated, but most tags were dislodged. Unprinted tags were shed at a significant but low rate immediately after release, but the long-term, annual shedding rate of these tags was only 0.013. Compared with unprinted tags, numbered tags dislodged at higher annual rates; we hypothesized that this was due to the greater frictional drag associated with the larger cross-sectional area of numbered tags.

Results from a new mathematical model of gastrointestinal transit that incorporates age and gender-dependent physiological parameters

International Nuclear Information System (INIS)

Stubbs, J.B.

1992-01-01

Recently published data on effects of age and gender-dependent GI physiology and motility have been used to develop a new mathematical model describing the transit and adsorption of substances through the GI tract. This mathematical description of GI tract kinetics utilises more physiologically accurate transit processes than the ICRP Report 30 GI model. The model uses a combination of zero and first-order kinetics to describe motility. Some of the physiological parameters that the new model uses are gender, age, phase of the menstrual cycle, meal composition and gastric phase (solid versus liquid). A computer algorithm based on this model has been derived and results for young males are compared to those of the ICRP 30 model. Comparison of gastrointestinal residence times for 99 Tc m and 111 In labelled compounds, as a function of gender and age, are also presented. (author)

Building damage assessment from PolSAR data using texture parameters of statistical model

Science.gov (United States)

Li, Linlin; Liu, Xiuguo; Chen, Qihao; Yang, Shuai

2018-04-01

Accurate building damage assessment is essential in providing decision support for disaster relief and reconstruction. Polarimetric synthetic aperture radar (PolSAR) has become one of the most effective means of building damage assessment, due to its all-day/all-weather ability and richer backscatter information of targets. However, intact buildings that are not parallel to the SAR flight pass (termed oriented buildings) and collapsed buildings share similar scattering mechanisms, both of which are dominated by volume scattering. This characteristic always leads to misjudgments between assessments of collapsed buildings and oriented buildings from PolSAR data. Because the collapsed buildings and the intact buildings (whether oriented or parallel buildings) have different textures, a novel building damage assessment method is proposed in this study to address this problem by introducing texture parameters of statistical models. First, the logarithms of the estimated texture parameters of different statistical models are taken as a new texture feature to describe the collapse of the buildings. Second, the collapsed buildings and intact buildings are distinguished using an appropriate threshold. Then, the building blocks are classified into three levels based on the building block collapse rate. Moreover, this paper also discusses the capability for performing damage assessment using texture parameters from different statistical models or using different estimators. The RADARSAT-2 and ALOS-1 PolSAR images are used to present and analyze the performance of the proposed method. The results show that using the texture parameters avoids the problem of confusing collapsed and oriented buildings and improves the assessment accuracy. The results assessed by using the K/G0 distribution texture parameters estimated based on the second moment obtain the highest extraction accuracies. For the RADARSAT-2 and ALOS-1 data, the overall accuracy (OA) for these three types of

Modelling tourists arrival using time varying parameter

Science.gov (United States)

Suciptawati, P.; Sukarsa, K. G.; Kencana, Eka N.

2017-06-01

The importance of tourism and its related sectors to support economic development and poverty reduction in many countries increase researchers’ attentions to study and model tourists’ arrival. This work is aimed to demonstrate time varying parameter (TVP) technique to model the arrival of Korean’s tourists to Bali. The number of Korean tourists whom visiting Bali for period January 2010 to December 2015 were used to model the number of Korean’s tourists to Bali (KOR) as dependent variable. The predictors are the exchange rate of Won to IDR (WON), the inflation rate in Korea (INFKR), and the inflation rate in Indonesia (INFID). Observing tourists visit to Bali tend to fluctuate by their nationality, then the model was built by applying TVP and its parameters were approximated using Kalman Filter algorithm. The results showed all of predictor variables (WON, INFKR, INFID) significantly affect KOR. For in-sample and out-of-sample forecast with ARIMA’s forecasted values for the predictors, TVP model gave mean absolute percentage error (MAPE) as much as 11.24 percent and 12.86 percent, respectively.

SPOTting model parameters using a ready-made Python package

Science.gov (United States)

Houska, Tobias; Kraft, Philipp; Breuer, Lutz

2015-04-01

The selection and parameterization of reliable process descriptions in ecological modelling is driven by several uncertainties. The procedure is highly dependent on various criteria, like the used algorithm, the likelihood function selected and the definition of the prior parameter distributions. A wide variety of tools have been developed in the past decades to optimize parameters. Some of the tools are closed source. Due to this, the choice for a specific parameter estimation method is sometimes more dependent on its availability than the performance. A toolbox with a large set of methods can support users in deciding about the most suitable method. Further, it enables to test and compare different methods. We developed the SPOT (Statistical Parameter Optimization Tool), an open source python package containing a comprehensive set of modules, to analyze and optimize parameters of (environmental) models. SPOT comes along with a selected set of algorithms for parameter optimization and uncertainty analyses (Monte Carlo, MC; Latin Hypercube Sampling, LHS; Maximum Likelihood, MLE; Markov Chain Monte Carlo, MCMC; Scuffled Complex Evolution, SCE-UA; Differential Evolution Markov Chain, DE-MCZ), together with several likelihood functions (Bias, (log-) Nash-Sutcliff model efficiency, Correlation Coefficient, Coefficient of Determination, Covariance, (Decomposed-, Relative-, Root-) Mean Squared Error, Mean Absolute Error, Agreement Index) and prior distributions (Binomial, Chi-Square, Dirichlet, Exponential, Laplace, (log-, multivariate-) Normal, Pareto, Poisson, Cauchy, Uniform, Weibull) to sample from. The model-independent structure makes it suitable to analyze a wide range of applications. We apply all algorithms of the SPOT package in three different case studies. Firstly, we investigate the response of the Rosenbrock function, where the MLE algorithm shows its strengths. Secondly, we study the Griewank function, which has a challenging response surface for

Model parameters estimation and sensitivity by genetic algorithms

International Nuclear Information System (INIS)

Marseguerra, Marzio; Zio, Enrico; Podofillini, Luca

2003-01-01

In this paper we illustrate the possibility of extracting qualitative information on the importance of the parameters of a model in the course of a Genetic Algorithms (GAs) optimization procedure for the estimation of such parameters. The Genetic Algorithms' search of the optimal solution is performed according to procedures that resemble those of natural selection and genetics: an initial population of alternative solutions evolves within the search space through the four fundamental operations of parent selection, crossover, replacement, and mutation. During the search, the algorithm examines a large amount of solution points which possibly carries relevant information on the underlying model characteristics. A possible utilization of this information amounts to create and update an archive with the set of best solutions found at each generation and then to analyze the evolution of the statistics of the archive along the successive generations. From this analysis one can retrieve information regarding the speed of convergence and stabilization of the different control (decision) variables of the optimization problem. In this work we analyze the evolution strategy followed by a GA in its search for the optimal solution with the aim of extracting information on the importance of the control (decision) variables of the optimization with respect to the sensitivity of the objective function. The study refers to a GA search for optimal estimates of the effective parameters in a lumped nuclear reactor model of literature. The supporting observation is that, as most optimization procedures do, the GA search evolves towards convergence in such a way to stabilize first the most important parameters of the model and later those which influence little the model outputs. In this sense, besides estimating efficiently the parameters values, the optimization approach also allows us to provide a qualitative ranking of their importance in contributing to the model output. The

WATGIS: A GIS-Based Lumped Parameter Water Quality Model

Science.gov (United States)

Glenn P. Fernandez; George M. Chescheir; R. Wayne Skaggs; Devendra M. Amatya

2002-01-01

A Geographic Information System (GIS)­based, lumped parameter water quality model was developed to estimate the spatial and temporal nitrogen­loading patterns for lower coastal plain watersheds in eastern North Carolina. The model uses a spatially distributed delivery ratio (DR) parameter to account for nitrogen retention or loss along a drainage network. Delivery...

ZZ-CENPL, Chinese Evaluated Nuclear Parameter Library. ZZ CENPL-DLS, Discrete Level Schemes and Gamma Branching Ratios Library; ZZ CENPL-FBP, Fission Barrier Parameter Library; ZZ CENPL-GDRP, Giant Dipole Resonance Parameter Library; ZZ CENPL-NLD, Nuclear Level Density Parameter Library; ZZ CENPL-MCC, Nuclear Ground State Atomic Masses Library; ZZ CENPL-OMP, Optical Model Parameter Library

International Nuclear Information System (INIS)

Su Zongdi

1995-01-01

Description of program or function: CENPL - GDRP (Giant Dipole Resonance Parameters for Gamma-Ray): - Format: special format described in documentation; - Nuclides: V, Mn, Co, Ni, Cu, Zn, Ga, Ge, As, Se, Rb, Sr, Y, Zr, Nb, Mo, Rh, Pd, Ag, Cd, In, Sn, Sb, Te, I, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Ho, Er, Lu, Ta, W, Re, Os, Ir, Pt, Au, Hg, Pb, Bi, Th, U, Np, Pu. - Origin: Experimental values offered by S.S. Dietrich and B.L. Berman. CENPL - FBP (Fission Barrier Parameter Sub-Library): - Format: special format described in documentation; - Nuclides: (1) 51 nuclei region from Th-230 to Cf-255, (2) 46 nuclei region from Th-229 to Cf-253, (3) 24 nuclei region from Pa-232 to Cf-253; - Origin: (1) Lynn, (2) Analysis of experimental data by Back et al., (3) Ohsawa. CENPL - DLS (Discrete level scheme and branch ratio of gamma decay: - Format: Special format described in documentation; - Origin: ENSDF - BNL. CENPL - NLD (Nuclear Level Density): - Format: Special format described in documentation; - Origin: Huang Zhongfu et al. CENPL - OMP (Optical model parameter sub-library): - Format: special format described in documentation ; - Origin: CENDL, ENDF/B-VI, JENDL-3. CENPL - MC (I) and (II) (Atomic masses and characteristic constants for nuclear ground states) : - Format: Brief table format; - Nuclides: 4760 nuclides ranging from Z=0 A=1 to Z=122 A=318. - Origin: Experimental data and systematic results evaluated by Wapstra, theoretical results calculated by Moller, ENSDF - BNL and Nuclear Wallet Cards. CENPL contains the following six sub-libraries: 1. Atomic Masses and Characteristic Constants for nuclear ground states (MCC). This data consists of calculated and in most cases also measured mass excesses, atomic masses, total binding energies, spins, parities, and half-lives of nuclear ground states, abundances, etc. for 4800 nuclides. 2. Discrete Level Schemes and branching ratios of gamma decay (DLS). The data on nuclear discrete levels are based on the Evaluated

Practical expressions describing detective quantum efficiency in flat-panel detectors

Science.gov (United States)

Kim, H. K.

2011-11-01

In radiology, image quality excellence is a balance between system performance and patient dose, hence x-ray systems must be designed to ensure the maximum image quality is obtained for the lowest consistent dose. The concept of detective quantum efficiency (DQE) is widely used to quantify, understand, measure, and predict the performance of x-ray detectors and imaging systems. Cascaded linear-systems theory can be used to estimate DQE based on the system design parameters and this theoretical DQE can be utilized for determining the impact of various physical processes, such as secondary quantum sinks, noise aliasing, reabsorption noise, and others. However, the prediction of DQE usually requires tremendous efforts to determine each parameter consisting of the cascaded linear-systems model. In this paper, practical DQE formalisms assessing both the photoconductor- and scintillator-based flat-panel detectors under quantum-noise-limited operation are described. The developed formalisms are experimentally validated and discussed for their limits. The formalisms described in this paper would be helpful for the rapid prediction of the DQE performances of developing systems as well as the optimal design of systems.

Standard model parameters determination and validity tests in Z{sup 0} hadronic disintegrations; Determination des parametres du modele standard et tests de sa validite dans les desintegrations hadroniques du Z{sup 0}

Energy Technology Data Exchange (ETDEWEB)

Todorov, T

1993-05-01

This thesis describes the determination of the electroweak parameters from the measurements of the total hadronic cross-section by the DELPHI experiment at LEP-I. The analysed data was taken in the years 1991 and 1992; a previous analysis of the data taken in 1990 is included in the final fits. The first part of the thesis describes the interest of the measurement of the Z{sup 0} resonance parameters in the framework of the Standard Model as well as their implications for alternative models. The Standard Model predictions are described in some detail, and their precision is estimated. Then follows a brief description of the LEP collider, of the measurement of the collision energy, and of the experimental setup. A chapter is devoted to the description of the luminosity measurement, essential for the determination of total cross-sections. The measurement of the hadronic cross-section (event selection, study of backgrounds, study of sources of systematic uncertainties) is described in detail in the next chapter. Then follows a description of the method of the extraction of the resonance parameters, and a discussion of the uncertainties in their determination. The values obtained are interpreted in the framework of the Standard Model, as well as in the framework of some more general theories. Finally, the event generator for hadron production in two-photon collisions is described in the appendix. (author). 69 refs., 51 figs., 9 tabs., 1 ann.

Consensus conference on core radiological parameters to describe lumbar stenosis - an initiative for structured reporting

Energy Technology Data Exchange (ETDEWEB)

Andreisek, Gustav; Winklhofer, Sebastian F.X. [University Hospital Zurich, Department of Radiology, Zurich (Switzerland); Deyo, Richard A. [Oregon Health and Science University, Portland, OR (United States); Jarvik, Jeffrey G. [University of Washington, Seattle, WA (United States); Porchet, Francois [Schulthess Klinik, Zuerich (Switzerland); Steurer, Johann [University Hospital Zurich, Horten Center for patient oriented research and knowledge transfer, Zurich (Switzerland); Collaboration: On behalf of the LSOS working group

2014-12-15

To define radiological criteria and parameters as a minimum standard in a structured radiological report for patients with lumbar spinal stenosis (LSS) and to identify criteria and parameters for research purposes. All available radiological criteria and parameters for LSS were identified using systematic literature reviews and a Delphi survey. We invited to the consensus meeting, and provided data, to 15 internationally renowned experts from different countries. During the meeting, these experts reached consensus in a structured and systematic discussion about a core list of radiological criteria and parameters for standard reporting. We identified a total of 27 radiological criteria and parameters for LSS. During the meeting, the experts identified five of these as core items for a structured report. For central stenosis, these were ''compromise of the central zone'' and ''relation between fluid and cauda equina''. For lateral stenosis, the group agreed that ''nerve root compression in the lateral recess'' was a core item. For foraminal stenosis, we included ''nerve root impingement'' and ''compromise of the foraminal zone''. As a minimum standard, five radiological criteria should be used in a structured radiological report in LSS. Other parameters are well suited for research. (orig.)

NONLINEAR PLANT PIECEWISE-CONTINUOUS MODEL MATRIX PARAMETERS ESTIMATION

Directory of Open Access Journals (Sweden)

Roman L. Leibov

2017-09-01

Full Text Available This paper presents a nonlinear plant piecewise-continuous model matrix parameters estimation technique using nonlinear model time responses and random search method. One of piecewise-continuous model application areas is defined. The results of proposed approach application for aircraft turbofan engine piecewisecontinuous model formation are presented