Optimizing Methods of Obtaining Stellar Parameters for the H3 Survey

Science.gov (United States)

Ivory, KeShawn; Conroy, Charlie; Cargile, Phillip

2018-01-01

The Stellar Halo at High Resolution with Hectochelle Survey (H3) is in the process of observing and collecting stellar parameters for stars in the Milky Way's halo. With a goal of measuring radial velocities for fainter stars, it is crucial that we have optimal methods of obtaining this and other parameters from the data from these stars.The method currently developed is The Payne, named after Cecilia Payne-Gaposchkin, a code that uses neural networks and Markov Chain Monte Carlo methods to utilize both spectra and photometry to obtain values for stellar parameters. This project was to investigate the benefit of fitting both spectra and spectral energy distributions (SED). Mock spectra using the parameters of the Sun were created and noise was inserted at various signal to noise values. The Payne then fit each mock spectrum with and without a mock SED also generated from solar parameters. The result was that at high signal to noise, the spectrum dominated and the effect of fitting the SED was minimal. But at low signal to noise, the addition of the SED greatly decreased the standard deviation of the data and resulted in more accurate values for temperature and metallicity.

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.

Dynamics of a neuron model in different two-dimensional parameter-spaces

International Nuclear Information System (INIS)

Rech, Paulo C.

2011-01-01

We report some two-dimensional parameter-space diagrams numerically obtained for the multi-parameter Hindmarsh-Rose neuron model. Several different parameter planes are considered, and we show that regardless of the combination of parameters, a typical scenario is preserved: for all choice of two parameters, the parameter-space presents a comb-shaped chaotic region immersed in a large periodic region. We also show that exist regions close these chaotic region, separated by the comb teeth, organized themselves in period-adding bifurcation cascades. - Research highlights: → We report parameter-spaces obtained for the Hindmarsh-Rose neuron model. → Regardless of the combination of parameters, a typical scenario is preserved. → The scenario presents a comb-shaped chaotic region immersed in a periodic region. → Periodic regions near the chaotic region are in period-adding bifurcation cascades.

Four-parameter analytical local model potential for atoms

International Nuclear Information System (INIS)

Fei, Yu; Jiu-Xun, Sun; Rong-Gang, Tian; Wei, Yang

2009-01-01

Analytical local model potential for modeling the interaction in an atom reduces the computational effort in electronic structure calculations significantly. A new four-parameter analytical local model potential is proposed for atoms Li through Lr, and the values of four parameters are shell-independent and obtained by fitting the results of X a method. At the same time, the energy eigenvalues, the radial wave functions and the total energies of electrons are obtained by solving the radial Schrödinger equation with a new form of potential function by Numerov's numerical method. The results show that our new form of potential function is suitable for high, medium and low Z atoms. A comparison among the new potential function and other analytical potential functions shows the greater flexibility and greater accuracy of the present new potential function. (atomic and molecular physics)

Optimization-Based Inverse Identification of the Parameters of a Concrete Cap Material Model

Science.gov (United States)

Král, Petr; Hokeš, Filip; Hušek, Martin; Kala, Jiří; Hradil, Petr

2017-10-01

Issues concerning the advanced numerical analysis of concrete building structures in sophisticated computing systems currently require the involvement of nonlinear mechanics tools. The efforts to design safer, more durable and mainly more economically efficient concrete structures are supported via the use of advanced nonlinear concrete material models and the geometrically nonlinear approach. The application of nonlinear mechanics tools undoubtedly presents another step towards the approximation of the real behaviour of concrete building structures within the framework of computer numerical simulations. However, the success rate of this application depends on having a perfect understanding of the behaviour of the concrete material models used and having a perfect understanding of the used material model parameters meaning. The effective application of nonlinear concrete material models within computer simulations often becomes very problematic because these material models very often contain parameters (material constants) whose values are difficult to obtain. However, getting of the correct values of material parameters is very important to ensure proper function of a concrete material model used. Today, one possibility, which permits successful solution of the mentioned problem, is the use of optimization algorithms for the purpose of the optimization-based inverse material parameter identification. Parameter identification goes hand in hand with experimental investigation while it trying to find parameter values of the used material model so that the resulting data obtained from the computer simulation will best approximate the experimental data. This paper is focused on the optimization-based inverse identification of the parameters of a concrete cap material model which is known under the name the Continuous Surface Cap Model. Within this paper, material parameters of the model are identified on the basis of interaction between nonlinear computer simulations

Personalization of models with many model parameters: an efficient sensitivity analysis approach.

Science.gov (United States)

Donders, W P; Huberts, W; van de Vosse, F N; Delhaas, T

2015-10-01

Uncertainty quantification and global sensitivity analysis are indispensable for patient-specific applications of models that enhance diagnosis or aid decision-making. Variance-based sensitivity analysis methods, which apportion each fraction of the output uncertainty (variance) to the effects of individual input parameters or their interactions, are considered the gold standard. The variance portions are called the Sobol sensitivity indices and can be estimated by a Monte Carlo (MC) approach (e.g., Saltelli's method [1]) or by employing a metamodel (e.g., the (generalized) polynomial chaos expansion (gPCE) [2, 3]). All these methods require a large number of model evaluations when estimating the Sobol sensitivity indices for models with many parameters [4]. To reduce the computational cost, we introduce a two-step approach. In the first step, a subset of important parameters is identified for each output of interest using the screening method of Morris [5]. In the second step, a quantitative variance-based sensitivity analysis is performed using gPCE. Efficient sampling strategies are introduced to minimize the number of model runs required to obtain the sensitivity indices for models considering multiple outputs. The approach is tested using a model that was developed for predicting post-operative flows after creation of a vascular access for renal failure patients. We compare the sensitivity indices obtained with the novel two-step approach with those obtained from a reference analysis that applies Saltelli's MC method. The two-step approach was found to yield accurate estimates of the sensitivity indices at two orders of magnitude lower computational cost. Copyright © 2015 John Wiley & Sons, Ltd.

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.

Electro-optical parameters of bond polarizability model for aluminosilicates.

Science.gov (United States)

Smirnov, Konstantin S; Bougeard, Daniel; Tandon, Poonam

2006-04-06

Electro-optical parameters (EOPs) of bond polarizability model (BPM) for aluminosilicate structures were derived from quantum-chemical DFT calculations of molecular models. The tensor of molecular polarizability and the derivatives of the tensor with respect to the bond length are well reproduced with the BPM, and the EOPs obtained are in a fair agreement with available experimental data. The parameters derived were found to be transferable to larger molecules. This finding suggests that the procedure used can be applied to systems with partially ionic chemical bonds. The transferability of the parameters to periodic systems was tested in molecular dynamics simulation of the polarized Raman spectra of alpha-quartz. It appeared that the molecular Si-O bond EOPs failed to reproduce the intensity of peaks in the spectra. This limitation is due to large values of the longitudinal components of the bond polarizability and its derivative found in the molecular calculations as compared to those obtained from periodic DFT calculations of crystalline silica polymorphs by Umari et al. (Phys. Rev. B 2001, 63, 094305). It is supposed that the electric field of the solid is responsible for the difference of the parameters. Nevertheless, the EOPs obtained can be used as an initial set of parameters for calculations of polarizability related characteristics of relevant systems in the framework of BPM.

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.

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)

ENDF-6 File 30: Data covariances obtained from parameter covariances and sensitivities

International Nuclear Information System (INIS)

Muir, D.W.

1989-01-01

File 30 is provided as a means of describing the covariances of tabulated cross sections, multiplicities, and energy-angle distributions that result from propagating the covariances of a set of underlying parameters (for example, the input parameters of a nuclear-model code), using an evaluator-supplied set of parameter covariances and sensitivities. Whenever nuclear data are evaluated primarily through the application of nuclear models, the covariances of the resulting data can be described very adequately, and compactly, by specifying the covariance matrix for the underlying nuclear parameters, along with a set of sensitivity coefficients giving the rate of change of each nuclear datum of interest with respect to each of the model parameters. Although motivated primarily by these applications of nuclear theory, use of File 30 is not restricted to any one particular evaluation methodology. It can be used to describe data covariances of any origin, so long as they can be formally separated into a set of parameters with specified covariances and a set of data sensitivities

Rain storm models and the relationship between their parameters

NARCIS (Netherlands)

Stol, P.T.

1977-01-01

Rainfall interstation correlation functions can be obtained with the aid of analytic rainfall or storm models. Since alternative storm models have different mathematical formulas, comparison should be based on equallity of parameters like storm diameter, mean rainfall amount, storm maximum or total

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...

One parameter model potential for noble metals

International Nuclear Information System (INIS)

Idrees, M.; Khwaja, F.A.; Razmi, M.S.K.

1981-08-01

A phenomenological one parameter model potential which includes s-d hybridization and core-core exchange contributions is proposed for noble metals. A number of interesting properties like liquid metal resistivities, band gaps, thermoelectric powers and ion-ion interaction potentials are calculated for Cu, Ag and Au. The results obtained are in better agreement with experiment than the ones predicted by the other model potentials in the literature. (author)

Geometry parameters for musculoskeletal modelling of the shoulder system

NARCIS (Netherlands)

Van der Helm, F C; Veeger, DirkJan (H. E. J.); Pronk, G M; Van der Woude, L H; Rozendal, R H

A dynamical finite-element model of the shoulder mechanism consisting of thorax, clavicula, scapula and humerus is outlined. The parameters needed for the model are obtained in a cadaver experiment consisting of both shoulders of seven cadavers. In this paper, in particular, the derivation of

Dynamics in the Parameter Space of a Neuron Model

Science.gov (United States)

Paulo, C. Rech

2012-06-01

Some two-dimensional parameter-space diagrams are numerically obtained by considering the largest Lyapunov exponent for a four-dimensional thirteen-parameter Hindmarsh—Rose neuron model. Several different parameter planes are considered, and it is shown that depending on the combination of parameters, a typical scenario can be preserved: for some choice of two parameters, the parameter plane presents a comb-shaped chaotic region embedded in a large periodic region. It is also shown that there exist regions close to these comb-shaped chaotic regions, separated by the comb teeth, organizing themselves in period-adding bifurcation cascades.

Model-based verification method for solving the parameter uncertainty in the train control system

International Nuclear Information System (INIS)

Cheng, Ruijun; Zhou, Jin; Chen, Dewang; Song, Yongduan

2016-01-01

This paper presents a parameter analysis method to solve the parameter uncertainty problem for hybrid system and explore the correlation of key parameters for distributed control system. For improving the reusability of control model, the proposed approach provides the support for obtaining the constraint sets of all uncertain parameters in the abstract linear hybrid automata (LHA) model when satisfying the safety requirements of the train control system. Then, in order to solve the state space explosion problem, the online verification method is proposed to monitor the operating status of high-speed trains online because of the real-time property of the train control system. Furthermore, we construct the LHA formal models of train tracking model and movement authority (MA) generation process as cases to illustrate the effectiveness and efficiency of the proposed method. In the first case, we obtain the constraint sets of uncertain parameters to avoid collision between trains. In the second case, the correlation of position report cycle and MA generation cycle is analyzed under both the normal and the abnormal condition influenced by packet-loss factor. Finally, considering stochastic characterization of time distributions and real-time feature of moving block control system, the transient probabilities of wireless communication process are obtained by stochastic time petri nets. - Highlights: • We solve the parameters uncertainty problem by using model-based method. • We acquire the parameter constraint sets by verifying linear hybrid automata models. • Online verification algorithms are designed to monitor the high-speed trains. • We analyze the correlation of key parameters and uncritical parameters. • The transient probabilities are obtained by using reliability analysis.

A New Five-Parameter Fréchet Model for Extreme Values

Directory of Open Access Journals (Sweden)

Muhammad Ahsan ul Haq

2017-09-01

Full Text Available A new five parameter Fréchet model for Extreme Values was proposed and studied. Various mathematical properties including moments, quantiles, and moment generating function were derived. Incomplete moments and probability weighted moments were also obtained. The maximum likelihood method was used to estimate the model parameters. The flexibility of the derived model was accessed using two real data set applications.

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.

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...

Determination of the Corona model parameters with artificial neural networks

International Nuclear Information System (INIS)

Ahmet, Nayir; Bekir, Karlik; Arif, Hashimov

2005-01-01

Full text : The aim of this study is to calculate new model parameters taking into account the corona of electrical transmission line wires. For this purpose, a neural network modeling proposed for the corona frequent characteristics modeling. Then this model was compared with the other model developed at the Polytechnic Institute of Saint Petersburg. The results of development of the specified corona model for calculation of its influence on the wave processes in multi-wires line and determination of its parameters are submitted. Results of obtained calculation equations are brought for electrical transmission line with allowance for superficial effect in the ground and wires with reference to developed corona model

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

photosynthesis irradiance functions and parameters for modeling in situ production profiles. In light of the results obtained in this work we argue that the choice of the primary production model should reflect the available data and these models should be data driven regarding parameter estimation.

Determining extreme parameter correlation in ground water models

DEFF Research Database (Denmark)

Hill, Mary Cole; Østerby, Ole

2003-01-01

can go undetected even by experienced modelers. Extreme parameter correlation can be detected using parameter correlation coefficients, but their utility depends on the presence of sufficient, but not excessive, numerical imprecision of the sensitivities, such as round-off error. This work...... investigates the information that can be obtained from parameter correlation coefficients in the presence of different levels of numerical imprecision, and compares it to the information provided by an alternative method called the singular value decomposition (SVD). Results suggest that (1) calculated...... correlation coefficients with absolute values that round to 1.00 were good indicators of extreme parameter correlation, but smaller values were not necessarily good indicators of lack of correlation and resulting unique parameter estimates; (2) the SVD may be more difficult to interpret than parameter...

Regionalising Parameters of a Conceptual Rainfall-Runoff Model for ...

African Journals Online (AJOL)

IHACRES, a lumped conceptual rainfall-runoff model, was calibrated to six catchments ranging in size from 49km2 to 600 km2 within the upper Tana River basin to obtain a set of model parameters that characterise the hydrological behaviour within the region. Physical catchment attributes indexing topography, soil and ...

Condition Parameter Modeling for Anomaly Detection in Wind Turbines

Directory of Open Access Journals (Sweden)

Yonglong Yan

2014-05-01

Full Text Available Data collected from the supervisory control and data acquisition (SCADA system, used widely in wind farms to obtain operational and condition information about wind turbines (WTs, is of important significance for anomaly detection in wind turbines. The paper presents a novel model for wind turbine anomaly detection mainly based on SCADA data and a back-propagation neural network (BPNN for automatic selection of the condition parameters. The SCADA data sets are determined through analysis of the cumulative probability distribution of wind speed and the relationship between output power and wind speed. The automatic BPNN-based parameter selection is for reduction of redundant parameters for anomaly detection in wind turbines. Through investigation of cases of WT faults, the validity of the automatic parameter selection-based model for WT anomaly detection is verified.

Possibilities of paired comparison of receptor binding parameters obtained in a single experiment

International Nuclear Information System (INIS)

Mashilova, K.V.; Kiriakov, G.V.; Malin, K.M.; Rozhanets, V.V.

1987-01-01

The authors explore the use of comparing control and experimental groups on the basis of extrapolation parameters obtained from a single pair of experiments. One of the experiments study the parameters of specific binding of 3 H-D-ala-2-enkephelin-5-D-leucine in the striatum of different groups of rats. The analysis was done by the Cornish-Bowden method

Dynamics of a neuron model in different two-dimensional parameter-spaces

Science.gov (United States)

Rech, Paulo C.

2011-03-01

We report some two-dimensional parameter-space diagrams numerically obtained for the multi-parameter Hindmarsh-Rose neuron model. Several different parameter planes are considered, and we show that regardless of the combination of parameters, a typical scenario is preserved: for all choice of two parameters, the parameter-space presents a comb-shaped chaotic region immersed in a large periodic region. We also show that exist regions close these chaotic region, separated by the comb teeth, organized themselves in period-adding bifurcation cascades.

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.

Determination of milling parameters to obtain mechanosynthesized ZnFe2O4

International Nuclear Information System (INIS)

Jean, Malick; Nachbaur, Virginie

2008-01-01

In this work, the mechanosynthesis of zinc ferrite in WC vials is studied. Millings are performed under air, with a planetary ball-milling, starting from elemental oxides α-Fe 2 O 3 and ZnO. As-milled powders are structurally and magnetically characterized by X-ray diffraction and Moessbauer spectroscopy. Milling parameters as rotational speeds of main disc and vials are particularly discussed in terms of influence on the obtaining of a pure zinc ferrite phase. These parameters have a strong influence on injected power, on radial and tangential components of the impact force. Friction phenomenon, associated with injected power, have been found to be the governing parameters of the end product

An automatic scaling method for obtaining the trace and parameters from oblique ionogram based on hybrid genetic algorithm

Science.gov (United States)

Song, Huan; Hu, Yaogai; Jiang, Chunhua; Zhou, Chen; Zhao, Zhengyu; Zou, Xianjian

2016-12-01

Scaling oblique ionogram plays an important role in obtaining ionospheric structure at the midpoint of oblique sounding path. The paper proposed an automatic scaling method to extract the trace and parameters of oblique ionogram based on hybrid genetic algorithm (HGA). The extracted 10 parameters come from F2 layer and Es layer, such as maximum observation frequency, critical frequency, and virtual height. The method adopts quasi-parabolic (QP) model to describe F2 layer's electron density profile that is used to synthesize trace. And it utilizes secant theorem, Martyn's equivalent path theorem, image processing technology, and echoes' characteristics to determine seven parameters' best fit values, and three parameter's initial values in QP model to set up their searching spaces which are the needed input data of HGA. Then HGA searches the three parameters' best fit values from their searching spaces based on the fitness between the synthesized trace and the real trace. In order to verify the performance of the method, 240 oblique ionograms are scaled and their results are compared with manual scaling results and the inversion results of the corresponding vertical ionograms. The comparison results show that the scaling results are accurate or at least adequate 60-90% of the time.

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

Analytical calculation of detailed model parameters of cast resin dry-type transformers

International Nuclear Information System (INIS)

Eslamian, M.; Vahidi, B.; Hosseinian, S.H.

2011-01-01

Highlights: → In this paper high frequency behavior of cast resin dry-type transformers was simulated. → Parameters of detailed model were calculated using analytical method and compared with FEM results. → A lab transformer was constructed in order to compare theoretical and experimental results. -- Abstract: Non-flammable characteristic of cast resin dry-type transformers make them suitable for different kind of usages. This paper presents an analytical method of how to obtain parameters of detailed model of these transformers. The calculated parameters are compared and verified with the corresponding FEM results and if it was necessary, correction factors are introduced for modification of the analytical solutions. Transient voltages under full and chopped test impulses are calculated using the obtained detailed model. In order to validate the model, a setup was constructed for testing on high-voltage winding of cast resin dry-type transformer. The simulation results were compared with the experimental data measured from FRA and impulse tests.

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.

Interpretation of acoustic parameters obtained by EMAR measurement for non-destructive hydrogen concentration measurement in Zr alloy

International Nuclear Information System (INIS)

Nakatsuka, Masafumi; Uchida, Katsuya; Miyazaki, Akihiro; Ishii, Yoshiaki

2007-01-01

An obvious quantitative relation between hydrogen concentrations in zirconium alloy and acoustic anisotropy parameters obtained by the electromagnetic acoustic resonance (EMAR) method was reported. To elucidate the mechanism, the acoustic parameters were calculated based on the elastic theory and the equation of motion. The acoustic parameters of obtained by the EMAR method were interpreted quantitatively using the anisotropic elastic constants of the specimen, and value calculated from texture data for non-hydrogen charged specimens showed good agreement with those obtained by the EMAR method. Calculated temperature dependence of the acoustic anisotropy for the non-hydrogen charged specimen also agreed well with that by the EMAR method. The consistencies demonstrated that the absolute values of the acoustic parameters for non-hydrogen charged specimen can be calculated from both the texture data of (0002) pole figure and the elastic constants of the specimen. Hydrogen addition up to approximately 650ppm was found not to change the original (0002) pole figure and, correspondingly, no hydrogen concentration dependence of the acoustic parameters was obtained from the calculation. These results implied that the zirconium hydride itself played an important role for the change in the acoustic parameters of the hydrogen charged specimens, and the importance of obtaining the information on the elastic constants of the zirconium hydride was pointed out. (author)

Parameters Calculation of ZnO Surge Arrester Models by Genetic Algorithms

Directory of Open Access Journals (Sweden)

A. Bayadi

2006-09-01

Full Text Available This paper proposes to provide a new technique based on the genetic algorithm to obtain the best possible series of values of the parameters of the ZnO surge arresters models. The validity of the predicted parameters is then checked by comparing the results predicted with the experimental results available in the literature. Using the ATP-EMTP package an application of the arrester model on network system studies is presented and discussed.

Repetitive Identification of Structural Systems Using a Nonlinear Model Parameter Refinement Approach

Directory of Open Access Journals (Sweden)

Jeng-Wen Lin

2009-01-01

Full Text Available This paper proposes a statistical confidence interval based nonlinear model parameter refinement approach for the health monitoring of structural systems subjected to seismic excitations. The developed model refinement approach uses the 95% confidence interval of the estimated structural parameters to determine their statistical significance in a least-squares regression setting. When the parameters' confidence interval covers the zero value, it is statistically sustainable to truncate such parameters. The remaining parameters will repetitively undergo such parameter sifting process for model refinement until all the parameters' statistical significance cannot be further improved. This newly developed model refinement approach is implemented for the series models of multivariable polynomial expansions: the linear, the Taylor series, and the power series model, leading to a more accurate identification as well as a more controllable design for system vibration control. Because the statistical regression based model refinement approach is intrinsically used to process a “batch” of data and obtain an ensemble average estimation such as the structural stiffness, the Kalman filter and one of its extended versions is introduced to the refined power series model for structural health monitoring.

A new approach to the extraction of single exponential diode model parameters

Science.gov (United States)

Ortiz-Conde, Adelmo; García-Sánchez, Francisco J.

2018-06-01

A new integration method is presented for the extraction of the parameters of a single exponential diode model with series resistance from the measured forward I-V characteristics. The extraction is performed using auxiliary functions based on the integration of the data which allow to isolate the effects of each of the model parameters. A differentiation method is also presented for data with low level of experimental noise. Measured and simulated data are used to verify the applicability of both proposed method. Physical insight about the validity of the model is also obtained by using the proposed graphical determinations of the parameters.

Parameter extraction of different fuel cell models with transferred adaptive differential evolution

International Nuclear Information System (INIS)

Gong, Wenyin; Yan, Xuesong; Liu, Xiaobo; Cai, Zhihua

2015-01-01

To improve the design and control of FC (fuel cell) models, it is important to extract their unknown parameters. Generally, the parameter extraction problems of FC models can be transformed as nonlinear and multi-variable optimization problems. To extract the parameters of different FC models exactly and fast, in this paper, we propose a transferred adaptive DE (differential evolution) framework, in which the successful parameters of the adaptive DE solving previous problems are properly transferred to solve new optimization problems in the similar problem-domains. Based on this framework, an improved adaptive DE method (TRADE, in short) is presented as an illustration. To verify the performance of our proposal, TRADE is used to extract the unknown parameters of two types of fuel cell models, i.e., PEMFC (proton exchange membrane fuel cell) and SOFC (solid oxide fuel cell). The results of TRADE are also compared with those of other state-of-the-art EAs (evolutionary algorithms). Even though the modification is very simple, the results indicate that TRADE can extract the parameters of both PEMFC and SOFC models exactly and fast. Moreover, the V–I characteristics obtained by TRADE agree well with the simulated and experimental data in all cases for both types of fuel cell models. Also, it improves the performance of the original adaptive DE significantly in terms of both the quality of final solutions and the convergence speed in all cases. Additionally, TRADE is able to provide better results compared with other EAs. - Highlights: • A framework of transferred adaptive differential evolution is proposed. • Based on the framework, an improved differential evolution (TRADE) is presented. • TRADE obtains very promising results to extract the parameters of PEMFC and SOFC models

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.)

Evaluation of ACE inhibitors lipophilicity using in silico and chromatographically obtained hydrophobicity parameters

Directory of Open Access Journals (Sweden)

Odović Jadranka V.

2013-01-01

Full Text Available The aim of this study was to compare different calculation methods to determine lipophilicity, expressed as logP value, of seven ACE inhibitors (enalapril, quinapril, fosinopril, lisinopril, cilazapril, ramipril, and benazapril with significantly different structure. Experimentally determined n-octanol/water partition coefficients, logPO/W values, were obtained from relevant literature. The correlations between all collected logP values were studied and the best agreements between calculated logP and experimentally determined logPO/W values, were observed for KOWWINlogP or MilogP values (r = 0.999 or r = 0.974, respectively. The correlations between all collected logP values and chromatographically (reversed-phase thin-layer chromatography obtained hydrophobicity parameters, RM0 and C0, were established. The good correlations (r > 0.90 were obtained in the majority of relationships. The KOWWINlogP was established as the most suitable hydrophobicity parameter of investigated group of ACE inhibitors with r = 0.981 for correlation with RM0 and r = 0.977 for correlation with C0 parameters (water-methanol mobile phase. Using multiple linear regressions, it was established that application of two selected logP, calculated by different mathematical approaches, led to very good correlation due to the benefits of both calculation methods. The good relationships indicate that the computed logP, with careful selection of method calculation, can be useful in ACE inhibitors lipophilicity evaluation, as high-throughput screening technique.

Oyster Creek cycle 10 nodal model parameter optimization study using PSMS

International Nuclear Information System (INIS)

Dougher, J.D.

1987-01-01

The power shape monitoring system (PSMS) is an on-line core monitoring system that uses a three-dimensional nodal code (NODE-B) to perform nodal power calculations and compute thermal margins. The PSMS contains a parameter optimization function that improves the ability of NODE-B to accurately monitor core power distributions. This functions iterates on the model normalization parameters (albedos and mixing factors) to obtain the best agreement between predicted and measured traversing in-core probe (TIP) reading on a statepoint-by-statepoint basis. Following several statepoint optimization runs, an average set of optimized normalization parameters can be determined and can be implemented into the current or subsequent cycle core model for on-line core monitoring. A statistical analysis of 19 high-power steady-state state-points throughout Oyster Creek cycle 10 operation has shown a consistently poor virgin model performance. The normalization parameters used in the cycle 10 NODE-B model were based on a cycle 8 study, which evaluated only Exxon fuel types. The introduction of General Electric (GE) fuel into cycle 10 (172 assemblies) was a significant fuel/core design change that could have altered the optimum set of normalization parameters. Based on the need to evaluate a potential change in the model normalization parameters for cycle 11 and in an attempt to account for the poor cycle 10 model performance, a parameter optimization study was performed

Lunar tidal acceleration obtained from satellite-derived ocean tide parameters

Science.gov (United States)

Goad, C. C.; Douglas, B. C.

1978-01-01

One hundred sets of mean elements of GEOS-3 computed at 2-day intervals yielded observation equations for the M sub 2 ocean tide from the long periodic variations of the inclination and node of the orbit. The 2nd degree Love number was given the value k sub 2 = 0.30 and the solid tide phase angle was taken to be zero. Combining obtained equations with results for the satellite 1967-92A gives the M sub 2 ocean tide parameter values. Under the same assumption of zero solid tide phase lag, the lunar tidal acceleration was found mostly due to the C sub 22 term in the expansion of the M sub 2 tide with additional small contributions from the 0 sub 1 and N sub 2 tides. Using Lambeck's (1975) estimates for the latter, the obtained acceleration in lunar longitudal in excellent agreement with the most recent determinations from ancient and modern astronomical data.

Non centered minor hysteresis loops evaluation based on exponential parameters transforms of the modified inverse Jiles–Atherton model

International Nuclear Information System (INIS)

Hamimid, M.; Mimoune, S.M.; Feliachi, M.; Atallah, K.

2014-01-01

In this present work, a non centered minor hysteresis loops evaluation is performed using the exponential transforms (ET) of the modified inverse Jiles–Atherton model parameters. This model improves the non centered minor hysteresis loops representation. The parameters of the non centered minor hysteresis loops are obtained from exponential expressions related to the major ones. The parameters of minor loops are obtained by identification using the stochastic optimization method “simulated annealing”. The four parameters of JA model (a,α, k and c) obtained by this transformation are applied only in both ascending and descending branches of the non centered minor hysteresis loops while the major ones are applied to the rest of the cycle. This proposal greatly improves both branches and consequently the minor loops. To validate this model, calculated non-centered minor hysteresis loops are compared with measured ones and good agreements are obtained

A lumped parameter, low dimension model of heat exchanger

International Nuclear Information System (INIS)

Kanoh, Hideaki; Furushoo, Junji; Masubuchi, Masami

1980-01-01

This paper reports on the results of investigation of the distributed parameter model, the difference model, and the model of the method of weighted residuals for heat exchangers. By the method of weighted residuals (MWR), the opposite flow heat exchanger system is approximated by low dimension, lumped parameter model. By assuming constant specific heat, constant density, the same form of tube cross-section, the same form of the surface of heat exchange, uniform flow velocity, the linear relation of heat transfer to flow velocity, liquid heat carrier, and the thermal insulation of liquid from outside, fundamental equations are obtained. The experimental apparatus was made of acrylic resin. The response of the temperature at the exit of first liquid to the variation of the flow rate of second liquid was measured and compared with the models. The MWR model shows good approximation for the low frequency region, and as the number of division increases, good approximation spreads to higher frequency region. (Kato, T.)

Errors and parameter estimation in precipitation-runoff modeling: 1. Theory

Science.gov (United States)

Troutman, Brent M.

1985-01-01

Errors in complex conceptual precipitation-runoff models may be analyzed by placing them into a statistical framework. This amounts to treating the errors as random variables and defining the probabilistic structure of the errors. By using such a framework, a large array of techniques, many of which have been presented in the statistical literature, becomes available to the modeler for quantifying and analyzing the various sources of error. A number of these techniques are reviewed in this paper, with special attention to the peculiarities of hydrologic models. Known methodologies for parameter estimation (calibration) are particularly applicable for obtaining physically meaningful estimates and for explaining how bias in runoff prediction caused by model error and input error may contribute to bias in parameter estimation.

Is there any correlation between model-based perfusion parameters and model-free parameters of time-signal intensity curve on dynamic contrast enhanced MRI in breast cancer patients?

Energy Technology Data Exchange (ETDEWEB)

Yi, Boram; Kang, Doo Kyoung; Kim, Tae Hee [Ajou University School of Medicine, Department of Radiology, Suwon, Gyeonggi-do (Korea, Republic of); Yoon, Dukyong [Ajou University School of Medicine, Department of Biomedical Informatics, Suwon (Korea, Republic of); Jung, Yong Sik; Kim, Ku Sang [Ajou University School of Medicine, Department of Surgery, Suwon (Korea, Republic of); Yim, Hyunee [Ajou University School of Medicine, Department of Pathology, Suwon (Korea, Republic of)

2014-05-15

To find out any correlation between dynamic contrast-enhanced (DCE) model-based parameters and model-free parameters, and evaluate correlations between perfusion parameters with histologic prognostic factors. Model-based parameters (Ktrans, Kep and Ve) of 102 invasive ductal carcinomas were obtained using DCE-MRI and post-processing software. Correlations between model-based and model-free parameters and between perfusion parameters and histologic prognostic factors were analysed. Mean Kep was significantly higher in cancers showing initial rapid enhancement (P = 0.002) and a delayed washout pattern (P = 0.001). Ve was significantly lower in cancers showing a delayed washout pattern (P = 0.015). Kep significantly correlated with time to peak enhancement (TTP) (ρ = -0.33, P < 0.001) and washout slope (ρ = 0.39, P = 0.002). Ve was significantly correlated with TTP (ρ = 0.33, P = 0.002). Mean Kep was higher in tumours with high nuclear grade (P = 0.017). Mean Ve was lower in tumours with high histologic grade (P = 0.005) and in tumours with negative oestrogen receptor status (P = 0.047). TTP was shorter in tumours with negative oestrogen receptor status (P = 0.037). We could acquire general information about the tumour vascular physiology, interstitial space volume and pathologic prognostic factors by analyzing time-signal intensity curve without a complicated acquisition process for the model-based parameters. (orig.)

Assessment of the impact of a parameter estimation method for the Nash Model on selected parameters of a catchment discharge hydrograph

Directory of Open Access Journals (Sweden)

Kołodziejczyk Katarzyna

2017-01-01

Full Text Available An analysis of the usefulness of two parameter calculation methods (N and k parameters for the Nash Model was performed to transform effective rainfall into discharge based on two rainfall episodes gauged at the Kostrze gauging station as well as urban development data for the city of Cracow for 2014 and data obtained from a soil and agriculture map. The methods were the Rao et al. method and the Bajkiewicz-Grabowska method for regression relationships between instantaneous unit hydrograph model parameters and the physiographic parameters of a catchment. Effective rainfall was calculated for each rainfall episode using the SCS-CN method. A direct discharge hydrograph was calculated based on an effective rainfall hyetograph and using the Nash Model. Research has found that both studied methods yield comparable results, which indicates that both methods of effective rainfall transformation into discharge are useful. In addition, it has been shown that the impact of the Nash Model parameter estimation method on discharge hydrographs is minimal.

Updated climatological model predictions of ionospheric and HF propagation parameters

International Nuclear Information System (INIS)

Reilly, M.H.; Rhoads, F.J.; Goodman, J.M.; Singh, M.

1991-01-01

The prediction performances of several climatological models, including the ionospheric conductivity and electron density model, RADAR C, and Ionospheric Communications Analysis and Predictions Program, are evaluated for different regions and sunspot number inputs. Particular attention is given to the near-real-time (NRT) predictions associated with single-station updates. It is shown that a dramatic improvement can be obtained by using single-station ionospheric data to update the driving parameters for an ionospheric model for NRT predictions of f(0)F2 and other ionospheric and HF circuit parameters. For middle latitudes, the improvement extends out thousands of kilometers from the update point to points of comparable corrected geomagnetic latitude. 10 refs

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.

Comparison of three-dimensional parameters of Halo CMEs using three cone models

Science.gov (United States)

Na, H.; Moon, Y.; Jang, S.; Lee, K.

2012-12-01

Halo coronal mass ejections (HCMEs) are a major cause of geomagnetic storms and their three dimensional structures are important for space weather. In this study, we compare three cone models: an elliptical cone model, an ice-cream cone model, and an asymmetric cone model. These models allow us to determine the three dimensional parameters of HCMEs such as radial speed, angular width, and the angle (γ) between sky plane and cone axis. We compare these parameters obtained from three models using 62 well-observed HCMEs observed by SOHO/LASCO from 2001 to 2002. Then we obtain the root mean square error (RMS error) between maximum measured projection speeds and their calculated projection speeds from the cone models. As a result, we find that the radial speeds obtained from the models are well correlated with one another (R > 0.84). The correlation coefficients between angular widths are ranges from 0.04 to 0.53 and those between γ values are from -0.15 to 0.47, which are much smaller than expected. The reason may be due to different assumptions and methods. The RMS errors between the maximum measured projection speeds and the maximum estimated projection speeds of the elliptical cone model, the ice-cream cone model, and the asymmetric cone model are 213 km/s, 254 km/s, and 267 km/s, respectively. And we obtain the correlation coefficients between the location from the models and the flare location (R > 0.75). Finally, we discuss strengths and weaknesses of these models in terms of space weather application.

Application of an Evolutionary Algorithm for Parameter Optimization in a Gully Erosion Model

Energy Technology Data Exchange (ETDEWEB)

Rengers, Francis; Lunacek, Monte; Tucker, Gregory

2016-06-01

Herein we demonstrate how to use model optimization to determine a set of best-fit parameters for a landform model simulating gully incision and headcut retreat. To achieve this result we employed the Covariance Matrix Adaptation Evolution Strategy (CMA-ES), an iterative process in which samples are created based on a distribution of parameter values that evolve over time to better fit an objective function. CMA-ES efficiently finds optimal parameters, even with high-dimensional objective functions that are non-convex, multimodal, and non-separable. We ran model instances in parallel on a high-performance cluster, and from hundreds of model runs we obtained the best parameter choices. This method is far superior to brute-force search algorithms, and has great potential for many applications in earth science modeling. We found that parameters representing boundary conditions tended to converge toward an optimal single value, whereas parameters controlling geomorphic processes are defined by a range of optimal values.

Parameter identification of ZnO surge arrester models based on genetic algorithms

Energy Technology Data Exchange (ETDEWEB)

Bayadi, Abdelhafid [Laboratoire d' Automatique de Setif, Departement d' Electrotechnique, Faculte des Sciences de l' Ingenieur, Universite Ferhat ABBAS de Setif, Route de Bejaia Setif 19000 (Algeria)

2008-07-15

The correct and adequate modelling of ZnO surge arresters characteristics is very important for insulation coordination studies and systems reliability. In this context many researchers addressed considerable efforts to the development of surge arresters models to reproduce the dynamic characteristics observed in their behaviour when subjected to fast front impulse currents. The difficulties with these models reside essentially in the calculation and the adjustment of their parameters. This paper proposes a new technique based on genetic algorithm to obtain the best possible series of parameter values of ZnO surge arresters models. The validity of the predicted parameters is then checked by comparing the predicted results with the experimental results available in the literature. Using the ATP-EMTP package, an application of the arrester model on network system studies is presented and discussed. (author)

Low-dimensional modeling of a driven cavity flow with two free parameters

DEFF Research Database (Denmark)

Jørgensen, Bo Hoffmann; Sørensen, Jens Nørkær; Brøns, Morten

2003-01-01

. By carrying out such a procedure one obtains a low-dimensional model consisting of a reduced set of Ordinary Differential Equations (ODEs) which models the original equations. A technique called Sequential Proper Orthogonal Decomposition (SPOD) is developed to perform decompositions suitable for low...... parameters to appear in the inhomogeneous boundary conditions without the addition of any constraints. This is necessary because both the driving lid and the rotating rod are controlled simultaneously. Apparently, the results reported for this model are the first to be obtained for a low-dimensional model...

Innovation of Methods for Measurement and Modelling of Twisted Pair Parameters

Directory of Open Access Journals (Sweden)

Lukas Cepa

2011-01-01

Full Text Available The goal of this paper is to optimize a measurement methodology for the most accurate broadband modelling of characteristic impedance and other parameters for twisted pairs. Measured values and theirs comparison is presented in this article. Automated measurement facility was implemented at the Department of telecommunication of Faculty of electrical engineering of Czech technical university in Prague. Measurement facility contains RF switches allowing measurements up to 300 MHz or 1GHz. Measured twisted pair’s parameters can be obtained by measurement but for purposes of fundamental characteristics modelling is useful to define functions that model the properties of the twisted pair. Its primary and secondary parameters depend mostly on the frequency. For twisted pair deployment, we are interested in a frequency band range from 1 MHz to 100 MHz.

Genetic Learning of Fuzzy Parameters in Predictive and Decision Support Modelling

Directory of Open Access Journals (Sweden)

Nebot

2012-04-01

Full Text Available In this research a genetic fuzzy system (GFS is proposed that performs discretization parameter learning in the context of the Fuzzy Inductive Reasoning (FIR methodology and the Linguistic Rule FIR (LR-FIR algorithm. The main goal of the GFS is to take advantage of the potentialities of GAs to learn the fuzzification parameters of the FIR and LR-FIR approaches in order to obtain reliable and useful predictive (FIR models and decision support (LR-FIR models. The GFS is evaluated in an e-learning context.

Stochastic Mixed-Effects Parameters Bertalanffy Process, with Applications to Tree Crown Width Modeling

Directory of Open Access Journals (Sweden)

Petras Rupšys

2015-01-01

Full Text Available A stochastic modeling approach based on the Bertalanffy law gained interest due to its ability to produce more accurate results than the deterministic approaches. We examine tree crown width dynamic with the Bertalanffy type stochastic differential equation (SDE and mixed-effects parameters. In this study, we demonstrate how this simple model can be used to calculate predictions of crown width. We propose a parameter estimation method and computational guidelines. The primary goal of the study was to estimate the parameters by considering discrete sampling of the diameter at breast height and crown width and by using maximum likelihood procedure. Performance statistics for the crown width equation include statistical indexes and analysis of residuals. We use data provided by the Lithuanian National Forest Inventory from Scots pine trees to illustrate issues of our modeling technique. Comparison of the predicted crown width values of mixed-effects parameters model with those obtained using fixed-effects parameters model demonstrates the predictive power of the stochastic differential equations model with mixed-effects parameters. All results were implemented in a symbolic algebra system MAPLE.

Identifiability and error minimization of receptor model parameters with PET

International Nuclear Information System (INIS)

Delforge, J.; Syrota, A.; Mazoyer, B.M.

1989-01-01

The identifiability problem and the general framework for experimental design optimization are presented. The methodology is applied to the problem of the receptor-ligand model parameter estimation with dynamic positron emission tomography data. The first attempts to identify the model parameters from data obtained with a single tracer injection led to disappointing numerical results. The possibility of improving parameter estimation using a new experimental design combining an injection of the labelled ligand and an injection of the cold ligand (displacement experiment) has been investigated. However, this second protocol led to two very different numerical solutions and it was necessary to demonstrate which solution was biologically valid. This has been possible by using a third protocol including both a displacement and a co-injection experiment. (authors). 16 refs.; 14 figs

Robust nonlinear autoregressive moving average model parameter estimation using stochastic recurrent artificial neural networks

DEFF Research Database (Denmark)

Chon, K H; Hoyer, D; Armoundas, A A

1999-01-01

In this study, we introduce a new approach for estimating linear and nonlinear stochastic autoregressive moving average (ARMA) model parameters, given a corrupt signal, using artificial recurrent neural networks. This new approach is a two-step approach in which the parameters of the deterministic...... part of the stochastic ARMA model are first estimated via a three-layer artificial neural network (deterministic estimation step) and then reestimated using the prediction error as one of the inputs to the artificial neural networks in an iterative algorithm (stochastic estimation step). The prediction...... error is obtained by subtracting the corrupt signal of the estimated ARMA model obtained via the deterministic estimation step from the system output response. We present computer simulation examples to show the efficacy of the proposed stochastic recurrent neural network approach in obtaining accurate...

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.

Plasma parameters of the cathode spot explosive electron emission cell obtained from the model of liquid-metal jet tearing and electrical explosion

Science.gov (United States)

Tsventoukh, M. M.

2018-05-01

A model has been developed for the explosive electron emission cell pulse of a vacuum discharge cathode spot that describes the ignition and extinction of the explosive pulse. The pulse is initiated due to hydrodynamic tearing of a liquid-metal jet which propagates from the preceding cell crater boundary and draws the ion current from the plasma produced by the preceding explosion. Once the jet neck has been resistively heated to a critical temperature (˜1 eV), the plasma starts expanding and decreasing in density, which corresponds to the extinction phase. Numerical and analytical solutions have been obtained that describe both the time behavior of the pulse plasma parameters and their average values. For the cell plasma, the momentum per transferred charge has been estimated to be some tens of g cm/(s C), which is consistent with the known measurements of ion velocity, ion erosion rate, and specific recoil force. This supports the model of the pressure-gradient-driven plasma acceleration mechanism for the explosive cathode spot cells. The ohmic electric field within the explosive current-carrying plasma has been estimated to be some tens of kV/cm, which is consistent with the known experimental data on cathode potential fall and explosive cell plasma size. This supports the model that assumes the ohmic nature of the cathode potential fall in a vacuum discharge.

Effect of calibration data series length on performance and optimal parameters of hydrological model

Directory of Open Access Journals (Sweden)

Chuan-zhe Li

2010-12-01

Full Text Available In order to assess the effects of calibration data series length on the performance and optimal parameter values of a hydrological model in ungauged or data-limited catchments (data are non-continuous and fragmental in some catchments, we used non-continuous calibration periods for more independent streamflow data for SIMHYD (simple hydrology model calibration. Nash-Sutcliffe efficiency and percentage water balance error were used as performance measures. The particle swarm optimization (PSO method was used to calibrate the rainfall-runoff models. Different lengths of data series ranging from one year to ten years, randomly sampled, were used to study the impact of calibration data series length. Fifty-five relatively unimpaired catchments located all over Australia with daily precipitation, potential evapotranspiration, and streamflow data were tested to obtain more general conclusions. The results show that longer calibration data series do not necessarily result in better model performance. In general, eight years of data are sufficient to obtain steady estimates of model performance and parameters for the SIMHYD model. It is also shown that most humid catchments require fewer calibration data to obtain a good performance and stable parameter values. The model performs better in humid and semi-humid catchments than in arid catchments. Our results may have useful and interesting implications for the efficiency of using limited observation data for hydrological model calibration in different climates.

Patient-specific parameter estimation in single-ventricle lumped circulation models under uncertainty

Science.gov (United States)

Schiavazzi, Daniele E.; Baretta, Alessia; Pennati, Giancarlo; Hsia, Tain-Yen; Marsden, Alison L.

2017-01-01

Summary Computational models of cardiovascular physiology can inform clinical decision-making, providing a physically consistent framework to assess vascular pressures and flow distributions, and aiding in treatment planning. In particular, lumped parameter network (LPN) models that make an analogy to electrical circuits offer a fast and surprisingly realistic method to reproduce the circulatory physiology. The complexity of LPN models can vary significantly to account, for example, for cardiac and valve function, respiration, autoregulation, and time-dependent hemodynamics. More complex models provide insight into detailed physiological mechanisms, but their utility is maximized if one can quickly identify patient specific parameters. The clinical utility of LPN models with many parameters will be greatly enhanced by automated parameter identification, particularly if parameter tuning can match non-invasively obtained clinical data. We present a framework for automated tuning of 0D lumped model parameters to match clinical data. We demonstrate the utility of this framework through application to single ventricle pediatric patients with Norwood physiology. Through a combination of local identifiability, Bayesian estimation and maximum a posteriori simplex optimization, we show the ability to automatically determine physiologically consistent point estimates of the parameters and to quantify uncertainty induced by errors and assumptions in the collected clinical data. We show that multi-level estimation, that is, updating the parameter prior information through sub-model analysis, can lead to a significant reduction in the parameter marginal posterior variance. We first consider virtual patient conditions, with clinical targets generated through model solutions, and second application to a cohort of four single-ventricle patients with Norwood physiology. PMID:27155892

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.

Statistical approach for uncertainty quantification of experimental modal model parameters

DEFF Research Database (Denmark)

Luczak, M.; Peeters, B.; Kahsin, M.

2014-01-01

Composite materials are widely used in manufacture of aerospace and wind energy structural components. These load carrying structures are subjected to dynamic time-varying loading conditions. Robust structural dynamics identification procedure impose tight constraints on the quality of modal models...... represent different complexity levels ranging from coupon, through sub-component up to fully assembled aerospace and wind energy structural components made of composite materials. The proposed method is demonstrated on two application cases of a small and large wind turbine blade........ This paper aims at a systematic approach for uncertainty quantification of the parameters of the modal models estimated from experimentally obtained data. Statistical analysis of modal parameters is implemented to derive an assessment of the entire modal model uncertainty measure. Investigated structures...

Sensitivity of tumor motion simulation accuracy to lung biomechanical modeling approaches and parameters.

Science.gov (United States)

Tehrani, Joubin Nasehi; Yang, Yin; Werner, Rene; Lu, Wei; Low, Daniel; Guo, Xiaohu; Wang, Jing

2015-11-21

Finite element analysis (FEA)-based biomechanical modeling can be used to predict lung respiratory motion. In this technique, elastic models and biomechanical parameters are two important factors that determine modeling accuracy. We systematically evaluated the effects of lung and lung tumor biomechanical modeling approaches and related parameters to improve the accuracy of motion simulation of lung tumor center of mass (TCM) displacements. Experiments were conducted with four-dimensional computed tomography (4D-CT). A Quasi-Newton FEA was performed to simulate lung and related tumor displacements between end-expiration (phase 50%) and other respiration phases (0%, 10%, 20%, 30%, and 40%). Both linear isotropic and non-linear hyperelastic materials, including the neo-Hookean compressible and uncoupled Mooney-Rivlin models, were used to create a finite element model (FEM) of lung and tumors. Lung surface displacement vector fields (SDVFs) were obtained by registering the 50% phase CT to other respiration phases, using the non-rigid demons registration algorithm. The obtained SDVFs were used as lung surface displacement boundary conditions in FEM. The sensitivity of TCM displacement to lung and tumor biomechanical parameters was assessed in eight patients for all three models. Patient-specific optimal parameters were estimated by minimizing the TCM motion simulation errors between phase 50% and phase 0%. The uncoupled Mooney-Rivlin material model showed the highest TCM motion simulation accuracy. The average TCM motion simulation absolute errors for the Mooney-Rivlin material model along left-right, anterior-posterior, and superior-inferior directions were 0.80 mm, 0.86 mm, and 1.51 mm, respectively. The proposed strategy provides a reliable method to estimate patient-specific biomechanical parameters in FEM for lung tumor motion simulation.

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

Uncertainty in dual permeability model parameters for structured soils

Science.gov (United States)

Arora, B.; Mohanty, B. P.; McGuire, J. T.

2012-01-01

Successful application of dual permeability models (DPM) to predict contaminant transport is contingent upon measured or inversely estimated soil hydraulic and solute transport parameters. The difficulty in unique identification of parameters for the additional macropore- and matrix-macropore interface regions, and knowledge about requisite experimental data for DPM has not been resolved to date. Therefore, this study quantifies uncertainty in dual permeability model parameters of experimental soil columns with different macropore distributions (single macropore, and low- and high-density multiple macropores). Uncertainty evaluation is conducted using adaptive Markov chain Monte Carlo (AMCMC) and conventional Metropolis-Hastings (MH) algorithms while assuming 10 out of 17 parameters to be uncertain or random. Results indicate that AMCMC resolves parameter correlations and exhibits fast convergence for all DPM parameters while MH displays large posterior correlations for various parameters. This study demonstrates that the choice of parameter sampling algorithms is paramount in obtaining unique DPM parameters when information on covariance structure is lacking, or else additional information on parameter correlations must be supplied to resolve the problem of equifinality of DPM parameters. This study also highlights the placement and significance of matrix-macropore interface in flow experiments of soil columns with different macropore densities. Histograms for certain soil hydraulic parameters display tri-modal characteristics implying that macropores are drained first followed by the interface region and then by pores of the matrix domain in drainage experiments. Results indicate that hydraulic properties and behavior of the matrix-macropore interface is not only a function of saturated hydraulic conductivity of the macroporematrix interface (Ksa) and macropore tortuosity (lf) but also of other parameters of the matrix and macropore domains.

The application of model with lumped parameters for transient condition analyses of NPP

International Nuclear Information System (INIS)

Stankovic, B.; Stevanovic, V.

1985-01-01

The transient behaviour of NPP Krsko during the accident of pressurizer spray valve stuck open has been simulated y lumped parameters model of the PWR coolant system components, developed at the faculty of Mechanical Engineering, University of Belgrade. The elementary volumes which are characterised by the process and state parameters, and by junctions which are characterised by the geometrical and flow parameters are basic structure of physical model. The process parameters obtained by the model RESI, show qualitative agreement with the measured valves, in a degree in which the actions of reactor safety engineered system and emergency core cooling system are adequately modelled; in spite of the elementary physical model structure and only the modelling of thermal process in reactor core and equilibrium conditions of pressurizer and steam generator. The pressurizer pressure and liquid level predicted by the non-equilibrium pressurizer model SOP show good agreement until the HIPS (high pressure pumps) is activated. (author)

Parameter estimation of Monod model by the Least-Squares method for microalgae Botryococcus Braunii sp

Science.gov (United States)

See, J. J.; Jamaian, S. S.; Salleh, R. M.; Nor, M. E.; Aman, F.

2018-04-01

This research aims to estimate the parameters of Monod model of microalgae Botryococcus Braunii sp growth by the Least-Squares method. Monod equation is a non-linear equation which can be transformed into a linear equation form and it is solved by implementing the Least-Squares linear regression method. Meanwhile, Gauss-Newton method is an alternative method to solve the non-linear Least-Squares problem with the aim to obtain the parameters value of Monod model by minimizing the sum of square error ( SSE). As the result, the parameters of the Monod model for microalgae Botryococcus Braunii sp can be estimated by the Least-Squares method. However, the estimated parameters value obtained by the non-linear Least-Squares method are more accurate compared to the linear Least-Squares method since the SSE of the non-linear Least-Squares method is less than the linear Least-Squares method.

[Temporal and spatial heterogeneity analysis of optimal value of sensitive parameters in ecological process model: The BIOME-BGC model as an example.

Science.gov (United States)

Li, Yi Zhe; Zhang, Ting Long; Liu, Qiu Yu; Li, Ying

2018-01-01

The ecological process models are powerful tools for studying terrestrial ecosystem water and carbon cycle at present. However, there are many parameters for these models, and weather the reasonable values of these parameters were taken, have important impact on the models simulation results. In the past, the sensitivity and the optimization of model parameters were analyzed and discussed in many researches. But the temporal and spatial heterogeneity of the optimal parameters is less concerned. In this paper, the BIOME-BGC model was used as an example. In the evergreen broad-leaved forest, deciduous broad-leaved forest and C3 grassland, the sensitive parameters of the model were selected by constructing the sensitivity judgment index with two experimental sites selected under each vegetation type. The objective function was constructed by using the simulated annealing algorithm combined with the flux data to obtain the monthly optimal values of the sensitive parameters at each site. Then we constructed the temporal heterogeneity judgment index, the spatial heterogeneity judgment index and the temporal and spatial heterogeneity judgment index to quantitatively analyze the temporal and spatial heterogeneity of the optimal values of the model sensitive parameters. The results showed that the sensitivity of BIOME-BGC model parameters was different under different vegetation types, but the selected sensitive parameters were mostly consistent. The optimal values of the sensitive parameters of BIOME-BGC model mostly presented time-space heterogeneity to different degrees which varied with vegetation types. The sensitive parameters related to vegetation physiology and ecology had relatively little temporal and spatial heterogeneity while those related to environment and phenology had generally larger temporal and spatial heterogeneity. In addition, the temporal heterogeneity of the optimal values of the model sensitive parameters showed a significant linear correlation

Using a scalar parameter to trace dislocation evolution in atomistic modeling

Energy Technology Data Exchange (ETDEWEB)

Yang, Jinbo [ORNL; Zhang, Z F [Shenyang National Laboratory for Materials Science; Osetskiy, Yury N [ORNL; Stoller, Roger E [ORNL

2015-01-01

A scalar gamma-parameter is proposed from the Nye tensor. Its maximum value occurs along a dislocation line, either straight or curved, when the coordinate system is purposely chosen. This parameter can be easily obtained from the Nye tensor calculated at each atom in atomistic modeling. Using the gamma-parameter, a fully automated approach is developed to determine core atoms and the Burgers vectors of dislocations simultaneously. The approach is validated by revealing the smallest dislocation loop and by tracing the whole formation process of complicated dislocation networks on the fly.

The influence of technological parameters on the dynamic behavior of "liquid wood" samples obtained by injection molding

Science.gov (United States)

Plavanescu Mazurchevici, Simona; Carausu, Constantin; Comaneci, Radu; Nedelcu, Dumitru

2017-10-01

The plastic products contribute to environmental pollution. Replacing the plastic materials with biodegradable materials with superior properties is an absolute necessity and important research direction for the near future. The first steps in this regard were the creation of composite materials containing natural fibers with positive effects on the environment that have penetrated in different fields. The bioplastics and biocomposites made from natural fibers is a topical solution. The next step was made towards obtaining biodegradable and recyclable materials based on cellulose, lignin and no carcinogens. In this category fall the "liquid wood" with a use up to five times without affecting the mechanical properties. "Liquid wood" is a high quality thermoplastic biocomposite. "Liquid wood" is a biopolymer composite divided in three categories, ARBOFORM®, ARBOBLEND® and ARBOFILL®, which have differed composition in terms of lignin percentage, being delivered by Tecnaro, as granules, [1]. The paper's research was focus on Arboform L V3 Nature and Arboform L V3 Nature reinforced with aramid fiber. In the experimental plan were taken into account six parameters (Dinj - direction of injection [°]; Ttop - melting temperature [°C]; Pinj - injection pressure [MPa] Ss - speed [m/min]; tinj - injection time [s] and tc - cooling time [s]) each with two levels, research carried on by Taguchi methodology. Processing Taguchi method allowed both Taguchi setting work parameters influence on storage modulus and damping as the size and influence their ranking. Experimental research concerning the influence technological parameters on storage modulus of samples obtained by injection from Arboform L V3 Nature yielded an average of 6055MPa and descending order as follows: Trac, Ss, Pinj, Dinj and Ttop. The average of model for reinforced material was 6419MPa and descending order of parameters influence such as: Dinj, Trac, Ttop, tinj, Ss and Pinj.

Derivation of potential model for LiAlO2 by simple and effective optimization of model parameters

International Nuclear Information System (INIS)

Tsuchihira, H.; Oda, T.; Tanaka, S.

2009-01-01

Interatomic potentials of LiAlO 2 were constructed by a simple and effective method. In this method, the model function consists of multiple inverse polynomial functions with an exponential truncation function, and parameters in the potential model can be optimized as a solution of simultaneous linear equations. Potential energies obtained by ab initio calculation are used as fitting targets for model parameter optimization. Lattice constants, elastic properties, defect-formation energy, thermal expansions and the melting point were calculated under the constructed potential models. The results showed good agreement with experimental values and ab initio calculation results, which underscores the validity of the presented method.

Sensitivity of tumor motion simulation accuracy to lung biomechanical modeling approaches and parameters

International Nuclear Information System (INIS)

Tehrani, Joubin Nasehi; Wang, Jing; Yang, Yin; Werner, Rene; Lu, Wei; Low, Daniel; Guo, Xiaohu

2015-01-01

Finite element analysis (FEA)-based biomechanical modeling can be used to predict lung respiratory motion. In this technique, elastic models and biomechanical parameters are two important factors that determine modeling accuracy. We systematically evaluated the effects of lung and lung tumor biomechanical modeling approaches and related parameters to improve the accuracy of motion simulation of lung tumor center of mass (TCM) displacements. Experiments were conducted with four-dimensional computed tomography (4D-CT). A Quasi-Newton FEA was performed to simulate lung and related tumor displacements between end-expiration (phase 50%) and other respiration phases (0%, 10%, 20%, 30%, and 40%). Both linear isotropic and non-linear hyperelastic materials, including the neo-Hookean compressible and uncoupled Mooney–Rivlin models, were used to create a finite element model (FEM) of lung and tumors. Lung surface displacement vector fields (SDVFs) were obtained by registering the 50% phase CT to other respiration phases, using the non-rigid demons registration algorithm. The obtained SDVFs were used as lung surface displacement boundary conditions in FEM. The sensitivity of TCM displacement to lung and tumor biomechanical parameters was assessed in eight patients for all three models. Patient-specific optimal parameters were estimated by minimizing the TCM motion simulation errors between phase 50% and phase 0%. The uncoupled Mooney–Rivlin material model showed the highest TCM motion simulation accuracy. The average TCM motion simulation absolute errors for the Mooney–Rivlin material model along left-right, anterior–posterior, and superior–inferior directions were 0.80 mm, 0.86 mm, and 1.51 mm, respectively. The proposed strategy provides a reliable method to estimate patient-specific biomechanical parameters in FEM for lung tumor motion simulation. (paper)

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.

Lumped-parameter modeling of PWR downcomer and pressurizer for LOCA conditions

International Nuclear Information System (INIS)

Rohatgi, U.S.; Saha, P.; Dubow, A.A.

1978-01-01

Two lumped-parameter models, one for a PWR downcomer and the other for a pressurizer, are presented. The models are based on the transient, nonhomogeneous, drift-flux description of two-phase flow, and are suitable for simulating a hypothetical LOCA condition. Effects of thermal nonequilibrium are incorporated in the downcomer model, whereas the pressurizer model can track the interfaces among various flow regimes. Semiimplicit numerical schemes are used for solution. Encouraging results have been obtained for both the models. (author)

Parameter Estimation

DEFF Research Database (Denmark)

Sales-Cruz, Mauricio; Heitzig, Martina; Cameron, Ian

2011-01-01

of optimisation techniques coupled with dynamic solution of the underlying model. Linear and nonlinear approaches to parameter estimation are investigated. There is also the application of maximum likelihood principles in the estimation of parameters, as well as the use of orthogonal collocation to generate a set......In this chapter the importance of parameter estimation in model development is illustrated through various applications related to reaction systems. In particular, rate constants in a reaction system are obtained through parameter estimation methods. These approaches often require the application...... of algebraic equations as the basis for parameter estimation.These approaches are illustrated using estimations of kinetic constants from reaction system models....

Set-base dynamical parameter estimation and model invalidation for biochemical reaction networks.

Science.gov (United States)

Rumschinski, Philipp; Borchers, Steffen; Bosio, Sandro; Weismantel, Robert; Findeisen, Rolf

2010-05-25

Mathematical modeling and analysis have become, for the study of biological and cellular processes, an important complement to experimental research. However, the structural and quantitative knowledge available for such processes is frequently limited, and measurements are often subject to inherent and possibly large uncertainties. This results in competing model hypotheses, whose kinetic parameters may not be experimentally determinable. Discriminating among these alternatives and estimating their kinetic parameters is crucial to improve the understanding of the considered process, and to benefit from the analytical tools at hand. In this work we present a set-based framework that allows to discriminate between competing model hypotheses and to provide guaranteed outer estimates on the model parameters that are consistent with the (possibly sparse and uncertain) experimental measurements. This is obtained by means of exact proofs of model invalidity that exploit the polynomial/rational structure of biochemical reaction networks, and by making use of an efficient strategy to balance solution accuracy and computational effort. The practicability of our approach is illustrated with two case studies. The first study shows that our approach allows to conclusively rule out wrong model hypotheses. The second study focuses on parameter estimation, and shows that the proposed method allows to evaluate the global influence of measurement sparsity, uncertainty, and prior knowledge on the parameter estimates. This can help in designing further experiments leading to improved parameter estimates.

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.

Modelling of slaughterhouse solid waste anaerobic digestion: determination of parameters and continuous reactor simulation.

Science.gov (United States)

López, Iván; Borzacconi, Liliana

2010-10-01

A model based on the work of Angelidaki et al. (1993) was applied to simulate the anaerobic biodegradation of ruminal contents. In this study, two fractions of solids with different biodegradation rates were considered. A first-order kinetic was used for the easily biodegradable fraction and a kinetic expression that is function of the extracellular enzyme concentration was used for the slowly biodegradable fraction. Batch experiments were performed to obtain an accumulated methane curve that was then used to obtain the model parameters. For this determination, a methodology derived from the "multiple-shooting" method was successfully used. Monte Carlo simulations allowed a confidence range to be obtained for each parameter. Simulations of a continuous reactor were performed using the optimal set of model parameters. The final steady-states were determined as functions of the operational conditions (solids load and residence time). The simulations showed that methane flow peaked at a flow rate of 0.5-0.8 Nm(3)/d/m(reactor)(3) at a residence time of 10-20 days. Simulations allow the adequate selection of operating conditions of a continuous reactor. (c) 2010 Elsevier Ltd. All rights reserved.

Inverse modeling for the determination of hydrogeological parameters of a two-phase system

International Nuclear Information System (INIS)

Finsterle, S.

1993-02-01

Investigations related to the disposal of radioactive wastes in Switzerland consider formations containing natural gas as potential rocks for a repository. Moreover, gas generation in the repository itself may lead to an unsaturated zone of significant extent and impact on the system's performance. The site characterization procedure requires the estimation of hydraulic properties being used as input parameters for a two-phase two-component numerical simulator. In this study, estimates of gas-related formation parameters are obtained by inverse modeling. Based on discrete observations of the system's state, model parameters can be estimated within the framework of a given conceptual model by means of optimization techniques. This study presents the theoretical background that related field data to the model parameters. A parameter estimation procedure is proposed and implemented in a computer code for automatic model calibration. This tool allows identification of key parameters affecting flow of water and gas in porous media. The inverse modeling approach is verified using data from a synthetic laboratory experiment. In addition, the Gas test performed at the Grimsel Test Site is analyzed in order to demonstrate the applicability of the proposed procedure when used with data from a real geologic environment. Estimation of hydrogeologic parameters by automatic model calibration improves the understanding of the two-phase flow processes and therefore increases the reliability of the subsequent simulation runs. (author) figs., tabs., refs

Inverse modeling for the determination of hydrogeological parameters of a two-phase system

International Nuclear Information System (INIS)

Finsterle, S.

1993-01-01

Investigations related to the disposal of radioactive wastes in Switzerland are dealing with formations containing natural gas as potential host rock for a repository. Moreover, gas generation in the repository itself may lead to an unsaturated zone of significant extent and impact on the system's performance. The site characterization procedure requires the estimation of hydraulic properties being used as input parameters for a two-phase two-component numerical simulator. In this study, estimates of gas related formation parameters are obtained by inverse modeling. Based on discrete observations of the system's state, model parameters can be estimated within the framework of a given conceptual model by means of optimization techniques. This study presents the theoretical background that relates field data to the model parameters. A parameter estimation procedure is proposed and implemented in a computer code for automatic model calibration. This tool allows to identify key parameters affecting flow of water and gas in porous media. The inverse modeling approach is verified using data from a synthetic laboratory experiment. In addition, the Gastest performed at the Grimsel Test Site is analyzed in order to demonstrate the applicability of the proposed procedure when used with data from a real geologic environment. Estimation of hydrogeologic parameters by automatic model calibration improves the understanding of the two-phase flow processes and therefore increases the reliability of the subsequent simulation runs. (author) figs., tabs., 100 refs

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.

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)

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.

HIV Model Parameter Estimates from Interruption Trial Data including Drug Efficacy and Reservoir Dynamics

Science.gov (United States)

Luo, Rutao; Piovoso, Michael J.; Martinez-Picado, Javier; Zurakowski, Ryan

2012-01-01

Mathematical models based on ordinary differential equations (ODE) have had significant impact on understanding HIV disease dynamics and optimizing patient treatment. A model that characterizes the essential disease dynamics can be used for prediction only if the model parameters are identifiable from clinical data. Most previous parameter identification studies for HIV have used sparsely sampled data from the decay phase following the introduction of therapy. In this paper, model parameters are identified from frequently sampled viral-load data taken from ten patients enrolled in the previously published AutoVac HAART interruption study, providing between 69 and 114 viral load measurements from 3–5 phases of viral decay and rebound for each patient. This dataset is considerably larger than those used in previously published parameter estimation studies. Furthermore, the measurements come from two separate experimental conditions, which allows for the direct estimation of drug efficacy and reservoir contribution rates, two parameters that cannot be identified from decay-phase data alone. A Markov-Chain Monte-Carlo method is used to estimate the model parameter values, with initial estimates obtained using nonlinear least-squares methods. The posterior distributions of the parameter estimates are reported and compared for all patients. PMID:22815727

Analytical and numerical comparisons of the α parameter obtained from several expressions versus reactivity

International Nuclear Information System (INIS)

Minguez, E.

1983-01-01

A review of the analytical expressions between α-parameter and reactivity (rho) has been done in this article. Finally, several lineal approximations have been obtained for two important points of values. The analytical expansion has been made using only one family of precursors from delayed neutrons; diffusion equation in one energy groups; neutron flux separability with known spatial function. Numerical results have been obtained using typical data from fast and thermal reactors; considering and homogeneous media of Pu 239 in the first case, and U-235 in the second one. (author)

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.

Parameter Estimation in Probit Model for Multivariate Multinomial Response Using SMLE

Directory of Open Access Journals (Sweden)

Jaka Nugraha

2012-02-01

Full Text Available In  the  research  field  of  transportation,  market  research and  politics,  often involving  the  response  of  the multinomial multivariate  observations.  In  this  paper, we discused  a  modeling  of  multivariate  multinomial  responses  using  probit  model.  The estimated  parameters  were  calculated  using Maximum  Likelihood  Estimations  (MLE based  on  the  GHK  simulation.  method  known  as Simulated  Maximum  Likelihood Estimations (SMLE. Likelihood function on the Probit model contains probability values that must be resolved by simulation. By using  the GHK simulation algorithm,  the estimator equation has been obtained for the parameters in the model Probit  Keywords : Probit Model, Newton-Raphson Iteration,  GHK simulator, MLE, simulated log-likelihood

Model atmospheres and parameters of central stars of planetary nebulae

International Nuclear Information System (INIS)

Patriarchi, P.; Cerruti-sola, M.; Perinotto, M.

1989-01-01

Non-LTE hydrogen and helium model atmospheres have been obtained for temperatures and gravities relevant to the central stars of planetary nebulae. Low-resolution and high-resolution observations obtained by the IUE satellite have been used along with optical data to determine Zanstra temperatures of the central stars of NGC 1535, NGC 6210, NGC 7009, IC 418, and IC 4593. Comparison of the observed stellar continuum of these stars with theoretical results allowed further information on the stellar temperature to be derived. The final temperatures are used to calculate accurate stellar parameters. 62 refs

Selecting Sensitive Parameter Subsets in Dynamical Models With Application to Biomechanical System Identification.

Science.gov (United States)

Ramadan, Ahmed; Boss, Connor; Choi, Jongeun; Peter Reeves, N; Cholewicki, Jacek; Popovich, John M; Radcliffe, Clark J

2018-07-01

Estimating many parameters of biomechanical systems with limited data may achieve good fit but may also increase 95% confidence intervals in parameter estimates. This results in poor identifiability in the estimation problem. Therefore, we propose a novel method to select sensitive biomechanical model parameters that should be estimated, while fixing the remaining parameters to values obtained from preliminary estimation. Our method relies on identifying the parameters to which the measurement output is most sensitive. The proposed method is based on the Fisher information matrix (FIM). It was compared against the nonlinear least absolute shrinkage and selection operator (LASSO) method to guide modelers on the pros and cons of our FIM method. We present an application identifying a biomechanical parametric model of a head position-tracking task for ten human subjects. Using measured data, our method (1) reduced model complexity by only requiring five out of twelve parameters to be estimated, (2) significantly reduced parameter 95% confidence intervals by up to 89% of the original confidence interval, (3) maintained goodness of fit measured by variance accounted for (VAF) at 82%, (4) reduced computation time, where our FIM method was 164 times faster than the LASSO method, and (5) selected similar sensitive parameters to the LASSO method, where three out of five selected sensitive parameters were shared by FIM and LASSO methods.

Parameter Identification of the 2-Chlorophenol Oxidation Model Using Improved Differential Search Algorithm

Directory of Open Access Journals (Sweden)

Guang-zhou Chen

2015-01-01

Full Text Available Parameter identification plays a crucial role for simulating and using model. This paper firstly carried out the sensitivity analysis of the 2-chlorophenol oxidation model in supercritical water using the Monte Carlo method. Then, to address the nonlinearity of the model, two improved differential search (DS algorithms were proposed to carry out the parameter identification of the model. One strategy is to adopt the Latin hypercube sampling method to replace the uniform distribution of initial population; the other is to combine DS with simplex method. The results of sensitivity analysis reveal the sensitivity and the degree of difficulty identified for every model parameter. Furthermore, the posteriori probability distribution of parameters and the collaborative relationship between any two parameters can be obtained. To verify the effectiveness of the improved algorithms, the optimization performance of improved DS in kinetic parameter estimation is studied and compared with that of the basic DS algorithm, differential evolution, artificial bee colony optimization, and quantum-behaved particle swarm optimization. And the experimental results demonstrate that the DS with the Latin hypercube sampling method does not present better performance, while the hybrid methods have the advantages of strong global search ability and local search ability and are more effective than the other algorithms.

METHODOLOGY FOR THE ESTIMATION OF PARAMETERS, OF THE MODIFIED BOUC-WEN MODEL

Directory of Open Access Journals (Sweden)

Tomasz HANISZEWSKI

2015-03-01

Full Text Available Bouc-Wen model is theoretical formulation that allows to reflect real hysteresis loop of modeled object. Such object is for example a wire rope, which is present on equipment of crane lifting mechanism. Where adopted modified version of the model has nine parameters. Determination of such a number of parameters is complex and problematic issue. In this article are shown the methodology to identify and sample results of numerical simulations. The results were compared with data obtained on the basis of laboratory tests of ropes [3] and on their basis it was found that there is compliance between results and there is possibility to apply in dynamic systems containing in their structures wire ropes [4].

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.

Application of Powell's optimization method to surge arrester circuit models' parameters

Energy Technology Data Exchange (ETDEWEB)

Christodoulou, C.A.; Stathopulos, I.A. [National Technical University of Athens, School of Electrical and Computer Engineering, 9 Iroon Politechniou St., Zografou Campus, 157 80 Athens (Greece); Vita, V.; Ekonomou, L.; Chatzarakis, G.E. [A.S.PE.T.E. - School of Pedagogical and Technological Education, Department of Electrical Engineering Educators, N. Heraklion, 141 21 Athens (Greece)

2010-08-15

Powell's optimization method has been used for the evaluation of the surge arrester models parameters. The proper modelling of metal-oxide surge arresters and the right selection of equivalent circuit parameters are very significant issues, since quality and reliability of lightning performance studies can be improved with the more efficient representation of the arresters' dynamic behavior. The proposed approach selects optimum arrester model equivalent circuit parameter values, minimizing the error between the simulated peak residual voltage value and this given by the manufacturer. Application of the method in performed on a 120 kV metal oxide arrester. The use of the obtained optimum parameter values reduces significantly the relative error between the simulated and manufacturer's peak residual voltage value, presenting the effectiveness of the method. (author)

Incorporating Parameter Uncertainty in Bayesian Segmentation Models: Application to Hippocampal Subfield Volumetry

DEFF Research Database (Denmark)

Iglesias, J. E.; Sabuncu, M. R.; Van Leemput, Koen

2012-01-01

Many successful segmentation algorithms are based on Bayesian models in which prior anatomical knowledge is combined with the available image information. However, these methods typically have many free parameters that are estimated to obtain point estimates only, whereas a faithful Bayesian anal...

Spatial extrapolation of light use efficiency model parameters to predict gross primary production

Directory of Open Access Journals (Sweden)

Karsten Schulz

2011-12-01

Full Text Available To capture the spatial and temporal variability of the gross primary production as a key component of the global carbon cycle, the light use efficiency modeling approach in combination with remote sensing data has shown to be well suited. Typically, the model parameters, such as the maximum light use efficiency, are either set to a universal constant or to land class dependent values stored in look-up tables. In this study, we employ the machine learning technique support vector regression to explicitly relate the model parameters of a light use efficiency model calibrated at several FLUXNET sites to site-specific characteristics obtained by meteorological measurements, ecological estimations and remote sensing data. A feature selection algorithm extracts the relevant site characteristics in a cross-validation, and leads to an individual set of characteristic attributes for each parameter. With this set of attributes, the model parameters can be estimated at sites where a parameter calibration is not possible due to the absence of eddy covariance flux measurement data. This will finally allow a spatially continuous model application. The performance of the spatial extrapolation scheme is evaluated with a cross-validation approach, which shows the methodology to be well suited to recapture the variability of gross primary production across the study sites.

A Modified Penalty Parameter Approach for Optimal Estimation of UH with Simultaneous Estimation of Infiltration Parameters

Science.gov (United States)

Bhattacharjya, Rajib Kumar

2018-05-01

The unit hydrograph and the infiltration parameters of a watershed can be obtained from observed rainfall-runoff data by using inverse optimization technique. This is a two-stage optimization problem. In the first stage, the infiltration parameters are obtained and the unit hydrograph ordinates are estimated in the second stage. In order to combine this two-stage method into a single stage one, a modified penalty parameter approach is proposed for converting the constrained optimization problem to an unconstrained one. The proposed approach is designed in such a way that the model initially obtains the infiltration parameters and then searches the optimal unit hydrograph ordinates. The optimization model is solved using Genetic Algorithms. A reduction factor is used in the penalty parameter approach so that the obtained optimal infiltration parameters are not destroyed during subsequent generation of genetic algorithms, required for searching optimal unit hydrograph ordinates. The performance of the proposed methodology is evaluated by using two example problems. The evaluation shows that the model is superior, simple in concept and also has the potential for field application.

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)

Assessing Interval Estimation Methods for Hill Model Parameters in a High-Throughput Screening Context (IVIVE meeting)

Science.gov (United States)

The Hill model of concentration-response is ubiquitous in toxicology, perhaps because its parameters directly relate to biologically significant metrics of toxicity such as efficacy and potency. Point estimates of these parameters obtained through least squares regression or maxi...

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...

Two-parameter double-oscillator model of Mathews-Lakshmanan type: Series solutions and supersymmetric partners

International Nuclear Information System (INIS)

Schulze-Halberg, Axel; Wang, Jie

2015-01-01

We obtain series solutions, the discrete spectrum, and supersymmetric partners for a quantum double-oscillator system. Its potential features a superposition of the one-parameter Mathews-Lakshmanan interaction and a one-parameter harmonic or inverse harmonic oscillator contribution. Furthermore, our results are transferred to a generalized Pöschl-Teller model that is isospectral to the double-oscillator system

Two-parameter double-oscillator model of Mathews-Lakshmanan type: Series solutions and supersymmetric partners

Energy Technology Data Exchange (ETDEWEB)

Schulze-Halberg, Axel, E-mail: axgeschu@iun.edu, E-mail: xbataxel@gmail.com [Department of Mathematics and Actuarial Science and Department of Physics, Indiana University Northwest, 3400 Broadway, Gary, Indiana 46408 (United States); Wang, Jie, E-mail: wangjie@iun.edu [Department of Computer Information Systems, Indiana University Northwest, 3400 Broadway, Gary, Indiana 46408 (United States)

2015-07-15

We obtain series solutions, the discrete spectrum, and supersymmetric partners for a quantum double-oscillator system. Its potential features a superposition of the one-parameter Mathews-Lakshmanan interaction and a one-parameter harmonic or inverse harmonic oscillator contribution. Furthermore, our results are transferred to a generalized Pöschl-Teller model that is isospectral to the double-oscillator system.

Parameters identification of photovoltaic models using an improved JAYA optimization algorithm

International Nuclear Information System (INIS)

Yu, Kunjie; Liang, J.J.; Qu, B.Y.; Chen, Xu; Wang, Heshan

2017-01-01

Highlights: • IJAYA algorithm is proposed to identify the PV model parameters efficiently. • A self-adaptive weight is introduced to purposefully adjust the search process. • Experience-based learning strategy is developed to enhance the population diversity. • Chaotic learning method is proposed to refine the quality of the best solution. • IJAYA features the superior performance in identifying parameters of PV models. - Abstract: Parameters identification of photovoltaic (PV) models based on measured current-voltage characteristic curves is significant for the simulation, evaluation, and control of PV systems. To accurately and reliably identify the parameters of different PV models, an improved JAYA (IJAYA) optimization algorithm is proposed in the paper. In IJAYA, a self-adaptive weight is introduced to adjust the tendency of approaching the best solution and avoiding the worst solution at different search stages, which enables the algorithm to approach the promising area at the early stage and implement the local search at the later stage. Furthermore, an experience-based learning strategy is developed and employed randomly to maintain the population diversity and enhance the exploration ability. A chaotic elite learning method is proposed to refine the quality of the best solution in each generation. The proposed IJAYA is used to solve the parameters identification problems of different PV models, i.e., single diode, double diode, and PV module. Comprehensive experiment results and analyses indicate that IJAYA can obtain a highly competitive performance compared with other state-of-the-state algorithms, especially in terms of accuracy and reliability.

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

Lithium-ion Battery Electrothermal Model, Parameter Estimation, and Simulation Environment

Directory of Open Access Journals (Sweden)

Simone Orcioni

2017-03-01

Full Text Available The market for lithium-ion batteries is growing exponentially. The performance of battery cells is growing due to improving production technology, but market request is growing even more rapidly. Modeling and characterization of single cells and an efficient simulation environment is fundamental for the development of an efficient battery management system. The present work is devoted to defining a novel lumped electrothermal circuit of a single battery cell, the extraction procedure of the parameters of the single cell from experiments, and a simulation environment in SystemC-WMS for the simulation of a battery pack. The electrothermal model of the cell was validated against experimental measurements obtained in a climatic chamber. The model is then used to simulate a 48-cell battery, allowing statistical variations among parameters. The different behaviors of the cells in terms of state of charge, current, voltage, or heat flow rate can be observed in the results of the simulation environment.

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...

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.

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.

Hydrological model performance and parameter estimation in the wavelet-domain

Directory of Open Access Journals (Sweden)

B. Schaefli

2009-10-01

Full Text Available This paper proposes a method for rainfall-runoff model calibration and performance analysis in the wavelet-domain by fitting the estimated wavelet-power spectrum (a representation of the time-varying frequency content of a time series of a simulated discharge series to the one of the corresponding observed time series. As discussed in this paper, calibrating hydrological models so as to reproduce the time-varying frequency content of the observed signal can lead to different results than parameter estimation in the time-domain. Therefore, wavelet-domain parameter estimation has the potential to give new insights into model performance and to reveal model structural deficiencies. We apply the proposed method to synthetic case studies and a real-world discharge modeling case study and discuss how model diagnosis can benefit from an analysis in the wavelet-domain. The results show that for the real-world case study of precipitation – runoff modeling for a high alpine catchment, the calibrated discharge simulation captures the dynamics of the observed time series better than the results obtained through calibration in the time-domain. In addition, the wavelet-domain performance assessment of this case study highlights the frequencies that are not well reproduced by the model, which gives specific indications about how to improve the model structure.

Recalibrating disease parameters for increasing realism in modeling epidemics in closed settings

Directory of Open Access Journals (Sweden)

Livio Bioglio

2016-11-01

Full Text Available Abstract Background The homogeneous mixing assumption is widely adopted in epidemic modelling for its parsimony and represents the building block of more complex approaches, including very detailed agent-based models. The latter assume homogeneous mixing within schools, workplaces and households, mostly for the lack of detailed information on human contact behaviour within these settings. The recent data availability on high-resolution face-to-face interactions makes it now possible to assess the goodness of this simplified scheme in reproducing relevant aspects of the infection dynamics. Methods We consider empirical contact networks gathered in different contexts, as well as synthetic data obtained through realistic models of contacts in structured populations. We perform stochastic spreading simulations on these contact networks and in populations of the same size under a homogeneous mixing hypothesis. We adjust the epidemiological parameters of the latter in order to fit the prevalence curve of the contact epidemic model. We quantify the agreement by comparing epidemic peak times, peak values, and epidemic sizes. Results Good approximations of the peak times and peak values are obtained with the homogeneous mixing approach, with a median relative difference smaller than 20 % in all cases investigated. Accuracy in reproducing the peak time depends on the setting under study, while for the peak value it is independent of the setting. Recalibration is found to be linear in the epidemic parameters used in the contact data simulations, showing changes across empirical settings but robustness across groups and population sizes. Conclusions An adequate rescaling of the epidemiological parameters can yield a good agreement between the epidemic curves obtained with a real contact network and a homogeneous mixing approach in a population of the same size. The use of such recalibrated homogeneous mixing approximations would enhance the accuracy and

Recalibrating disease parameters for increasing realism in modeling epidemics in closed settings.

Science.gov (United States)

Bioglio, Livio; Génois, Mathieu; Vestergaard, Christian L; Poletto, Chiara; Barrat, Alain; Colizza, Vittoria

2016-11-14

The homogeneous mixing assumption is widely adopted in epidemic modelling for its parsimony and represents the building block of more complex approaches, including very detailed agent-based models. The latter assume homogeneous mixing within schools, workplaces and households, mostly for the lack of detailed information on human contact behaviour within these settings. The recent data availability on high-resolution face-to-face interactions makes it now possible to assess the goodness of this simplified scheme in reproducing relevant aspects of the infection dynamics. We consider empirical contact networks gathered in different contexts, as well as synthetic data obtained through realistic models of contacts in structured populations. We perform stochastic spreading simulations on these contact networks and in populations of the same size under a homogeneous mixing hypothesis. We adjust the epidemiological parameters of the latter in order to fit the prevalence curve of the contact epidemic model. We quantify the agreement by comparing epidemic peak times, peak values, and epidemic sizes. Good approximations of the peak times and peak values are obtained with the homogeneous mixing approach, with a median relative difference smaller than 20 % in all cases investigated. Accuracy in reproducing the peak time depends on the setting under study, while for the peak value it is independent of the setting. Recalibration is found to be linear in the epidemic parameters used in the contact data simulations, showing changes across empirical settings but robustness across groups and population sizes. An adequate rescaling of the epidemiological parameters can yield a good agreement between the epidemic curves obtained with a real contact network and a homogeneous mixing approach in a population of the same size. The use of such recalibrated homogeneous mixing approximations would enhance the accuracy and realism of agent-based simulations and limit the intrinsic biases of

A One-Step-Ahead Smoothing-Based Joint Ensemble Kalman Filter for State-Parameter Estimation of Hydrological Models

KAUST Repository

El Gharamti, Mohamad; Ait-El-Fquih, Boujemaa; Hoteit, Ibrahim

2015-01-01

The ensemble Kalman filter (EnKF) recursively integrates field data into simulation models to obtain a better characterization of the model’s state and parameters. These are generally estimated following a state-parameters joint augmentation

Innovative parameters obtained for digital analysis of microscopic images to evaluate in vitro hemorheological action of anesthetics

Science.gov (United States)

Alet, Analía. I.; Basso, Sabrina; Delannoy, Marcela; Alet, Nicolás. A.; D'Arrigo, Mabel; Castellini, Horacio V.; Riquelme, Bibiana D.

2015-06-01

Drugs used during anesthesia could enhance microvascular flow disturbance, not only for their systemic cardiovascular actions but also by a direct effect on the microcirculation and in particular on hemorheology. This is particularly important in high-risk surgical patients such as those with vascular disease (diabetes, hypertension, etc.). Therefore, in this work we propose a set of innovative parameters obtained by digital analysis of microscopic images to study the in vitro hemorheological effect of propofol and vecuronium on red blood cell from type 2 diabetic patients compared to healthy donors. Obtained innovative parameters allow quantifying alterations in erythrocyte aggregation, which can increase the in vivo risk of microcapillary obstruction.

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.

Models for the prediction of the cetane index of biofuels obtained from different vegetable oils using their fatty acid composition

International Nuclear Information System (INIS)

Sanchez Borroto, Yisel; Piloto Rodriguez, Ramon; Goyos Perez, Leonardo

2011-01-01

The objective of the present work is to obtain a physical-mathematical model that establishes a relationship between the cetane index of biofuels obtained from different vegetable oils and its composition of essential fatty acid. This model is based on experimental data obtained by the authors of the present work and an experimental data reported by different extracted authors of indexed databases. The adjustment of the coefficients of the model is based on the obtaining of residual minima in the capacity of prediction of the model. Starting from these results it is established a very useful tool for the determination of such an important parameter for the fuel diesel as it is the cetane index obtained from an analysis of chemical composition and not obtained from tests in engines banks, to save time and economic resources. (author)

Inference of reactive transport model parameters using a Bayesian multivariate approach

Science.gov (United States)

Carniato, Luca; Schoups, Gerrit; van de Giesen, Nick

2014-08-01

Parameter estimation of subsurface transport models from multispecies data requires the definition of an objective function that includes different types of measurements. Common approaches are weighted least squares (WLS), where weights are specified a priori for each measurement, and weighted least squares with weight estimation (WLS(we)) where weights are estimated from the data together with the parameters. In this study, we formulate the parameter estimation task as a multivariate Bayesian inference problem. The WLS and WLS(we) methods are special cases in this framework, corresponding to specific prior assumptions about the residual covariance matrix. The Bayesian perspective allows for generalizations to cases where residual correlation is important and for efficient inference by analytically integrating out the variances (weights) and selected covariances from the joint posterior. Specifically, the WLS and WLS(we) methods are compared to a multivariate (MV) approach that accounts for specific residual correlations without the need for explicit estimation of the error parameters. When applied to inference of reactive transport model parameters from column-scale data on dissolved species concentrations, the following results were obtained: (1) accounting for residual correlation between species provides more accurate parameter estimation for high residual correlation levels whereas its influence for predictive uncertainty is negligible, (2) integrating out the (co)variances leads to an efficient estimation of the full joint posterior with a reduced computational effort compared to the WLS(we) method, and (3) in the presence of model structural errors, none of the methods is able to identify the correct parameter values.

Prediction of the intensity and diversity of day-to-day activities among people with schizophrenia using parameters obtained during acute hospitalization.

Science.gov (United States)

Lipskaya-Velikovsky, Lena; Jarus, Tal; Kotler, Moshe

2017-06-01

Participation in day-to-day activities of people with schizophrenia is restricted, causing concern to them, their families, service providers and the communities at large. Participation is a significant component of health and recovery; however, factors predicting participation are still not well established. This study examines whether the parameters obtained during acute hospitalization can predict the intensity and diversity of participation in day-to-day activities six months after discharge. In-patients with chronic schizophrenia (N = 104) were enrolled into the study and assessed for cognitive functioning, functional capacity in instrumental activities of daily living (IADL), and symptoms. Six months after discharge, the intensity and diversity of participation in day-to-day activities were evaluated (N = 70). Multiple correlations were found between parameters obtained during hospitalization and participation diversity, but not participation intensity. The model that is better suited to the prediction of participation diversity contains cognitive ability of construction, negative symptoms and number of previous hospitalizations. The total explained variance is 37.8% (F 3,66  =   14.99, p process for the prediction of participation diversity in day-to-day activities six months after discharge. Participation diversity is best predicted through a set of factors reflecting personal and environmental indicators. Implications for rehabilitation Results of in-patient evaluations can predict the diversity of participation in day-to-day activities six months after discharge. Higher prediction of participation diversity is obtained using a holistic evaluation model that includes assessments for cognitive abilities, negative symptoms severity and number of hospitalizations.

Stability switches, Hopf bifurcation and chaos of a neuron model with delay-dependent parameters

International Nuclear Information System (INIS)

Xu, X.; Hu, H.Y.; Wang, H.L.

2006-01-01

It is very common that neural network systems usually involve time delays since the transmission of information between neurons is not instantaneous. Because memory intensity of the biological neuron usually depends on time history, some of the parameters may be delay dependent. Yet, little attention has been paid to the dynamics of such systems. In this Letter, a detailed analysis on the stability switches, Hopf bifurcation and chaos of a neuron model with delay-dependent parameters is given. Moreover, the direction and the stability of the bifurcating periodic solutions are obtained by the normal form theory and the center manifold theorem. It shows that the dynamics of the neuron model with delay-dependent parameters is quite different from that of systems with delay-independent parameters only

Uncertainty analyses of the calibrated parameter values of a water quality model

Science.gov (United States)

Rode, M.; Suhr, U.; Lindenschmidt, K.-E.

2003-04-01

For river basin management water quality models are increasingly used for the analysis and evaluation of different management measures. However substantial uncertainties exist in parameter values depending on the available calibration data. In this paper an uncertainty analysis for a water quality model is presented, which considers the impact of available model calibration data and the variance of input variables. The investigation was conducted based on four extensive flowtime related longitudinal surveys in the River Elbe in the years 1996 to 1999 with varying discharges and seasonal conditions. For the model calculations the deterministic model QSIM of the BfG (Germany) was used. QSIM is a one dimensional water quality model and uses standard algorithms for hydrodynamics and phytoplankton dynamics in running waters, e.g. Michaelis Menten/Monod kinetics, which are used in a wide range of models. The multi-objective calibration of the model was carried out with the nonlinear parameter estimator PEST. The results show that for individual flow time related measuring surveys very good agreements between model calculation and measured values can be obtained. If these parameters are applied to deviating boundary conditions, substantial errors in model calculation can occur. These uncertainties can be decreased with an increased calibration database. More reliable model parameters can be identified, which supply reasonable results for broader boundary conditions. The extension of the application of the parameter set on a wider range of water quality conditions leads to a slight reduction of the model precision for the specific water quality situation. Moreover the investigations show that highly variable water quality variables like the algal biomass always allow a smaller forecast accuracy than variables with lower coefficients of variation like e.g. nitrate.

Building X-ray pulsar timing model without the use of radio parameters

Science.gov (United States)

Sun, Hai-feng; Sun, Xiong; Fang, Hai-yan; Shen, Li-rong; Cong, Shao-peng; Liu, Yan-ming; Li, Xiao-ping; Bao, Wei-min

2018-02-01

This paper develops a timing solution for the X-ray pulsar timing model without the use of the initial radio model parameters. First, we address the problem of phase ambiguities for the pre-fit residuals in the construction of pulsar timing model. To improve the estimation accuracy of the pulse time of arrival (TOA), we have deduced the general form of test statistics in Fourier transform, and discussed their estimation performances. Meanwhile, a fast maximum likelihood (FML) technique is presented to estimate the pulse TOA, which outperforms cross correlation (CC) estimator and exhibits a performance comparable with maximum likelihood (ML) estimator in spite of a much less reduced computational complexity. Depending on the strategy of the difference minimum of pre-fit residuals, we present an effective forced phase-connected technique to achieve initial model parameters. Then, we use the observations with the Rossi X-Ray Timing Explorer (RXTE) and X-ray pulsar navigation-I (XPNAV-1) satellites for experimental studies, and discuss main differences for the root mean square (RMS) residuals calculated with the X-ray and radio ephemerides. Finally, a chi-square value (CSV) of pulse profiles is presented as a complementary indicator to the RMS residuals for evaluating the model parameters. The results show that the proposed timing solution is valid and effective, and the obtained model parameters can be a reasonable alternative to the radio ephemeris.

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.

Summary of the DREAM8 Parameter Estimation Challenge: Toward Parameter Identification for Whole-Cell Models.

Directory of Open Access Journals (Sweden)

Jonathan R Karr

2015-05-01

Full Text Available Whole-cell models that explicitly represent all cellular components at the molecular level have the potential to predict phenotype from genotype. However, even for simple bacteria, whole-cell models will contain thousands of parameters, many of which are poorly characterized or unknown. New algorithms are needed to estimate these parameters and enable researchers to build increasingly comprehensive models. We organized the Dialogue for Reverse Engineering Assessments and Methods (DREAM 8 Whole-Cell Parameter Estimation Challenge to develop new parameter estimation algorithms for whole-cell models. We asked participants to identify a subset of parameters of a whole-cell model given the model's structure and in silico "experimental" data. Here we describe the challenge, the best performing methods, and new insights into the identifiability of whole-cell models. We also describe several valuable lessons we learned toward improving future challenges. Going forward, we believe that collaborative efforts supported by inexpensive cloud computing have the potential to solve whole-cell model parameter estimation.

Microscopic calculation of the Majorana parameters of the interacting boson model for the Hg isotopes

Energy Technology Data Exchange (ETDEWEB)

Druce, C.H.; Barrett, B.R. (Arizona Univ., Tucson (USA). Dept. of Physics); Pittel, S. (Delaware Univ., Newark (USA). Bartol Research Foundation); Duval, P.D. (BEERS Associates, Reston, VA (USA))

1985-07-11

The parameters of the Majorana interaction of the neutron-proton interacting boson model are calculated for the Hg isotopes. The calculations utilize the Otsuka-Arima-Iachello mapping procedure and also lead to predictions for the other boson parameters. The resulting spectra are compared with experimental spectra and those obtained from phenomenological fits.

Microscopic calculation of the Majorana parameters of the interacting boson model for the Hg isotopes

Science.gov (United States)

Druce, C. H.; Pittel, S.; Barrett, B. R.; Duval, P. D.

1985-07-01

The parameters of the Majorana interaction of the neutron-proton interacting boson model are calculated for the Hg isotopes. The calculations utilize the Otsuka-Arima-Iachello mapping procedure and also lead to predictions for the other boson parameters. The resulting spectra are compared with experimental spectra and those obtained from phenomenological fits.

Modeling phosphorus in the Lake Allatoona watershed using SWAT: I. Developing phosphorus parameter values.

Science.gov (United States)

Radcliffe, D E; Lin, Z; Risse, L M; Romeis, J J; Jackson, C R

2009-01-01

Lake Allatoona is a large reservoir north of Atlanta, GA, that drains an area of about 2870 km2 scheduled for a phosphorus (P) total maximum daily load (TMDL). The Soil and Water Assessment Tool (SWAT) model has been widely used for watershed-scale modeling of P, but there is little guidance on how to estimate P-related parameters, especially those related to in-stream P processes. In this paper, methods are demonstrated to individually estimate SWAT soil-related P parameters and to collectively estimate P parameters related to stream processes. Stream related parameters were obtained using the nutrient uptake length concept. In a manner similar to experiments conducted by stream ecologists, a small point source is simulated in a headwater sub-basin of the SWAT models, then the in-stream parameter values are adjusted collectively to get an uptake length of P similar to the values measured in the streams in the region. After adjusting the in-stream parameters, the P uptake length estimated in the simulations ranged from 53 to 149 km compared to uptake lengths measured by ecologists in the region of 11 to 85 km. Once the a priori P-related parameter set was developed, the SWAT models of main tributaries to Lake Allatoona were calibrated for daily transport. Models using SWAT P parameters derived from the methods in this paper outperformed models using default parameter values when predicting total P (TP) concentrations in streams during storm events and TP annual loads to Lake Allatoona.

Mathematical modeling and evaluation of radionuclide transport parameters from the ANL Laboratory Analog Program

International Nuclear Information System (INIS)

Chen, B.C.J.; Hull, J.R.; Seitz, M.G.; Sha, W.T.; Shah, V.L.; Soo, S.L.

1984-07-01

Computer model simulation is required to evaluate the performance of proposed or future high-level radioactive waste geological repositories. However, the accuracy of a model in predicting the real situation depends on how well the values of the transport properties are prescribed as input parameters. Knowledge of transport parameters is therefore essential. We have modeled ANL's Experiment Analog Program which was designed to simulate long-term radwaste migration process by groundwater flowing through a high-level radioactive waste repository. Using this model and experimental measurements, we have evaluated neptunium (actinide) deposition velocity and analyzed the complex phenomena of simultaneous deposition, erosion, and reentrainment of bentonite when groundwater is flowing through a narrow crack in a basalt rock. The present modeling demonstrates that we can obtain the values of transport parameters, as added information without any additional cost, from the available measurements of laboratory analog experiments. 8 figures, 3 tables

Coastal aquifer management under parameter uncertainty: Ensemble surrogate modeling based simulation-optimization

Science.gov (United States)

Janardhanan, S.; Datta, B.

2011-12-01

Surrogate models are widely used to develop computationally efficient simulation-optimization models to solve complex groundwater management problems. Artificial intelligence based models are most often used for this purpose where they are trained using predictor-predictand data obtained from a numerical simulation model. Most often this is implemented with the assumption that the parameters and boundary conditions used in the numerical simulation model are perfectly known. However, in most practical situations these values are uncertain. Under these circumstances the application of such approximation surrogates becomes limited. In our study we develop a surrogate model based coupled simulation optimization methodology for determining optimal pumping strategies for coastal aquifers considering parameter uncertainty. An ensemble surrogate modeling approach is used along with multiple realization optimization. The methodology is used to solve a multi-objective coastal aquifer management problem considering two conflicting objectives. Hydraulic conductivity and the aquifer recharge are considered as uncertain values. Three dimensional coupled flow and transport simulation model FEMWATER is used to simulate the aquifer responses for a number of scenarios corresponding to Latin hypercube samples of pumping and uncertain parameters to generate input-output patterns for training the surrogate models. Non-parametric bootstrap sampling of this original data set is used to generate multiple data sets which belong to different regions in the multi-dimensional decision and parameter space. These data sets are used to train and test multiple surrogate models based on genetic programming. The ensemble of surrogate models is then linked to a multi-objective genetic algorithm to solve the pumping optimization problem. Two conflicting objectives, viz, maximizing total pumping from beneficial wells and minimizing the total pumping from barrier wells for hydraulic control of

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.

Parameters Determination of Yoshida Uemori Model Through Optimization Process of Cyclic Tension-Compression Test and V-Bending Springback

Directory of Open Access Journals (Sweden)

Serkan Toros

Full Text Available Abstract In recent years, the studies on the enhancement of the prediction capability of the sheet metal forming simulations have increased remarkably. Among the used models in the finite element simulations, the yield criteria and hardening models have a great importance for the prediction of the formability and springback. The required model parameters are determined by using the several test results, i.e. tensile, compression, biaxial stretching tests (bulge test and cyclic tests (tension-compression. In this study, the Yoshida-Uemori (combined isotropic and kinematic hardening model is used to determine the performance of the springback prediction. The model parameters are determined by the optimization processes of the cyclic test by finite element simulations. However, in the study besides the cyclic tests, the model parameters are also evaluated by the optimization process of both cyclic and V-die bending simulations. The springback angle predictions with the model parameters obtained by the optimization of both cyclic and V-die bending simulations are found to mimic the experimental results in a better way than those obtained from only cyclic tests. However, the cyclic simulation results are found to be close enough to the experimental results.

Determination of the Nonlinearity Parameter in the TNM Model of Structural Recovery

Science.gov (United States)

Bari, Rozana; Simon, Sindee

Structural recovery of non-equilibrium glassy materials takes place by evolution of volume and enthalpy as the glass attempts to reach to equilibrium. Structural recovery is nonlinear, nonexponential, and depends on thermal history and the process can be described by phenomenological models of structural recovery, such as the Tool-Narayanaswamy-Moynihan (TNM) and the Kovacs-Aklonis-Hutchinson-Ramos (KAHR) models. The goal of the present work is to analyze methods to determine the nonlinearity parameter x and activation energy Δh/R. The methods to determine x includes the inflectional analysis, time-temperature superposition, and two-step temperature jump methods. The activation energy Δh/R can also be obtained by the first two methods. The TNM model is used to simulate structural recovery data, which are then used to test the accuracy of the methods to determine x and Δh/R, with a particular interest in data obtained after cooling at high rates as can be obtained in the Flash DSC. The nonlinearity parameter x by the inflectional analysis and two-step temperature methods are accurate for exponential recovery. However, for real systems with nonexponential relaxation, methods to determine x are not reliable. The activation energy is well estimated by both the time-temperature superposition and inflectional analysis methods, with the former being slightly better.

THE INFLUENCE OF CONVERSION MODEL CHOICE FOR EROSION RATE ESTIMATION AND THE SENSITIVITY OF THE RESULTS TO CHANGES IN THE MODEL PARAMETER

Directory of Open Access Journals (Sweden)

Nita Suhartini

2010-06-01

Full Text Available A study of soil erosion rates had been done on a slightly and long slope of cultivated area in Ciawi - Bogor, using 137Cs technique. The objective of the present study was to evaluate the applicability of the 137Cs technique in obtaining spatially distributed information of soil redistribution at small catchment. This paper reports the result of the choice of conversion model for erosion rate estimates and the sensitive of the changes in the model parameter. For this purpose, small site was selected, namely landuse I (LU-I. The top of a slope was chosen as a reference site. The erosion/deposit rate of individual sampling points was estimated using the conversion models, namely Proportional Model (PM, Mass Balance Model 1 (MBM1 and Mass Balance Model 2 (MBM2. A comparison of the conversion models showed that the lowest value is obtained by the PM. The MBM1 gave values closer to MBM2, but MBM2 gave a reliable values. In this study, a sensitivity analysis suggest that the conversion models are sensitive to changes in parameters that depend on the site conditions, but insensitive to changes in  parameters that interact to the onset of 137Cs fallout input.   Keywords: soil erosion, environmental radioisotope, cesium

Comparison of ionospheric profile parameters with IRI-2012 model over Jicamarca

Science.gov (United States)

Bello, S. A.; Abdullah, M.; Hamid, N. S. A.; Reinisch, B. W.

2017-05-01

We used the hourly ionogram data obtained from Jicamarca station (12° S, 76.9° W, dip latitude: 1.0° N) an equatorial region to study the variation of the electron density profile parameters: maximum height of F2-layer (hmF2), bottomside thickness (B0) and shape (B1) parameter of F-layer. The period of study is for the year 2010 (solar minimum period).The diurnal monthly averages of these parameters are compared with the updated IRI-2012 model. The results show that hmF2 is highest during the daytime than nighttime. The variation in hmF2 was observed to modulate the thickness of the bottomside F2-layer. The observed hmF2 and B0 post-sunset peak is as result of the upward drift velocity of ionospheric plasma. We found a close agreement between IRI-CCIR hmF2 model and observed hmF2 during 0000-0700 LT while outside this period the model predictions deviate significantly with the observational values. Significant discrepancies are observed between the IRI model options for B0 and the observed B0 values. Specifically, the modeled values do not show B0 post-sunset peak. A fairly good agreement was observed between the observed B1 and IRI model options (ABT-2009 and Bill 2000) for B1.

Comparison of ionospheric profile parameters with IRI-2012 model over Jicamarca

International Nuclear Information System (INIS)

Bello, S A; Abdullah, M; Hamid, N S A; Reinisch, B W

2017-01-01

We used the hourly ionogram data obtained from Jicamarca station (12° S, 76.9° W, dip latitude: 1.0° N) an equatorial region to study the variation of the electron density profile parameters: maximum height of F2-layer ( hm F2), bottomside thickness ( B0 ) and shape ( B1 ) parameter of F-layer. The period of study is for the year 2010 (solar minimum period).The diurnal monthly averages of these parameters are compared with the updated IRI-2012 model. The results show that hm F2 is highest during the daytime than nighttime. The variation in hmF2 was observed to modulate the thickness of the bottomside F2-layer. The observed hm F2 and B0 post-sunset peak is as result of the upward drift velocity of ionospheric plasma. We found a close agreement between IRI-CCIR hm F2 model and observed hm F2 during 0000-0700 LT while outside this period the model predictions deviate significantly with the observational values. Significant discrepancies are observed between the IRI model options for B0 and the observed B0 values. Specifically, the modeled values do not show B0 post-sunset peak. A fairly good agreement was observed between the observed B1 and IRI model options (ABT-2009 and Bill 2000) for B1 . (paper)

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.

Parameter Identification with the Random Perturbation Particle Swarm Optimization Method and Sensitivity Analysis of an Advanced Pressurized Water Reactor Nuclear Power Plant Model for Power Systems

Directory of Open Access Journals (Sweden)

Li Wang

2017-02-01

Full Text Available The ability to obtain appropriate parameters for an advanced pressurized water reactor (PWR unit model is of great significance for power system analysis. The attributes of that ability include the following: nonlinear relationships, long transition time, intercoupled parameters and difficult obtainment from practical test, posed complexity and difficult parameter identification. In this paper, a model and a parameter identification method for the PWR primary loop system were investigated. A parameter identification process was proposed, using a particle swarm optimization (PSO algorithm that is based on random perturbation (RP-PSO. The identification process included model variable initialization based on the differential equations of each sub-module and program setting method, parameter obtainment through sub-module identification in the Matlab/Simulink Software (Math Works Inc., Natick, MA, USA as well as adaptation analysis for an integrated model. A lot of parameter identification work was carried out, the results of which verified the effectiveness of the method. It was found that the change of some parameters, like the fuel temperature and coolant temperature feedback coefficients, changed the model gain, of which the trajectory sensitivities were not zero. Thus, obtaining their appropriate values had significant effects on the simulation results. The trajectory sensitivities of some parameters in the core neutron dynamic module were interrelated, causing the parameters to be difficult to identify. The model parameter sensitivity could be different, which would be influenced by the model input conditions, reflecting the parameter identifiability difficulty degree for various input conditions.

Revised Parameters for the AMOEBA Polarizable Atomic Multipole Water Model

Science.gov (United States)

Pande, Vijay S.; Head-Gordon, Teresa; Ponder, Jay W.

2016-01-01

A set of improved parameters for the AMOEBA polarizable atomic multipole water model is developed. The protocol uses an automated procedure, ForceBalance, to adjust model parameters to enforce agreement with ab initio-derived results for water clusters and experimentally obtained data for a variety of liquid phase properties across a broad temperature range. The values reported here for the new AMOEBA14 water model represent a substantial improvement over the previous AMOEBA03 model. The new AMOEBA14 water model accurately predicts the temperature of maximum density and qualitatively matches the experimental density curve across temperatures ranging from 249 K to 373 K. Excellent agreement is observed for the AMOEBA14 model in comparison to a variety of experimental properties as a function of temperature, including the 2nd virial coefficient, enthalpy of vaporization, isothermal compressibility, thermal expansion coefficient and dielectric constant. The viscosity, self-diffusion constant and surface tension are also well reproduced. In comparison to high-level ab initio results for clusters of 2 to 20 water molecules, the AMOEBA14 model yields results similar to the AMOEBA03 and the direct polarization iAMOEBA models. With advances in computing power, calibration data, and optimization techniques, we recommend the use of the AMOEBA14 water model for future studies employing a polarizable water model. PMID:25683601

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.

Impact of the calibration period on the conceptual rainfall-runoff model parameter estimates

Science.gov (United States)

Todorovic, Andrijana; Plavsic, Jasna

2015-04-01

. Correlation coefficients among optimised model parameters and total precipitation P, mean temperature T and mean flow Q are calculated to give an insight into parameter dependence on the hydrometeorological drivers. The results reveal high sensitivity of almost all model parameters towards calibration period. The highest variability is displayed by the refreezing coefficient, water holding capacity, and temperature gradient. The only statistically significant (decreasing) trend is detected in the evapotranspiration reduction threshold. Statistically significant correlation is detected between the precipitation gradient and precipitation depth, and between the time-area histogram base and flows. All other correlations are not statistically significant, implying that changes in optimised parameters cannot generally be linked to the changes in P, T or Q. As for the model performance, the model reproduces the observed runoff satisfactorily, though the runoff is slightly overestimated in wet periods. The Nash-Sutcliffe efficiency coefficient (NSE) ranges from 0.44 to 0.79. Higher NSE values are obtained over wetter periods, what is supported by statistically significant correlation between NSE and flows. Overall, no systematic variations in parameters or in model performance are detected. Parameter variability may therefore rather be attributed to errors in data or inadequacies in the model structure. Further research is required to examine the impact of the calibration strategy or model structure on the variability in optimised parameters in time.

Accelerating parameter identification of proton exchange membrane fuel cell model with ranking-based differential evolution

International Nuclear Information System (INIS)

Gong, Wenyin; Cai, Zhihua

2013-01-01

Parameter identification of PEM (proton exchange membrane) fuel cell model is a very active area of research. Generally, it can be treated as a numerical optimization problem with complex nonlinear and multi-variable features. DE (differential evolution), which has been successfully used in various fields, is a simple yet efficient evolutionary algorithm for global numerical optimization. In this paper, with the objective of accelerating the process of parameter identification of PEM fuel cell models and reducing the necessary computational efforts, we firstly present a generic and simple ranking-based mutation operator for the DE algorithm. Then, the ranking-based mutation operator is incorporated into five highly-competitive DE variants to solve the PEM fuel cell model parameter identification problems. The main contributions of this work are the proposed ranking-based DE variants and their application to the parameter identification problems of PEM fuel cell models. Experiments have been conducted by using both the simulated voltage–current data and the data obtained from the literature to validate the performance of our approach. The results indicate that the ranking-based DE methods provide better results with respect to the solution quality, the convergence rate, and the success rate compared with their corresponding original DE methods. In addition, the voltage–current characteristics obtained by our approach are in good agreement with the original voltage–current curves in all cases. - Highlights: • A simple and generic ranking-based mutation operator is presented in this paper. • Several DE (differential evolution) variants are used to solve the parameter identification of PEMFC (proton exchange membrane fuel cells) model. • Results show that our method accelerates the process of parameter identification. • The V–I characteristics are in very good agreement with experimental data

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.

A new analytical method for estimating lumped parameter constants of linear viscoelastic models from strain rate tests

Science.gov (United States)

Mattei, G.; Ahluwalia, A.

2018-04-01

We introduce a new function, the apparent elastic modulus strain-rate spectrum, E_{app} ( \\dot{ɛ} ), for the derivation of lumped parameter constants for Generalized Maxwell (GM) linear viscoelastic models from stress-strain data obtained at various compressive strain rates ( \\dot{ɛ}). The E_{app} ( \\dot{ɛ} ) function was derived using the tangent modulus function obtained from the GM model stress-strain response to a constant \\dot{ɛ} input. Material viscoelastic parameters can be rapidly derived by fitting experimental E_{app} data obtained at different strain rates to the E_{app} ( \\dot{ɛ} ) function. This single-curve fitting returns similar viscoelastic constants as the original epsilon dot method based on a multi-curve global fitting procedure with shared parameters. Its low computational cost permits quick and robust identification of viscoelastic constants even when a large number of strain rates or replicates per strain rate are considered. This method is particularly suited for the analysis of bulk compression and nano-indentation data of soft (bio)materials.

Incorporating model parameter uncertainty into inverse treatment planning

International Nuclear Information System (INIS)

Lian Jun; Xing Lei

2004-01-01

Radiobiological treatment planning depends not only on the accuracy of the models describing the dose-response relation of different tumors and normal tissues but also on the accuracy of tissue specific radiobiological parameters in these models. Whereas the general formalism remains the same, different sets of model parameters lead to different solutions and thus critically determine the final plan. Here we describe an inverse planning formalism with inclusion of model parameter uncertainties. This is made possible by using a statistical analysis-based frameset developed by our group. In this formalism, the uncertainties of model parameters, such as the parameter a that describes tissue-specific effect in the equivalent uniform dose (EUD) model, are expressed by probability density function and are included in the dose optimization process. We found that the final solution strongly depends on distribution functions of the model parameters. Considering that currently available models for computing biological effects of radiation are simplistic, and the clinical data used to derive the models are sparse and of questionable quality, the proposed technique provides us with an effective tool to minimize the effect caused by the uncertainties in a statistical sense. With the incorporation of the uncertainties, the technique has potential for us to maximally utilize the available radiobiology knowledge for better IMRT treatment

Modelling pesticide leaching under climate change: parameter vs. climate input uncertainty

Directory of Open Access Journals (Sweden)

K. Steffens

2014-02-01

Full Text Available Assessing climate change impacts on pesticide leaching requires careful consideration of different sources of uncertainty. We investigated the uncertainty related to climate scenario input and its importance relative to parameter uncertainty of the pesticide leaching model. The pesticide fate model MACRO was calibrated against a comprehensive one-year field data set for a well-structured clay soil in south-western Sweden. We obtained an ensemble of 56 acceptable parameter sets that represented the parameter uncertainty. Nine different climate model projections of the regional climate model RCA3 were available as driven by different combinations of global climate models (GCM, greenhouse gas emission scenarios and initial states of the GCM. The future time series of weather data used to drive the MACRO model were generated by scaling a reference climate data set (1970–1999 for an important agricultural production area in south-western Sweden based on monthly change factors for 2070–2099. 30 yr simulations were performed for different combinations of pesticide properties and application seasons. Our analysis showed that both the magnitude and the direction of predicted change in pesticide leaching from present to future depended strongly on the particular climate scenario. The effect of parameter uncertainty was of major importance for simulating absolute pesticide losses, whereas the climate uncertainty was relatively more important for predictions of changes of pesticide losses from present to future. The climate uncertainty should be accounted for by applying an ensemble of different climate scenarios. The aggregated ensemble prediction based on both acceptable parameterizations and different climate scenarios has the potential to provide robust probabilistic estimates of future pesticide losses.

Influence of sintering parameters in the ferroelectric properties os strontium bismuth tantalate samples obtained by oxide mixture

International Nuclear Information System (INIS)

Souza, R.R. de; Pereira, A.S.; Sousa, V.C.; Egea, J.R.J.

2012-01-01

The family of compounds layered-type perovskite, know as Aurivilius presents great alternative not only by the absence of lead in the composition, but because the polarization retention, replacing PZT in FeRAM devices. The strontium bismuth tantalate (SrBi 2 Ta 2 O 9 ) or SBT is ferroelectric material that has attracted considerable interest, since it has high fatigue resistance, supporting high hysteresis loops, with the change in polarization.Checking polarization and depolarization currents stimulated by temperature it is possible to obtain, for example, information about the nature of charges and about the activation energy for the process of dielectric relaxation. For analysis of ferroelectric properties of this compound, it is essential to obtain specimens with a relative density around 95%. Thus, it is important the optimization of the sintering process in order to obtain a ceramic body with a high densification. The influence of sintering parameters to obtain SrBi 2 Ta 2 O 9 in the polarization properties and in the microstructure of sintered samples was investigated by thermostimulated currents and electronic microscopy, respectively. Results show that variation of these parameters may cause changes in the ferroelectric properties of the material. (author)

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.

Study of Λ parameters and crossover phenomena in SU(N) x SU(N) sigma models in two dimensions

International Nuclear Information System (INIS)

Shigemitsu, J.; Kogut, J.B.

1981-01-01

The spin system analogues of recent studies of the string tension and Λ parameters of SU(N) gauge theories in 4 dimensions are carried out for the SU(N) x SU(N) and O(N) models in 2 dimensions. The relations between the Λ parameters of both the Euclidean and Hamiltonian formulation of the lattice models and the Λ parameter of the continuum models are obtained. The one loop finite renormalization of the speed of light in the lattice Hamiltonian formulations of the O(N) and SU(N) x SU(N) models is calculated. Strong coupling calculations of the mass gaps of these spin models are done for all N and the constants of proportionality between the gap and the Λ parameter of the continuum models are obtained. These results are contrasted with similar calculations for the SU(N) gauge models in 3+1 dimensions. Identifying suitable coupling constants for discussing the N → infinity limits, the numerical results suggest that the crossover from weak to strong coupling in the lattice O(N) models becomes less abrupt as N increases while the crossover for the SU(N) x SU(N) models becomes more abrupt. The crossover in SU(N) gauge theories also becomes more abrupt with increasing N, however, at an even greater rate than in the SU(N) x SU(N) spin models

Modelled basic parameters for semi-industrial irradiation plant design

International Nuclear Information System (INIS)

Mangussi, J.

2009-01-01

The basic parameters of an irradiation plant design are the total activity, the product uniformity ratio and the efficiency process. The target density, the minimum dose required and the throughput depends on the use to which the irradiator will be put at. In this work, a model for calculating the specific dose rate at several depths in an infinite homogeneous medium produced by a slab source irradiator is presented. The product minimum dose rate for a set of target thickness is obtained. The design method steps are detailed and an illustrative example is presented. (author)

A distributed parameter model of transmission line transformer for high voltage nanosecond pulse generation.

Science.gov (United States)

Li, Jiangtao; Zhao, Zheng; Li, Longjie; He, Jiaxin; Li, Chenjie; Wang, Yifeng; Su, Can

2017-09-01

A transmission line transformer has potential advantages for nanosecond pulse generation including excellent frequency response and no leakage inductance. The wave propagation process in a secondary mode line is indispensable due to an obvious inside transient electromagnetic transition in this scenario. The equivalent model of the transmission line transformer is crucial for predicting the output waveform and evaluating the effects of magnetic cores on output performance. However, traditional lumped parameter models are not sufficient for nanosecond pulse generation due to the natural neglect of wave propagations in secondary mode lines based on a lumped parameter assumption. In this paper, a distributed parameter model of transmission line transformer was established to investigate wave propagation in the secondary mode line and its influential factors through theoretical analysis and experimental verification. The wave propagation discontinuity in the secondary mode line induced by magnetic cores is emphasized. Characteristics of the magnetic core under a nanosecond pulse were obtained by experiments. Distribution and formation of the secondary mode current were determined for revealing essential wave propagation processes in secondary mode lines. The output waveform and efficiency were found to be affected dramatically by wave propagation discontinuity in secondary mode lines induced by magnetic cores. The proposed distributed parameter model was proved more suitable for nanosecond pulse generation in aspects of secondary mode current, output efficiency, and output waveform. In depth, comprehension of underlying mechanisms and a broader view of the working principle of the transmission line transformer for nanosecond pulse generation can be obtained through this research.

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

Grain Size and Parameter Recovery with TIMSS and the General Diagnostic Model

Science.gov (United States)

Skaggs, Gary; Wilkins, Jesse L. M.; Hein, Serge F.

2016-01-01

The purpose of this study was to explore the degree of grain size of the attributes and the sample sizes that can support accurate parameter recovery with the General Diagnostic Model (GDM) for a large-scale international assessment. In this resampling study, bootstrap samples were obtained from the 2003 Grade 8 TIMSS in Mathematics at varying…

Estimation of k-ε parameters using surrogate models and jet-in-crossflow data

Energy Technology Data Exchange (ETDEWEB)

Lefantzi, Sophia [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Ray, Jaideep [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Arunajatesan, Srinivasan [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Dechant, Lawrence [Sandia National Lab. (SNL-CA), Livermore, CA (United States)

2014-11-01

We demonstrate a Bayesian method that can be used to calibrate computationally expensive 3D RANS (Reynolds Av- eraged Navier Stokes) models with complex response surfaces. Such calibrations, conditioned on experimental data, can yield turbulence model parameters as probability density functions (PDF), concisely capturing the uncertainty in the parameter estimates. Methods such as Markov chain Monte Carlo (MCMC) estimate the PDF by sampling, with each sample requiring a run of the RANS model. Consequently a quick-running surrogate is used instead to the RANS simulator. The surrogate can be very difficult to design if the model's response i.e., the dependence of the calibration variable (the observable) on the parameter being estimated is complex. We show how the training data used to construct the surrogate can be employed to isolate a promising and physically realistic part of the parameter space, within which the response is well-behaved and easily modeled. We design a classifier, based on treed linear models, to model the "well-behaved region". This classifier serves as a prior in a Bayesian calibration study aimed at estimating 3 k - ε parameters ( C _μ, C _ε2 , C _ε1 ) from experimental data of a transonic jet-in-crossflow interaction. The robustness of the calibration is investigated by checking its predictions of variables not included in the cal- ibration data. We also check the limit of applicability of the calibration by testing at off-calibration flow regimes. We find that calibration yield turbulence model parameters which predict the flowfield far better than when the nomi- nal values of the parameters are used. Substantial improvements are still obtained when we use the calibrated RANS model to predict jet-in-crossflow at Mach numbers and jet strengths quite different from those used to generate the ex- perimental (calibration) data. Thus the primary reason for poor predictive skill of RANS, when using nominal

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.

Bayesian Inference on the Memory Parameter for Gamma-Modulated Regression Models

Directory of Open Access Journals (Sweden)

Plinio Andrade

2015-09-01

Full Text Available In this work, we propose a Bayesian methodology to make inferences for the memory parameter and other characteristics under non-standard assumptions for a class of stochastic processes. This class generalizes the Gamma-modulated process, with trajectories that exhibit long memory behavior, as well as decreasing variability as time increases. Different values of the memory parameter influence the speed of this decrease, making this heteroscedastic model very flexible. Its properties are used to implement an approximate Bayesian computation and MCMC scheme to obtain posterior estimates. We test and validate our method through simulations and real data from the big earthquake that occurred in 2010 in Chile.

Evaluation of Hydraulic Parameters Obtained by Different Measurement Methods for Heterogeneous Gravel Soil

Directory of Open Access Journals (Sweden)

Chen Zeng

2012-01-01

Full Text Available Knowledge of soil hydraulic parameters for the van Genuchten function is important to characterize soil water movement for watershed management. Accurate and rapid prediction of soil water flow in heterogeneous gravel soil has become a hot topic in recent years. However, it is difficult to precisely estimate hydraulic parameters in a heterogeneous soil with rock fragments. In this study, the HYDRUS-2D numerical model was used to evaluate hydraulic parameters for heterogeneous gravel soil that was irregularly embedded with rock fragments in a grape production base. The centrifugal method (CM, tensiometer method (TM and inverse solution method (ISM were compared for various parameters in the van Genuchten function. The soil core method (SCM, disc infiltration method (DIM and inverse solution method (ISM were also investigated for measuring saturated hydraulic conductivity. Simulation with the DIM approach revealed a problem of overestimating soil water infiltration whereas simulation with the SCM approach revealed a problem of underestimating water movement as compared to actual field observation. The ISM approach produced the best simulation result even though this approach slightly overestimated soil moisture by ignoring the impact of rock fragments. This study provides useful information on the overall evaluation of soil hydraulic parameters attained with different measurement methods for simulating soil water movement and distribution in heterogeneous gravel soil.

Positioning performance of the NTCM model driven by GPS Klobuchar model parameters

Science.gov (United States)

Hoque, Mohammed Mainul; Jakowski, Norbert; Berdermann, Jens

2018-03-01

Users of the Global Positioning System (GPS) utilize the Ionospheric Correction Algorithm (ICA) also known as Klobuchar model for correcting ionospheric signal delay or range error. Recently, we developed an ionosphere correction algorithm called NTCM-Klobpar model for single frequency GNSS applications. The model is driven by a parameter computed from GPS Klobuchar model and consecutively can be used instead of the GPS Klobuchar model for ionospheric corrections. In the presented work we compare the positioning solutions obtained using NTCM-Klobpar with those using the Klobuchar model. Our investigation using worldwide ground GPS data from a quiet and a perturbed ionospheric and geomagnetic activity period of 17 days each shows that the 24-hour prediction performance of the NTCM-Klobpar is better than the GPS Klobuchar model in global average. The root mean squared deviation of the 3D position errors are found to be about 0.24 and 0.45 m less for the NTCM-Klobpar compared to the GPS Klobuchar model during quiet and perturbed condition, respectively. The presented algorithm has the potential to continuously improve the accuracy of GPS single frequency mass market devices with only little software modification.

Obtaining of Analytical Relations for Hydraulic Parameters of Channels With Two Phase Flow Using Open CFD Toolbox

Science.gov (United States)

Varseev, E.

2017-11-01

The present work is dedicated to verification of numerical model in standard solver of open-source CFD code OpenFOAM for two-phase flow simulation and to determination of so-called “baseline” model parameters. Investigation of heterogeneous coolant flow parameters, which leads to abnormal friction increase of channel in two-phase adiabatic “water-gas” flows with low void fractions, presented.

An iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions, Addendum

Science.gov (United States)

Peters, B. C., Jr.; Walker, H. F.

1975-01-01

New results and insights concerning a previously published iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions were discussed. It was shown that the procedure converges locally to the consistent maximum likelihood estimate as long as a specified parameter is bounded between two limits. Bound values were given to yield optimal local convergence.

A parameters optimization method for planar joint clearance model and its application for dynamics simulation of reciprocating compressor

Science.gov (United States)

Hai-yang, Zhao; Min-qiang, Xu; Jin-dong, Wang; Yong-bo, Li

2015-05-01

In order to improve the accuracy of dynamics response simulation for mechanism with joint clearance, a parameter optimization method for planar joint clearance contact force model was presented in this paper, and the optimized parameters were applied to the dynamics response simulation for mechanism with oversized joint clearance fault. By studying the effect of increased clearance on the parameters of joint clearance contact force model, the relation of model parameters between different clearances was concluded. Then the dynamic equation of a two-stage reciprocating compressor with four joint clearances was developed using Lagrange method, and a multi-body dynamic model built in ADAMS software was used to solve this equation. To obtain a simulated dynamic response much closer to that of experimental tests, the parameters of joint clearance model, instead of using the designed values, were optimized by genetic algorithms approach. Finally, the optimized parameters were applied to simulate the dynamics response of model with oversized joint clearance fault according to the concluded parameter relation. The dynamics response of experimental test verified the effectiveness of this application.

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.

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

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...

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)

SDSS-II: Determination of shape and color parameter coefficients for SALT-II fit model

Energy Technology Data Exchange (ETDEWEB)

Dojcsak, L.; Marriner, J.; /Fermilab

2010-08-01

In this study we look at the SALT-II model of Type IA supernova analysis, which determines the distance moduli based on the known absolute standard candle magnitude of the Type IA supernovae. We take a look at the determination of the shape and color parameter coefficients, {alpha} and {beta} respectively, in the SALT-II model with the intrinsic error that is determined from the data. Using the SNANA software package provided for the analysis of Type IA supernovae, we use a standard Monte Carlo simulation to generate data with known parameters to use as a tool for analyzing the trends in the model based on certain assumptions about the intrinsic error. In order to find the best standard candle model, we try to minimize the residuals on the Hubble diagram by calculating the correct shape and color parameter coefficients. We can estimate the magnitude of the intrinsic errors required to obtain results with {chi}{sup 2}/degree of freedom = 1. We can use the simulation to estimate the amount of color smearing as indicated by the data for our model. We find that the color smearing model works as a general estimate of the color smearing, and that we are able to use the RMS distribution in the variables as one method of estimating the correct intrinsic errors needed by the data to obtain the correct results for {alpha} and {beta}. We then apply the resultant intrinsic error matrix to the real data and show our results.

Ann modeling of kerf transfer in Co2 laser cutting and optimization of cutting parameters using monte carlo method

Directory of Open Access Journals (Sweden)

Miloš Madić

2015-01-01

Full Text Available In this paper, an attempt has been made to develop a mathematical model in order to study the relationship between laser cutting parameters such as laser power, cutting speed, assist gas pressure and focus position, and kerf taper angle obtained in CO2 laser cutting of AISI 304 stainless steel. To this aim, a single hidden layer artificial neural network (ANN trained with gradient descent with momentum algorithm was used. To obtain an experimental database for the ANN training, laser cutting experiment was planned as per Taguchi’s L27 orthogonal array with three levels for each of the cutting parameters. Statistically assessed as adequate, ANN model was then used to investigate the effect of the laser cutting parameters on the kerf taper angle by generating 2D and 3D plots. It was observed that the kerf taper angle was highly sensitive to the selected laser cutting parameters, as well as their interactions. In addition to modeling, by applying the Monte Carlo method on the developed kerf taper angle ANN model, the near optimal laser cutting parameter settings, which minimize kerf taper angle, were determined.

Artificial Neural Network model for the determination of GSM Rxlevel from atmospheric parameters

Directory of Open Access Journals (Sweden)

Julia Ofure Eichie

2017-04-01

Full Text Available Accurate received signal level (Rxlevel values are useful for mobile telecommunication network planning. Rxlevel is affected by the dynamics of the atmosphere through which it propagates. Adequate knowledge of the prevailing atmospheric conditions in an environment is essential for proper network planning. However most of the existing GSM received signal determination model are function of distance between point of signal reception and transmitting site thus necessitating the development of a model that involve the use of atmospheric parameters in the determination of received GSM signal level. In this paper, a three stage approach was used in the development of the model using some atmospheric parameters such as atmospheric temperature, relative humidity and dew point. The selected and easily measurable atmospheric parameters were used as input parameters in developing two new models for computing the Rxlevel of GSM signal using a three-step approach. Data acquisition and pre-processing serves as the first stage and formulation of ANN design and the development of parametric model for the Rxlevel using ANN synaptic weights form the second stage of the proposed approach. The third stage involves the use of ANN weight and bias values, and network architecture in the development of the model equation. In evaluating the performance of the proposed models, network parameters were varied and the results obtained using mean squared error (MSE as performance measure showed the developed model with 33 neurons in the hidden layer and tansig activation, function in both the hidden and output layers as the optimal model with least MSE value of 0.056. Thus showing that the developed model has an acceptable accuracy value as demonstrated from comparison of results with actual measured values.

Inverse analyses of effective diffusion parameters relevant for a two-phase moisture model of cementitious materials

DEFF Research Database (Denmark)

Addassi, Mouadh; Johannesson, Björn; Wadsö, Lars

2018-01-01

Here we present an inverse analyses approach to determining the two-phase moisture transport properties relevant to concrete durability modeling. The purposed moisture transport model was based on a continuum approach with two truly separate equations for the liquid and gas phase being connected...... test, and, (iv) capillary suction test. Mass change over time, as obtained from the drying test, the two different cup test intervals and the capillary suction test, was used to obtain the effective diffusion parameters using the proposed inverse analyses approach. The moisture properties obtained...

SiC-VJFETs power switching devices: an improved model and parameter optimization technique

Science.gov (United States)

Ben Salah, T.; Lahbib, Y.; Morel, H.

2009-12-01

Silicon carbide junction field effect transistor (SiC-JFETs) is a mature power switch newly applied in several industrial applications. SiC-JFETs are often simulated by Spice model in order to predict their electrical behaviour. Although such a model provides sufficient accuracy for some applications, this paper shows that it presents serious shortcomings in terms of the neglect of the body diode model, among many others in circuit model topology. Simulation correction is then mandatory and a new model should be proposed. Moreover, this paper gives an enhanced model based on experimental dc and ac data. New devices are added to the conventional circuit model giving accurate static and dynamic behaviour, an effect not accounted in the Spice model. The improved model is implemented into VHDL-AMS language and steady-state dynamic and transient responses are simulated for many SiC-VJFETs samples. Very simple and reliable optimization algorithm based on the optimization of a cost function is proposed to extract the JFET model parameters. The obtained parameters are verified by comparing errors between simulations results and experimental data.

A high-throughput 2D-analytical technique to obtain single protein parameters from complex cell lysates for in silico process development of ion exchange chromatography.

Science.gov (United States)

Kröner, Frieder; Elsäßer, Dennis; Hubbuch, Jürgen

2013-11-29

The accelerating growth of the market for biopharmaceutical proteins, the market entry of biosimilars and the growing interest in new, more complex molecules constantly pose new challenges for bioseparation process development. In the presented work we demonstrate the application of a multidimensional, analytical separation approach to obtain the relevant physicochemical parameters of single proteins in a complex mixture for in silico chromatographic process development. A complete cell lysate containing a low titre target protein was first fractionated by multiple linear salt gradient anion exchange chromatography (AEC) with varying gradient length. The collected fractions were subsequently analysed by high-throughput capillary gel electrophoresis (HT-CGE) after being desalted and concentrated. From the obtained data of the 2D-separation the retention-volumes and the concentration of the single proteins were determined. The retention-volumes of the single proteins were used to calculate the related steric-mass action model parameters. In a final evaluation experiment the received parameters were successfully applied to predict the retention behaviour of the single proteins in salt gradient AEC. Copyright © 2013 Elsevier B.V. All rights reserved.

PARAMETER COORDINATION AND ROBUST OPTIMIZATION FOR MULTIDISCIPLINARY DESIGN

Institute of Scientific and Technical Information of China (English)

HU Jie; PENG Yinghong; XIONG Guangleng

2006-01-01

A new parameter coordination and robust optimization approach for multidisciplinary design is presented. Firstly, the constraints network model is established to support engineering change, coordination and optimization. In this model, interval boxes are adopted to describe the uncertainty of design parameters quantitatively to enhance the design robustness. Secondly, the parameter coordination method is presented to solve the constraints network model, monitor the potential conflicts due to engineering changes, and obtain the consistency solution space corresponding to the given product specifications. Finally, the robust parameter optimization model is established, and genetic arithmetic is used to obtain the robust optimization parameter. An example of bogie design is analyzed to show the scheme to be effective.

Customized Steady-State Constraints for Parameter Estimation in Non-Linear Ordinary Differential Equation Models.

Science.gov (United States)

Rosenblatt, Marcus; Timmer, Jens; Kaschek, Daniel

2016-01-01

Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g., by maximum-likelihood estimation. In particular, models of biological systems contain a large number of parameters. To reduce the dimensionality of the parameter space, steady-state information is incorporated in the parameter estimation process. For non-linear models, analytical steady-state calculation typically leads to higher-order polynomial equations for which no closed-form solutions can be obtained. This can be circumvented by solving the steady-state equations for kinetic parameters, which results in a linear equation system with comparatively simple solutions. At the same time multiplicity of steady-state solutions is avoided, which otherwise is problematic for optimization. When solved for kinetic parameters, however, steady-state constraints tend to become negative for particular model specifications, thus, generating new types of optimization problems. Here, we present an algorithm based on graph theory that derives non-negative, analytical steady-state expressions by stepwise removal of cyclic dependencies between dynamical variables. The algorithm avoids multiple steady-state solutions by construction. We show that our method is applicable to most common classes of biochemical reaction networks containing inhibition terms, mass-action and Hill-type kinetic equations. Comparing the performance of parameter estimation for different analytical and numerical methods of incorporating steady-state information, we show that our approach is especially well-tailored to guarantee a high success rate of optimization.

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 numerical method to estimate the parameters of the CEV model implied by American option prices: Evidence from NYSE

International Nuclear Information System (INIS)

Ballestra, Luca Vincenzo; Cecere, Liliana

2016-01-01

Highlights: • We develop a method to compute the parameters of the CEV model implied by American options. • This is the first procedure for calibrating the CEV model to American option prices. • The proposed approach is extensively tested on the NYSE market. • The novel method turns out to be very efficient in computing the CEV model parameters. • The CEV model provides only a marginal improvement over the lognormal model. - Abstract: We develop a highly efficient procedure to forecast the parameters of the constant elasticity of variance (CEV) model implied by American options. In particular, first of all, the American option prices predicted by the CEV model are calculated using an accurate and fast finite difference scheme. Then, the parameters of the CEV model are obtained by minimizing the distance between theoretical and empirical option prices, which yields an optimization problem that is solved using an ad-hoc numerical procedure. The proposed approach, which turns out to be very efficient from the computational standpoint, is used to test the goodness-of-fit of the CEV model in predicting the prices of American options traded on the NYSE. The results obtained reveal that the CEV model does not provide a very good agreement with real market data and yields only a marginal improvement over the more popular Black–Scholes model.

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.

Parameter Estimation of a Delay Time Model of Wearing Parts Based on Objective Data

Directory of Open Access Journals (Sweden)

Y. Tang

2015-01-01

Full Text Available The wearing parts of a system have a very high failure frequency, making it necessary to carry out continual functional inspections and maintenance to protect the system from unscheduled downtime. This allows for the collection of a large amount of maintenance data. Taking the unique characteristics of the wearing parts into consideration, we establish their respective delay time models in ideal inspection cases and nonideal inspection cases. The model parameters are estimated entirely using the collected maintenance data. Then, a likelihood function of all renewal events is derived based on their occurring probability functions, and the model parameters are calculated with the maximum likelihood function method, which is solved by the CRM. Finally, using two wearing parts from the oil and gas drilling industry as examples—the filter element and the blowout preventer rubber core—the parameters of the distribution function of the initial failure time and the delay time for each example are estimated, and their distribution functions are obtained. Such parameter estimation based on objective data will contribute to the optimization of the reasonable function inspection interval and will also provide some theoretical models to support the integrity management of equipment or systems.

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.

A Central Composite Face-Centered Design for Parameters Estimation of PEM Fuel Cell Electrochemical Model

Directory of Open Access Journals (Sweden)

Khaled MAMMAR

2013-11-01

Full Text Available In this paper, a new approach based on Experimental of design methodology (DoE is used to estimate the optimal of unknown model parameters proton exchange membrane fuel cell (PEMFC. This proposed approach combines the central composite face-centered (CCF and numerical PEMFC electrochemical. Simulation results obtained using electrochemical model help to predict the cell voltage in terms of inlet partial pressures of hydrogen and oxygen, stack temperature, and operating current. The value of the previous model and (CCF design methodology is used for parametric analysis of electrochemical model. Thus it is possible to evaluate the relative importance of each parameter to the simulation accuracy. However this methodology is able to define the exact values of the parameters from the manufacture data. It was tested for the BCS 500-W stack PEM Generator, a stack rated at 500 W, manufactured by American Company BCS Technologies FC.

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...

Parameter sensitivity analysis of a 1-D cold region lake model for land-surface schemes

Science.gov (United States)

Guerrero, José-Luis; Pernica, Patricia; Wheater, Howard; Mackay, Murray; Spence, Chris

2017-12-01

Lakes might be sentinels of climate change, but the uncertainty in their main feedback to the atmosphere - heat-exchange fluxes - is often not considered within climate models. Additionally, these fluxes are seldom measured, hindering critical evaluation of model output. Analysis of the Canadian Small Lake Model (CSLM), a one-dimensional integral lake model, was performed to assess its ability to reproduce diurnal and seasonal variations in heat fluxes and the sensitivity of simulated fluxes to changes in model parameters, i.e., turbulent transport parameters and the light extinction coefficient (Kd). A C++ open-source software package, Problem Solving environment for Uncertainty Analysis and Design Exploration (PSUADE), was used to perform sensitivity analysis (SA) and identify the parameters that dominate model behavior. The generalized likelihood uncertainty estimation (GLUE) was applied to quantify the fluxes' uncertainty, comparing daily-averaged eddy-covariance observations to the output of CSLM. Seven qualitative and two quantitative SA methods were tested, and the posterior likelihoods of the modeled parameters, obtained from the GLUE analysis, were used to determine the dominant parameters and the uncertainty in the modeled fluxes. Despite the ubiquity of the equifinality issue - different parameter-value combinations yielding equivalent results - the answer to the question was unequivocal: Kd, a measure of how much light penetrates the lake, dominates sensible and latent heat fluxes, and the uncertainty in their estimates is strongly related to the accuracy with which Kd is determined. This is important since accurate and continuous measurements of Kd could reduce modeling uncertainty.

How Sensitive Are Transdermal Transport Predictions by Microscopic Stratum Corneum Models to Geometric and Transport Parameter Input?

Science.gov (United States)

Wen, Jessica; Koo, Soh Myoung; Lape, Nancy

2018-02-01

While predictive models of transdermal transport have the potential to reduce human and animal testing, microscopic stratum corneum (SC) model output is highly dependent on idealized SC geometry, transport pathway (transcellular vs. intercellular), and penetrant transport parameters (e.g., compound diffusivity in lipids). Most microscopic models are limited to a simple rectangular brick-and-mortar SC geometry and do not account for variability across delivery sites, hydration levels, and populations. In addition, these models rely on transport parameters obtained from pure theory, parameter fitting to match in vivo experiments, and time-intensive diffusion experiments for each compound. In this work, we develop a microscopic finite element model that allows us to probe model sensitivity to variations in geometry, transport pathway, and hydration level. Given the dearth of experimentally-validated transport data and the wide range in theoretically-predicted transport parameters, we examine the model's response to a variety of transport parameters reported in the literature. Results show that model predictions are strongly dependent on all aforementioned variations, resulting in order-of-magnitude differences in lag times and permeabilities for distinct structure, hydration, and parameter combinations. This work demonstrates that universally predictive models cannot fully succeed without employing experimentally verified transport parameters and individualized SC structures. Copyright © 2018 American Pharmacists Association®. Published by Elsevier Inc. All rights reserved.

Assessing seasonality of biochemical CO2 exchange model parameters from micrometeorological flux observations at boreal coniferous forest

Directory of Open Access Journals (Sweden)

T. Vesala

2008-12-01

Full Text Available The seasonality of the NEE of the northern boreal coniferous forests was investigated by means of inversion modelling using eddy covariance data. Eddy covariance data was used to optimize the biochemical model parameters. Our study sites consisted of three Scots pine (l. Pinus sylvestris forests and one Norway spruce (l. Picea abies forest that were located in Finland and Sweden. We obtained temperature and seasonal dependence for the biochemical model parameters: the maximum rate of carboxylation (Vc(max and the maximum rate of electron transport (Jmax. Both of the parameters were optimized without assumptions about their mutual magnitude. The values obtained for the biochemical model parameters were similar at all the sites during summer time. To describe seasonality, different temperature fits were made for the spring, summer and autumn periods. During summer, average Jmax across the sites was 54.0 μmol m−2 s−1 (variance 31.2 μmol m−2 s−1 and Vc(max was 12.0 μmol m−2 s−1 (variance 6.6 μmol m−2 s−1 at 17°C. The sensitivity of the model to LAI and atmospheric soil water stress was also studied. The impact of seasonality on annual GPP was 17% when only summertime parameterization was used throughout the year compared to seasonally changing parameterizations.

Performance assessment model development and parameter acquisition for analysis of the transport of natural radionuclides in a Mediterranean watershed

International Nuclear Information System (INIS)

Agueero, Almudena

2005-01-01

This paper describes the methodology developed to construct a model for predicting the behaviour of the natural radioisotopes of U, Th and Ra in a Mediterranean watershed. The methodology includes the development of the performance assessment model, obtaining water flow and radiological parameters based on experimental data and analysis of results. The model, which accounts for both water flows and mass balances of the radionuclides in a semi-natural environment, provides assessments of radionuclide behaviour in grassland and agricultural soils, rivers and reservoirs, including the processes of radionuclide migration through land and water and interactions between both. From field and laboratory data, it has been possible to obtain parameters for the driving processes considered in the model, water fluxes, source term definition, soil to plant transfer factors and distribution coefficient values. Ranges of parameter values obtained have shown good agreement with published literature data. This general methodological approach was developed to be extended to other radionuclides for the modelling of a biosphere watershed in the context of performance assessment of a High Level Waste (HLW) repository under Mediterranean climate conditions, as well as for forecasting radionuclide transport under similar Mediterranean conditions that will occur in the future in other areas. The application of sensitivity and uncertainty analysis was intended to identify key uncertainties with the aim of setting priorities for future research. The model results for the activity concentration in the reservoir indicate that for 238 U and 230 Th the most relevant parameter is the initial concentrations of the radionuclides in the reservoir sediments. However, for 226 Ra the most important parameter is the precipitation rate over the whole watershed

Significance of settling model structures and parameter subsets in modelling WWTPs under wet-weather flow and filamentous bulking conditions.

Science.gov (United States)

Ramin, Elham; Sin, Gürkan; Mikkelsen, Peter Steen; Plósz, Benedek Gy

2014-10-15

Current research focuses on predicting and mitigating the impacts of high hydraulic loadings on centralized wastewater treatment plants (WWTPs) under wet-weather conditions. The maximum permissible inflow to WWTPs depends not only on the settleability of activated sludge in secondary settling tanks (SSTs) but also on the hydraulic behaviour of SSTs. The present study investigates the impacts of ideal and non-ideal flow (dry and wet weather) and settling (good settling and bulking) boundary conditions on the sensitivity of WWTP model outputs to uncertainties intrinsic to the one-dimensional (1-D) SST model structures and parameters. We identify the critical sources of uncertainty in WWTP models through global sensitivity analysis (GSA) using the Benchmark simulation model No. 1 in combination with first- and second-order 1-D SST models. The results obtained illustrate that the contribution of settling parameters to the total variance of the key WWTP process outputs significantly depends on the influent flow and settling conditions. The magnitude of the impact is found to vary, depending on which type of 1-D SST model is used. Therefore, we identify and recommend potential parameter subsets for WWTP model calibration, and propose optimal choice of 1-D SST models under different flow and settling boundary conditions. Additionally, the hydraulic parameters in the second-order SST model are found significant under dynamic wet-weather flow conditions. These results highlight the importance of developing a more mechanistic based flow-dependent hydraulic sub-model in second-order 1-D SST models in the future. Copyright © 2014 Elsevier Ltd. All rights reserved.

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

A method to obtain the thermal parameters and the photothermal transduction efficiency in an optical hyperthermia device based on laser irradiation of gold nanoparticles.

Science.gov (United States)

Sánchez López de Pablo, Cristina; Olmedo, José Javier Serrano; Rosales, Alejandra Mina; Ramírez Hernández, Norma; Del Pozo Guerrero, Francisco

2014-01-01

Optical hyperthermia systems based on the laser irradiation of gold nanorods seem to be a promising tool in the development of therapies against cancer. After a proof of concept in which the authors demonstrated the efficiency of this kind of systems, a modeling process based on an equivalent thermal-electric circuit has been carried out to determine the thermal parameters of the system and an energy balance obtained from the time-dependent heating and cooling temperature curves of the irradiated samples in order to obtain the photothermal transduction efficiency. By knowing this parameter, it is possible to increase the effectiveness of the treatments, thanks to the possibility of predicting the response of the device depending on the working configuration. As an example, the thermal behavior of two different kinds of nanoparticles is compared. The results show that, under identical conditions, the use of PEGylated gold nanorods allows for a more efficient heating compared with bare nanorods, and therefore, it results in a more effective therapy.

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

Effect of Uncertainty Parameters in Blowdown and Reflood Models for OPR1000 LBLOCA Analysis

Energy Technology Data Exchange (ETDEWEB)

Huh, Byung Gil; Jin, Chang Yong; Seul, Kwangwon; Hwang, Taesuk [Korea Institute of Nuclear Safety, Daejeon (Korea, Republic of)

2014-05-15

KINS(Korea Institute of Nuclear Safety) has also performed the audit calculation with the KINS Realistic Evaluation Methodology(KINS-REM) to confirm the validity of licensee's calculation. In the BEPU method, it is very important to quantify the code and model uncertainty. It is referred in the following requirement: BE calculations in Regulatory Guide 1.157 - 'the code and models used are acceptable and applicable to the specific facility over the intended operating range and must quantify the uncertainty in the specific application'. In general, the uncertainty of model/code should be obtained through the data comparison with relevant integral- and separate-effect tests at different scales. However, it is not easy to determine these kinds of uncertainty because of the difficulty for evaluating accurately various experiments. Therefore, the expert judgment has been used in many cases even with the limitation that the uncertainty range of important parameters can be wide and inaccurate. In the KINS-REM, six heat transfer parameters in the blowdown phase have been used to consider the uncertainty of models. Recently, MARS-KS code was modified to consider the uncertainty of the five heat transfer parameters in the reflood phase. Accordingly, it is required that the uncertainty range for parameters of reflood models is determined and the effect of these ranges is evaluated. In this study, the large break LOCA (LBLOCA) analysis for OPR1000 was performed to identify the effect of uncertainty parameters in blowdown and reflood models.

Rotating Parabolic-Reflector Antenna Target in SAR Data: Model, Characteristics, and Parameter Estimation

Directory of Open Access Journals (Sweden)

Bin Deng

2013-01-01

Full Text Available Parabolic-reflector antennas (PRAs, usually possessing rotation, are a particular type of targets of potential interest to the synthetic aperture radar (SAR community. This paper is aimed to investigate PRA’s scattering characteristics and then to extract PRA’s parameters from SAR returns, for supporting image interpretation and target recognition. We at first obtain both closed-form and numeric solutions to PRA’s backscattering by geometrical optics (GO, physical optics, and graphical electromagnetic computation, respectively. Based on the GO solution, a migratory scattering center model is at first presented for representing the movement of the specular point with aspect angle, and then a hybrid model, named the migratory/micromotion scattering center (MMSC model, is proposed for characterizing a rotating PRA in the SAR geometry, which incorporates PRA’s rotation into its migratory scattering center model. Additionally, we in detail analyze PRA’s radar characteristics on radar cross-section, high-resolution range profiles, time-frequency distribution, and 2D images, which also confirm the models proposed. A maximal likelihood estimator is developed for jointly solving the MMSC model for PRA’s multiple parameters by optimization. By exploiting the aforementioned characteristics, the coarse parameter estimation guarantees convergency upon global minima. The signatures recovered can be favorably utilized for SAR image interpretation and target recognition.

Multivariate analysis and extraction of parameters in resistive RAMs using the Quantum Point Contact model

Science.gov (United States)

Roldán, J. B.; Miranda, E.; González-Cordero, G.; García-Fernández, P.; Romero-Zaliz, R.; González-Rodelas, P.; Aguilera, A. M.; González, M. B.; Jiménez-Molinos, F.

2018-01-01

A multivariate analysis of the parameters that characterize the reset process in Resistive Random Access Memory (RRAM) has been performed. The different correlations obtained can help to shed light on the current components that contribute in the Low Resistance State (LRS) of the technology considered. In addition, a screening method for the Quantum Point Contact (QPC) current component is presented. For this purpose, the second derivative of the current has been obtained using a novel numerical method which allows determining the QPC model parameters. Once the procedure is completed, a whole Resistive Switching (RS) series of thousands of curves is studied by means of a genetic algorithm. The extracted QPC parameter distributions are characterized in depth to get information about the filamentary pathways associated with LRS in the low voltage conduction regime.

A Novel Non-Iterative Method for Real-Time Parameter Estimation of the Fricke-Morse Model

Directory of Open Access Journals (Sweden)

SIMIC, M.

2016-11-01

Full Text Available Parameter estimation of Fricke-Morse model of biological tissue is widely used in bioimpedance data processing and analysis. Complex nonlinear least squares (CNLS data fitting is often used for parameter estimation of the model, but limitations such as high processing time, converging into local minimums, need for good initial guess of model parameters and non-convergence have been reported. Thus, there is strong motivation to develop methods which can solve these flaws. In this paper a novel real-time method for parameter estimation of Fricke-Morse model of biological cells is presented. The proposed method uses the value of characteristic frequency estimated from the measured imaginary part of bioimpedance, whereupon the Fricke-Morse model parameters are calculated using the provided analytical expressions. The proposed method is compared with CNLS in frequency ranges of 1 kHz to 10 MHz (beta-dispersion and 10 kHz to 100 kHz, which is more suitable for low-cost microcontroller-based bioimpedance measurement systems. The obtained results are promising, and in both frequency ranges, CNLS and the proposed method have accuracies suitable for most electrical bioimpedance (EBI applications. However, the proposed algorithm has significantly lower computation complexity, so it was 20-80 times faster than CNLS.

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

Modeling Coronal Mass Ejections with EUHFORIA: A Parameter Study of the Gibson-Low Flux Rope Model using Multi-Viewpoint Observations

Science.gov (United States)

Verbeke, C.; Asvestari, E.; Scolini, C.; Pomoell, J.; Poedts, S.; Kilpua, E.

2017-12-01

Coronal Mass Ejections (CMEs) are one of the big influencers on the coronal and interplanetary dynamics. Understanding their origin and evolution from the Sun to the Earth is crucial in order to determine the impact on our Earth and society. One of the key parameters that determine the geo-effectiveness of the coronal mass ejection is its internal magnetic configuration. We present a detailed parameter study of the Gibson-Low flux rope model. We focus on changes in the input parameters and how these changes affect the characteristics of the CME at Earth. Recently, the Gibson-Low flux rope model has been implemented into the inner heliosphere model EUHFORIA, a magnetohydrodynamics forecasting model of large-scale dynamics from 0.1 AU up to 2 AU. Coronagraph observations can be used to constrain the kinematics and morphology of the flux rope. One of the key parameters, the magnetic field, is difficult to determine directly from observations. In this work, we approach the problem by conducting a parameter study in which flux ropes with varying magnetic configurations are simulated. We then use the obtained dataset to look for signatures in imaging observations and in-situ observations in order to find an empirical way of constraining the parameters related to the magnetic field of the flux rope. In particular, we focus on events observed by at least two spacecraft (STEREO + L1) in order to discuss the merits of using observations from multiple viewpoints in constraining the parameters.

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.

A simple method for identifying parameter correlations in partially observed linear dynamic models.

Science.gov (United States)

Li, Pu; Vu, Quoc Dong

2015-12-14

Parameter estimation represents one of the most significant challenges in systems biology. This is because biological models commonly contain a large number of parameters among which there may be functional interrelationships, thus leading to the problem of non-identifiability. Although identifiability analysis has been extensively studied by analytical as well as numerical approaches, systematic methods for remedying practically non-identifiable models have rarely been investigated. We propose a simple method for identifying pairwise correlations and higher order interrelationships of parameters in partially observed linear dynamic models. This is made by derivation of the output sensitivity matrix and analysis of the linear dependencies of its columns. Consequently, analytical relations between the identifiability of the model parameters and the initial conditions as well as the input functions can be achieved. In the case of structural non-identifiability, identifiable combinations can be obtained by solving the resulting homogenous linear equations. In the case of practical non-identifiability, experiment conditions (i.e. initial condition and constant control signals) can be provided which are necessary for remedying the non-identifiability and unique parameter estimation. It is noted that the approach does not consider noisy data. In this way, the practical non-identifiability issue, which is popular for linear biological models, can be remedied. Several linear compartment models including an insulin receptor dynamics model are taken to illustrate the application of the proposed approach. Both structural and practical identifiability of partially observed linear dynamic models can be clarified by the proposed method. The result of this method provides important information for experimental design to remedy the practical non-identifiability if applicable. The derivation of the method is straightforward and thus the algorithm can be easily implemented into a

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

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

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.

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.

Genetic Algorithms for a Parameter Estimation of a Fermentation Process Model: A Comparison

Directory of Open Access Journals (Sweden)

Olympia Roeva

2005-12-01

Full Text Available In this paper the problem of a parameter estimation using genetic algorithms is examined. A case study considering the estimation of 6 parameters of a nonlinear dynamic model of E. coli fermentation is presented as a test problem. The parameter estimation problem is stated as a nonlinear programming problem subject to nonlinear differential-algebraic constraints. This problem is known to be frequently ill-conditioned and multimodal. Thus, traditional (gradient-based local optimization methods fail to arrive satisfied solutions. To overcome their limitations, the use of different genetic algorithms as stochastic global optimization methods is explored. These algorithms are proved to be very suitable for the optimization of highly non-linear problems with many variables. Genetic algorithms can guarantee global optimality and robustness. These facts make them advantageous in use for parameter identification of fermentation models. A comparison between simple, modified and multi-population genetic algorithms is presented. The best result is obtained using the modified genetic algorithm. The considered algorithms converged very closely to the cost value but the modified algorithm is in times faster than other two.

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.

Model Parameter Variability for Enhanced Anaerobic Bioremediation of DNAPL Source Zones

Science.gov (United States)

Mao, X.; Gerhard, J. I.; Barry, D. A.

2005-12-01

modelling parameters typically employed for simulating TCE dechlorination relevant for a range of system conditions (e.g, bioaugmented, high TCE concentrations, etc.). The significance of the obtained variability of parameters is illustrated with one-dimensional simulations of enhanced anaerobic bioremediation of residual TCE DNAPL.

Mistral project: identification and parameter adjustment. Theoretical part; Projet Mistral: identification et recalage des modeles. Etude theorique

Energy Technology Data Exchange (ETDEWEB)

Faille, D.; Codrons, B.; Gevers, M.

1996-03-01

This document belongs to the methodological part of the project MISTRAL, which builds a library of power plant models. The model equations are generally obtained from the first principles. The parameters are actually not always easily calculable (at least accurately) from the dimension data. We are therefore investigating the possibility of automatically adjusting the value of those parameters from experimental data. To do that, we must master the optimization algorithms and the techniques that are analyzing the model structure, like the identifiability theory. (authors). 7 refs., 1 fig., 1 append.

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)

Modeling of mouse eye and errors in ocular parameters affecting refractive state

Science.gov (United States)

Bawa, Gurinder

Rodents eye are particularly used to study refractive error state of an eye and development of refractive eye. Genetic organization of rodents is similar to that of humans, which makes them interesting candidates to be researched upon. From rodents family mice models are encouraged over rats because of availability of genetically engineered models. Despite of extensive work that has been performed on mice and rat models, still no one is able to quantify an optical model, due to variability in the reported ocular parameters. In this Dissertation, we have extracted ocular parameters and generated schematics of eye from the raw data from School of Medicine, Detroit. In order to see how the rays would travel through an eye and the defects associated with an eye; ray tracing has been performed using ocular parameters. Finally we have systematically evaluated the contribution of various ocular parameters, such as radii of curvature of ocular surfaces, thicknesses of ocular components, and refractive indices of ocular refractive media, using variational analysis and a computational model of the rodent eye. Variational analysis revealed that variation in all the ocular parameters does affect the refractive status of the eye, but depending upon the magnitude of the impact those parameters are listed as critical or non critical. Variation in the depth of the vitreous chamber, thickness of the lens, radius of the anterior surface of the cornea, radius of the anterior surface of the lens, as well as refractive indices for the lens and vitreous, appears to have the largest impact on the refractive error and thus are categorized as critical ocular parameters. The radii of the posterior surfaces of the cornea and lens have much smaller contributions to the refractive state, while the radii of the anterior and posterior surfaces of the retina have no effect on the refractive error. These data provide the framework for further refinement of the optical models of the rat and mouse

Application of parameters space analysis tools for empirical model validation

Energy Technology Data Exchange (ETDEWEB)

Paloma del Barrio, E. [LEPT-ENSAM UMR 8508, Talence (France); Guyon, G. [Electricite de France, Moret-sur-Loing (France)

2004-01-01

A new methodology for empirical model validation has been proposed in the framework of the Task 22 (Building Energy Analysis Tools) of the International Energy Agency. It involves two main steps: checking model validity and diagnosis. Both steps, as well as the underlying methods, have been presented in the first part of the paper. In this part, they are applied for testing modelling hypothesis in the framework of the thermal analysis of an actual building. Sensitivity analysis tools have been first used to identify the parts of the model that can be really tested on the available data. A preliminary diagnosis is then supplied by principal components analysis. Useful information for model behaviour improvement has been finally obtained by optimisation techniques. This example of application shows how model parameters space analysis is a powerful tool for empirical validation. In particular, diagnosis possibilities are largely increased in comparison with residuals analysis techniques. (author)

An Investigation on the Sensitivity of the Parameters of Urban Flood Model

Science.gov (United States)

M, A. B.; Lohani, B.; Jain, A.

2015-12-01

Global climatic change has triggered weather patterns which lead to heavy and sudden rainfall in different parts of world. The impact of heavy rainfall is severe especially on urban areas in the form of urban flooding. In order to understand the effect of heavy rainfall induced flooding, it is necessary to model the entire flooding scenario more accurately, which is now becoming possible with the availability of high resolution airborne LiDAR data and other real time observations. However, there is not much understanding on the optimal use of these data and on the effect of other parameters on the performance of the flood model. This study aims at developing understanding on these issues. In view of the above discussion, the aim of this study is to (i) understand that how the use of high resolution LiDAR data improves the performance of urban flood model, and (ii) understand the sensitivity of various hydrological parameters on urban flood modelling. In this study, modelling of flooding in urban areas due to heavy rainfall is carried out considering Indian Institute of Technology (IIT) Kanpur, India as the study site. The existing model MIKE FLOOD, which is accepted by Federal Emergency Management Agency (FEMA), is used along with the high resolution airborne LiDAR data. Once the model is setup it is made to run by changing the parameters such as resolution of Digital Surface Model (DSM), manning's roughness, initial losses, catchment description, concentration time, runoff reduction factor. In order to realize this, the results obtained from the model are compared with the field observations. The parametric study carried out in this work demonstrates that the selection of catchment description plays a very important role in urban flood modelling. Results also show the significant impact of resolution of DSM, initial losses and concentration time on urban flood model. This study will help in understanding the effect of various parameters that should be part of a

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.

Forecasting Models for Some Water Quality Parameters of Shatt Al-Hilla River, Iraq

Directory of Open Access Journals (Sweden)

Rafa H. Al-Suhili

2017-07-01

Full Text Available This paper provides Artificial Neural Networks model versions for forecasting the monthly averages of some chemical water quality parameters of Shatt Al-Hilla River, which is located at Hilla City, south of Iraq. The water quality parameters investigated were Sulphate, Magnesium, Calcium, Alkalinity, and Total Hardness. Results indicate that for Sulphate and Calcium high correlation coefficients models were observed to be (0.9 and 0.88, while for Magnesium, Alkalinity and Hardness low correlation coefficients model were observed to be (0.48,0.58, and 0.51 respectively. Serial correlation behavior of these variables indicate at that high lag time correlations sequences are observed for the first two variables and low ones for the last three water quality parameters. A serial correlation coefficient analysis was done and indicates that as the variable exhibited weak lag correlation structure, then a successful ANN forecasting model could not be obtained even if many trials were done to enhance it's performance, such as increasing the number of nodes, the lagged input variables, and/or changing the learning rate and the momentum term values, or the use of different types of activation functions. On the other hand, those variables that have a strong lag correlation structure can easily fit successful ANN forecasting models

Model parameter estimations from residual gravity anomalies due to simple-shaped sources using Differential Evolution Algorithm

Science.gov (United States)

Ekinci, Yunus Levent; Balkaya, Çağlayan; Göktürkler, Gökhan; Turan, Seçil

2016-06-01

An efficient approach to estimate model parameters from residual gravity data based on differential evolution (DE), a stochastic vector-based metaheuristic algorithm, has been presented. We have showed the applicability and effectiveness of this algorithm on both synthetic and field anomalies. According to our knowledge, this is a first attempt of applying DE for the parameter estimations of residual gravity anomalies due to isolated causative sources embedded in the subsurface. The model parameters dealt with here are the amplitude coefficient (A), the depth and exact origin of causative source (zo and xo, respectively) and the shape factors (q and ƞ). The error energy maps generated for some parameter pairs have successfully revealed the nature of the parameter estimation problem under consideration. Noise-free and noisy synthetic single gravity anomalies have been evaluated with success via DE/best/1/bin, which is a widely used strategy in DE. Additionally some complicated gravity anomalies caused by multiple source bodies have been considered, and the results obtained have showed the efficiency of the algorithm. Then using the strategy applied in synthetic examples some field anomalies observed for various mineral explorations such as a chromite deposit (Camaguey district, Cuba), a manganese deposit (Nagpur, India) and a base metal sulphide deposit (Quebec, Canada) have been considered to estimate the model parameters of the ore bodies. Applications have exhibited that the obtained results such as the depths and shapes of the ore bodies are quite consistent with those published in the literature. Uncertainty in the solutions obtained from DE algorithm has been also investigated by Metropolis-Hastings (M-H) sampling algorithm based on simulated annealing without cooling schedule. Based on the resulting histogram reconstructions of both synthetic and field data examples the algorithm has provided reliable parameter estimations being within the sampling limits of

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)

Bayesian estimation of parameters in a regional hydrological model

Directory of Open Access Journals (Sweden)

K. Engeland

2002-01-01

Full Text Available This study evaluates the applicability of the distributed, process-oriented Ecomag model for prediction of daily streamflow in ungauged basins. The Ecomag model is applied as a regional model to nine catchments in the NOPEX area, using Bayesian statistics to estimate the posterior distribution of the model parameters conditioned on the observed streamflow. The distribution is calculated by Markov Chain Monte Carlo (MCMC analysis. The Bayesian method requires formulation of a likelihood function for the parameters and three alternative formulations are used. The first is a subjectively chosen objective function that describes the goodness of fit between the simulated and observed streamflow, as defined in the GLUE framework. The second and third formulations are more statistically correct likelihood models that describe the simulation errors. The full statistical likelihood model describes the simulation errors as an AR(1 process, whereas the simple model excludes the auto-regressive part. The statistical parameters depend on the catchments and the hydrological processes and the statistical and the hydrological parameters are estimated simultaneously. The results show that the simple likelihood model gives the most robust parameter estimates. The simulation error may be explained to a large extent by the catchment characteristics and climatic conditions, so it is possible to transfer knowledge about them to ungauged catchments. The statistical models for the simulation errors indicate that structural errors in the model are more important than parameter uncertainties. Keywords: regional hydrological model, model uncertainty, Bayesian analysis, Markov Chain Monte Carlo analysis

A novel approach for epipolar resampling of cross-track linear pushbroom imagery using orbital parameters model

Science.gov (United States)

Jannati, Mojtaba; Valadan Zoej, Mohammad Javad; Mokhtarzade, Mehdi

2018-03-01

This paper presents a novel approach to epipolar resampling of cross-track linear pushbroom imagery using orbital parameters model (OPM). The backbone of the proposed method relies on modification of attitude parameters of linear array stereo imagery in such a way to parallelize the approximate conjugate epipolar lines (ACELs) with the instantaneous base line (IBL) of the conjugate image points (CIPs). Afterward, a complementary rotation is applied in order to parallelize all the ACELs throughout the stereo imagery. The new estimated attitude parameters are evaluated based on the direction of the IBL and the ACELs. Due to the spatial and temporal variability of the IBL (respectively changes in column and row numbers of the CIPs) and nonparallel nature of the epipolar lines in the stereo linear images, some polynomials in the both column and row numbers of the CIPs are used to model new attitude parameters. As the instantaneous position of sensors remains fix, the digital elevation model (DEM) of the area of interest is not required in the resampling process. According to the experimental results obtained from two pairs of SPOT and RapidEye stereo imagery with a high elevation relief, the average absolute values of remained vertical parallaxes of CIPs in the normalized images were obtained 0.19 and 0.28 pixels respectively, which confirm the high accuracy and applicability of the proposed method.

Mathematical Modelling for EOQ Inventory System with Advance Payment and Fuzzy Parameters

Directory of Open Access Journals (Sweden)

S Priyan

2014-11-01

Full Text Available This study considers an EOQ inventory model with advance payment policy in a fuzzy situation by employing two types of fuzzy numbers that are trapezoidal and triangular. Two fuzzy models are developed here. In the first model the cost parameters are fuzzified, but the demand rate is treated as crisp constant. In the second model, the demand rate is fuzzified but the cost parameters are treated as crisp constants. For each fuzzy model, we use signed distance method to defuzzify the fuzzy total cost and obtain an estimate of the total cost in the fuzzy sense. Numerical example is provided to ascertain the sensitiveness in the decision variables about fuzziness in the components. In practical situations, costs may be dependent on some foreign monetary unit. In such a case, due to a change in the exchange rates, the costs are often not known precisely. The first model can be used in this situation. In actual applications, demand is uncertain and must be predicted. Accordingly, the decision maker faces a fuzzy environment rather than a stochastic one in these cases. The second model can be used in this situation. Moreover, the proposed models can be expended for imperfect production process.

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

A Bayesian approach for parameter estimation and prediction using a computationally intensive model

International Nuclear Information System (INIS)

Higdon, Dave; McDonnell, Jordan D; Schunck, Nicolas; Sarich, Jason; Wild, Stefan M

2015-01-01

Bayesian methods have been successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model η(θ), where θ denotes the uncertain, best input setting. Hence the statistical model is of the form y=η(θ)+ϵ, where ϵ accounts for measurement, and possibly other, error sources. When nonlinearity is present in η(⋅), the resulting posterior distribution for the unknown parameters in the Bayesian formulation is typically complex and nonstandard, requiring computationally demanding computational approaches such as Markov chain Monte Carlo (MCMC) to produce multivariate draws from the posterior. Although generally applicable, MCMC requires thousands (or even millions) of evaluations of the physics model η(⋅). This requirement is problematic if the model takes hours or days to evaluate. To overcome this computational bottleneck, we present an approach adapted from Bayesian model calibration. This approach combines output from an ensemble of computational model runs with physical measurements, within a statistical formulation, to carry out inference. A key component of this approach is a statistical response surface, or emulator, estimated from the ensemble of model runs. We demonstrate this approach with a case study in estimating parameters for a density functional theory model, using experimental mass/binding energy measurements from a collection of atomic nuclei. We also demonstrate how this approach produces uncertainties in predictions for recent mass measurements obtained at Argonne National Laboratory. (paper)

Mathematical models for correction of images, obtained at radioisotope scan

International Nuclear Information System (INIS)

Glaz, A.; Lubans, A.

2002-01-01

The images, which obtained at radioisotope scintigraphy, contain distortions. Distortions appear as a result of absorption of radiation by patient's body's tissues. Two mathematical models for reducing of such distortions are proposed. Image obtained by only one gamma camera is used in the first mathematical model. Unfortunately, this model allows processing of the images only in case, when it can be assumed, that the investigated organ has a symmetric form. The images obtained by two gamma cameras are used in the second model. It gives possibility to assume that the investigated organ has non-symmetric form and to acquire more precise results. (authors)

Modeling of rheological characteristics of the fermented dairy products obtained by novel and traditional starter cultures.

Science.gov (United States)

Vukić, Dajana V; Vukić, Vladimir R; Milanović, Spasenija D; Ilicić, Mirela D; Kanurić, Katarina G

2018-06-01

Tree different fermented dairy products obtained by conventional and non-conventional starter cultures were investigated in this paper. Textural and rheological characteristics as well as chemical composition during 21 days of storage were analysed and subsequent data processing was performed by principal component analysis. The analysis of samples` flow behaviour was focused on their time dependent properties. Parameters of Power law model described flow behaviour of samples depended on used starter culture and days of storage. The Power law model was applied successfully to describe the flow of the fermented milk, which had characteristics of shear thinning and non-Newtonian fluid behaviour.

Conjugate heat transfer analysis of an energy conversion device with an updated numerical model obtained through inverse identification

International Nuclear Information System (INIS)

Hey, Jonathan; Malloy, Adam C.; Martinez-Botas, Ricardo; Lamperth, Michael

2015-01-01

Highlights: • Conjugate heat transfer analysis of an electric machine. • Inverse identification method for estimating the model parameters. • Experimentally determined thermal properties and electromagnetic losses. • Coupling of inverse identification method with a numerical model. • Improved modeling accuracy through introduction of interface material. - Abstract: Energy conversion devices undergo thermal loading during their operation as a result of inefficiencies in the energy conversion process. This will eventually lead to degradation and possible failure of the device if the heat generated is not properly managed. The ability to accurately predict the thermal behavior of such a device during the initial developmental stage is an important requirement. However, accurate predictions of critical temperature is challenging due to the variation of heat transfer parameters from one device to another. The ability to determine the model parameters is key to accurately representing the heat transfer in such a device. This paper presents the use of an inverse identification technique to estimate the model parameters of an energy conversion device designed for vehicular applications. To simulate the imperfect contact and the presence of insulating materials in the permanent magnet electric machine, thin material are introduced at the component interface of the numerical model. The proposed inverse identification method is used to estimate the equivalent thermal conductance of the thin material. In addition, the electromagnetic losses generated in the permanent magnet is also derived indirectly from the temperature measurement using the same method. With the thermal properties and input parameters of the numerical model obtained from the inverse identification method, the critical temperature of the device can be predicted more accurately. The deviation between the maximum measured and predicted winding temperature is less than 2.4%

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...

Application of the Elitist-Mutated PSO and an Improved GSA to Estimate Parameters of Linear and Nonlinear Muskingum Flood Routing Models.

Directory of Open Access Journals (Sweden)

Ling Kang

Full Text Available Heuristic search algorithms, which are characterized by faster convergence rates and can obtain better solutions than the traditional mathematical methods, are extensively used in engineering optimizations. In this paper, a newly developed elitist-mutated particle swarm optimization (EMPSO technique and an improved gravitational search algorithm (IGSA are successively applied to parameter estimation problems of Muskingum flood routing models. First, the global optimization performance of the EMPSO and IGSA are validated by nine standard benchmark functions. Then, to further analyse the applicability of the EMPSO and IGSA for various forms of Muskingum models, three typical structures are considered: the basic two-parameter linear Muskingum model (LMM, a three-parameter nonlinear Muskingum model (NLMM and a four-parameter nonlinear Muskingum model which incorporates the lateral flow (NLMM-L. The problems are formulated as optimization procedures to minimize the sum of the squared deviations (SSQ or the sum of the absolute deviations (SAD between the observed and the estimated outflows. Comparative results of the selected numerical cases (Case 1-3 show that the EMPSO and IGSA not only rapidly converge but also obtain the same best optimal parameter vector in every run. The EMPSO and IGSA exhibit superior robustness and provide two efficient alternative approaches that can be confidently employed to estimate the parameters of both linear and nonlinear Muskingum models in engineering applications.

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

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

A One-Step-Ahead Smoothing-Based Joint Ensemble Kalman Filter for State-Parameter Estimation of Hydrological Models

KAUST Repository

El Gharamti, Mohamad

2015-11-26

The ensemble Kalman filter (EnKF) recursively integrates field data into simulation models to obtain a better characterization of the model’s state and parameters. These are generally estimated following a state-parameters joint augmentation strategy. In this study, we introduce a new smoothing-based joint EnKF scheme, in which we introduce a one-step-ahead smoothing of the state before updating the parameters. Numerical experiments are performed with a two-dimensional synthetic subsurface contaminant transport model. The improved performance of the proposed joint EnKF scheme compared to the standard joint EnKF compensates for the modest increase in the computational cost.

Estimations of parameters in Pareto reliability model in the presence of masked data

International Nuclear Information System (INIS)

Sarhan, Ammar M.

2003-01-01

Estimations of parameters included in the individual distributions of the life times of system components in a series system are considered in this paper based on masked system life test data. We consider a series system of two independent components each has a Pareto distributed lifetime. The maximum likelihood and Bayes estimators for the parameters and the values of the reliability of the system's components at a specific time are obtained. Symmetrical triangular prior distributions are assumed for the unknown parameters to be estimated in obtaining the Bayes estimators of these parameters. Large simulation studies are done in order: (i) explain how one can utilize the theoretical results obtained; (ii) compare the maximum likelihood and Bayes estimates obtained of the underlying parameters; and (iii) study the influence of the masking level and the sample size on the accuracy of the estimates obtained

Analysis of the spatial variation in the parameters of the SWAT model with application in Flanders, Northern Belgium

Directory of Open Access Journals (Sweden)

G. Heuvelmans

2004-01-01

Full Text Available Operational applications of a hydrological model often require the prediction of stream flow in (future time periods without stream flow observations or in ungauged catchments. Data for a case-specific optimisation of model parameters are not available for such applications, so parameters have to be derived from other catchments or time periods. It has been demonstrated that for applications of the SWAT in Northern Belgium, temporal transfers of the parameters have less influence than spatial transfers on the performance of the model. This study examines the spatial variation in parameter optima in more detail. The aim was to delineate zones wherein model parameters can be transferred without a significant loss of model performance. SWAT was calibrated for 25 catchments that are part of eight larger sub-basins of the Scheldt river basin. Two approaches are discussed for grouping these units in zones with a uniform set of parameters: a single parameter approach considering each parameter separately and a parameter set approach evaluating the parameterisation as a whole. For every catchment, the SWAT model was run with the local parameter optima, with the average parameter values for the entire study region (Flanders, with the zones delineated with the single parameter approach and with the zones obtained by the parameter set approach. Comparison of the model performances of these four parameterisation strategies indicates that both the single parameter and the parameter set zones lead to stream flow predictions that are more accurate than if the entire study region were treated as one single zone. On the other hand, the use of zonal average parameter values results in a considerably worse model fit compared to local parameter optima. Clustering of parameter sets gives a more accurate result than the single parameter approach and is, therefore, the preferred technique for use in the parameterisation of ungauged sub-catchments as part of the

Effect of x-ray tube current on the accuracy of cerebral perfusion parameters obtained by CT perfusion studies

International Nuclear Information System (INIS)

Murase, Kenya; Nanjo, Takafumi; Ii, Satoshi; Miyazaki, Shohei; Hirata, Masaaki; Sugawara, Yoshifumi; Kudo, Masayuki; Sasaki, Kousuke; Mochizuki, Teruhito

2005-01-01

The purpose of this study was to investigate the effect of x-ray tube current on the accuracy of cerebral perfusion parameters obtained by CT perfusion studies using multi-detector row CT (MDCT). Following the standard CT perfusion study protocol, continuous (cine) scans (1 s/rotation x 60 s) consisting of four 5 mm thick contiguous slices were performed using an MDCT scanner with a tube voltage of 80 kVp and a tube current of 200 mA. We generated the simulated images with tube currents of 50 mA, 100 mA and 150 mA by adding the corresponding noise to the raw scan data of the original image acquired above using a noise simulation tool. From the original and simulated images, we generated the functional images of cerebral blood flow (CBF), cerebral blood volume (CBV) and mean transit time (MTT) in seven patients with cerebrovascular disease, and compared the correlation coefficients (CCs) between the perfusion parameter values obtained from the original and simulated images. The coefficients of variation (CVs) in the white matter were also compared. The CC values deteriorated with decreasing tube current. There was a significant difference between 50 mA and 100 mA for all perfusion parameters. The CV values increased with decreasing tube current. There were significant differences between 50 mA and 100 mA and between 100 mA and 150 mA for CBF. For CBV and MTT, there was also a significant difference between 150 mA and 200 mA. This study will be useful for understanding the effect of x-ray tube current on the accuracy of cerebral perfusion parameters obtained by CT perfusion studies using MDCT, and for selecting the tube current

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

Neurofitter: a parameter tuning package for a wide range of electrophysiological neuron models

Directory of Open Access Journals (Sweden)

Werner Van Geit

2007-11-01

Full Text Available The increase in available computational power and the higher quality of experimental recordings have turned the tuning of neuron model parameters into a problem that can be solved by automatic global optimization algorithms. Neurofitter is a software tool that interfaces existing neural simulation software and sophisticated optimization algorithms with a new way to compute the error measure. This error measure represents how well a given parameter set is able to reproduce the experimental data. It is based on the phase-plane trajectory density method, which is insensitive to small phase differences between model and data. Neurofitter enables the effortless combination of many different time-dependent data traces into the error measure, allowing the neuroscientist to focus on what are the seminal properties of the model. We show results obtained by applying Neurofitter to a simple single compartmental model and a complex multi-compartmental Purkinje cell (PC model. These examples show that the method is able to solve a variety of tuning problems and demonstrate details of its practical application.

Reactive decontamination of absorbing thin film polymer coatings: model development and parameter determination

Science.gov (United States)

Varady, Mark; Mantooth, Brent; Pearl, Thomas; Willis, Matthew

2014-03-01

A continuum model of reactive decontamination in absorbing polymeric thin film substrates exposed to the chemical warfare agent O-ethyl S-[2-(diisopropylamino)ethyl] methylphosphonothioate (known as VX) was developed to assess the performance of various decontaminants. Experiments were performed in conjunction with an inverse analysis method to obtain the necessary model parameters. The experiments involved contaminating a substrate with a fixed VX exposure, applying a decontaminant, followed by a time-resolved, liquid phase extraction of the absorbing substrate to measure the residual contaminant by chromatography. Decontamination model parameters were uniquely determined using the Levenberg-Marquardt nonlinear least squares fitting technique to best fit the experimental time evolution of extracted mass. The model was implemented numerically in both a 2D axisymmetric finite element program and a 1D finite difference code, and it was found that the more computationally efficient 1D implementation was sufficiently accurate. The resulting decontamination model provides an accurate quantification of contaminant concentration profile in the material, which is necessary to assess exposure hazards.

The effect of uncertainty of reactor parameters obtained using k0-NAA on result of analysis

International Nuclear Information System (INIS)

Sasajima, Fumio

2006-01-01

Neutron Activation Analysis using the k 0 method is a useful method allowing convenient and accurate simultaneous analysis of plural elements, eliminating the need for the use of comparative reference samples. As already well known, it is essential for the correct result of an analysis to obtain the α-factor and f-factor for a neutron spectrum in an irradiation field accurately when an attempt is made to use the k 0 method. For this reason, based on data obtained from the experiment conducted in the JRR-3 PN-3 system, how uncertainty of the measured values for α-factor and f-factor affects the result of an analysis was evaluated. The process of evaluation involved intentionally varying the values for reactor parameters followed by making an analysis of environmental reference samples (NIST SRM-1632c) using the k 0 method to examine the effect of these factors on the concentrations of 19 elements. The result of the evaluation revealed that the degree of the effect of uncertainty on the concentrations of 19 elements was at best approx. 1% under the condition of this experiment assuming that the factor α, a reactor parameter, had uncertainty of approx. 200%. (author)

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.

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

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 ...

mr. A C++ library for the matching and running of the Standard Model parameters

International Nuclear Information System (INIS)

Kniehl, Bernd A.; Veretin, Oleg L.; Pikelner, Andrey F.; Joint Institute for Nuclear Research, Dubna

2016-01-01

We present the C++ program library mr that allows us to reliably calculate the values of the running parameters in the Standard Model at high energy scales. The initial conditions are obtained by relating the running parameters in the MS renormalization scheme to observables at lower energies with full two-loop precision. The evolution is then performed in accordance with the renormalization group equations with full three-loop precision. Pure QCD corrections to the matching and running are included through four loops. We also provide a Mathematica interface for this program library.

mr. A C++ library for the matching and running of the Standard Model parameters

Energy Technology Data Exchange (ETDEWEB)

Kniehl, Bernd A.; Veretin, Oleg L. [Hamburg Univ. (Germany). II. Inst. fuer Theoretische Physik; Pikelner, Andrey F. [Hamburg Univ. (Germany). II. Inst. fuer Theoretische Physik; Joint Institute for Nuclear Research, Dubna (Russian Federation). Bogoliubov Lab. of Theoretical Physics

2016-01-15

We present the C++ program library mr that allows us to reliably calculate the values of the running parameters in the Standard Model at high energy scales. The initial conditions are obtained by relating the running parameters in the MS renormalization scheme to observables at lower energies with full two-loop precision. The evolution is then performed in accordance with the renormalization group equations with full three-loop precision. Pure QCD corrections to the matching and running are included through four loops. We also provide a Mathematica interface for this program library.

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.

An empirical model for parameters affecting energy consumption in boron removal from boron-containing wastewaters by electrocoagulation.

Science.gov (United States)

Yilmaz, A Erdem; Boncukcuoğlu, Recep; Kocakerim, M Muhtar

2007-06-01

In this study, it was investigated parameters affecting energy consumption in boron removal from boron containing wastewaters prepared synthetically, via electrocoagulation method. The solution pH, initial boron concentration, dose of supporting electrolyte, current density and temperature of solution were selected as experimental parameters affecting energy consumption. The obtained experimental results showed that boron removal efficiency reached up to 99% under optimum conditions, in which solution pH was 8.0, current density 6.0 mA/cm(2), initial boron concentration 100mg/L and solution temperature 293 K. The current density was an important parameter affecting energy consumption too. High current density applied to electrocoagulation cell increased energy consumption. Increasing solution temperature caused to decrease energy consumption that high temperature decreased potential applied under constant current density. That increasing initial boron concentration and dose of supporting electrolyte caused to increase specific conductivity of solution decreased energy consumption. As a result, it was seen that energy consumption for boron removal via electrocoagulation method could be minimized at optimum conditions. An empirical model was predicted by statistically. Experimentally obtained values were fitted with values predicted from empirical model being as following; [formula in text]. Unfortunately, the conditions obtained for optimum boron removal were not the conditions obtained for minimum energy consumption. It was determined that support electrolyte must be used for increase boron removal and decrease electrical energy consumption.

Dynamic analysis of the 7-GeV APS experiment hall foundation based on equivalent lumped parameter modeling

International Nuclear Information System (INIS)

Wambsganss, M.W.

1989-01-01

In this technical note, mass-spring-dashpot, also referred to as equivalent lumped parameter, models are employed to model the soil-foundation interaction of two typical floor segments from the 7-GeV APS experiment hall. Equivalent lumped parameter models have the advantage of being easy to apply and of readily allowing for parameter studies. Analysis requires knowledge of certain properties of the soil including density, shear wave velocity, and Poisson's ratio, as well as knowledge of the degree of homogeneity of the underlying soil stratum. These data for the APS site were determined by a geotechnical investigation. A soil profile and pertinent data, obtained from crosshole seismic testing, are given. Natural frequencies and damping are calculated for the vertical, sliding, rocking, and coupled rocking/sliding modes of vibration. Subsequently, various corrections to account for modeling ''deficiencies'' are considered and their influences evaluated. The equivalent lumped parameter models were developed for machine foundations which, compared with the APS foundation, are smaller in plan dimension. Therefore, the applicability of these models in the analysis of the dynamic characteristics of the APS foundation must be established. The modeling is evaluated by applying the equivalent lumped parameter models in the analysis of large foundations for which test data exists. A comparison of theoretical and test results establishes the basis for an assessment of the applicability and accuracy of the modeling

The x rescaling parameter' formula of the extended x rescaling model and the nuclear effect l-A DIS process

International Nuclear Information System (INIS)

Gao Yonghua; He Mingzhong; Duan Chungui

2003-01-01

The authors present an x rescaling parameters' formula for partons in the extended x rescaling model, where we have established the connection between the rescaling parameter and the mean binding energy. Using x rescaling parameters obtained by the x rescaling parameters' formula, we calculate the average nucleon structure function ratio of DIS process in l-A (Al, Ca, Pb) collision to l-C collision respectively. The result is in good agreement with experimental data

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.

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

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.

Pre-equilibrium assumptions and statistical model parameters effects on reaction cross-section calculations

International Nuclear Information System (INIS)

Avrigeanu, M.; Avrigeanu, V.

1992-02-01

A systematic study on effects of statistical model parameters and semi-classical pre-equilibrium emission models has been carried out for the (n,p) reactions on the 56 Fe and 60 Co target nuclei. The results obtained by using various assumptions within a given pre-equilibrium emission model differ among them more than the ones of different models used under similar conditions. The necessity of using realistic level density formulas is emphasized especially in connection with pre-equilibrium emission models (i.e. with the exciton state density expression), while a basic support could be found only by replacement of the Williams exciton state density formula with a realistic one. (author). 46 refs, 12 figs, 3 tabs

The Effect of Process and Model Parameters in Temperature Prediction for Hot Stamping of Boron Steel

Directory of Open Access Journals (Sweden)

Chaoyang Sun

2013-01-01

Full Text Available Finite element models of the hot stamping and cold die quenching process for boron steel sheet were developed using either rigid or elastic tools. The effect of tool elasticity and process parameters on workpiece temperature was investigated. Heat transfer coefficient between blank and tools was modelled as a function of gap and contact pressure. Temperature distribution and thermal history in the blank were predicted, and thickness distribution of the blank was obtained. Tests were carried out and the test results are used for the validation of numerical predictions. The effect of holding load and the size of cooling ducts on temperature distribution during the forming and the cool die quenching process was also studied by using two models. The results show that higher accuracy predictions of blank thickness and temperature distribution during deformation were obtained using the elastic tool model. However, temperature results obtained using the rigid tool model were close to those using the elastic tool model for a range of holding load.

Sensitivity analysis and parameter estimation for distributed hydrological modeling: potential of variational methods

Directory of Open Access Journals (Sweden)

W. Castaings

2009-04-01

Full Text Available Variational methods are widely used for the analysis and control of computationally intensive spatially distributed systems. In particular, the adjoint state method enables a very efficient calculation of the derivatives of an objective function (response function to be analysed or cost function to be optimised with respect to model inputs.

In this contribution, it is shown that the potential of variational methods for distributed catchment scale hydrology should be considered. A distributed flash flood model, coupling kinematic wave overland flow and Green Ampt infiltration, is applied to a small catchment of the Thoré basin and used as a relatively simple (synthetic observations but didactic application case.

It is shown that forward and adjoint sensitivity analysis provide a local but extensive insight on the relation between the assigned model parameters and the simulated hydrological response. Spatially distributed parameter sensitivities can be obtained for a very modest calculation effort (~6 times the computing time of a single model run and the singular value decomposition (SVD of the Jacobian matrix provides an interesting perspective for the analysis of the rainfall-runoff relation.

For the estimation of model parameters, adjoint-based derivatives were found exceedingly efficient in driving a bound-constrained quasi-Newton algorithm. The reference parameter set is retrieved independently from the optimization initial condition when the very common dimension reduction strategy (i.e. scalar multipliers is adopted.

Furthermore, the sensitivity analysis results suggest that most of the variability in this high-dimensional parameter space can be captured with a few orthogonal directions. A parametrization based on the SVD leading singular vectors was found very promising but should be combined with another regularization strategy in order to prevent overfitting.

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

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)

Parameter estimation of multivariate multiple regression model using bayesian with non-informative Jeffreys’ prior distribution

Science.gov (United States)

Saputro, D. R. S.; Amalia, F.; Widyaningsih, P.; Affan, R. C.

2018-05-01

Bayesian method is a method that can be used to estimate the parameters of multivariate multiple regression model. Bayesian method has two distributions, there are prior and posterior distributions. Posterior distribution is influenced by the selection of prior distribution. Jeffreys’ prior distribution is a kind of Non-informative prior distribution. This prior is used when the information about parameter not available. Non-informative Jeffreys’ prior distribution is combined with the sample information resulting the posterior distribution. Posterior distribution is used to estimate the parameter. The purposes of this research is to estimate the parameters of multivariate regression model using Bayesian method with Non-informative Jeffreys’ prior distribution. Based on the results and discussion, parameter estimation of β and Σ which were obtained from expected value of random variable of marginal posterior distribution function. The marginal posterior distributions for β and Σ are multivariate normal and inverse Wishart. However, in calculation of the expected value involving integral of a function which difficult to determine the value. Therefore, approach is needed by generating of random samples according to the posterior distribution characteristics of each parameter using Markov chain Monte Carlo (MCMC) Gibbs sampling algorithm.

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…

A fundamental parameter-based calibration model for an intrinsic germanium X-ray fluorescence spectrometer

International Nuclear Information System (INIS)

Christensen, L.H.; Pind, N.

1982-01-01

A matrix-independent fundamental parameter-based calibration model for an energy-dispersive X-ray fluorescence spectrometer has been developed. This model, which is part of a fundamental parameter approach quantification method, accounts for both the excitation and detection probability. For each secondary target a number of relative calibration constants are calculated on the basis of knowledge of the irradiation geometry, the detector specifications, and tabulated fundamental physical parameters. The absolute calibration of the spectrometer is performed by measuring one pure element standard per secondary target. For sample systems where all elements can be analyzed by means of the same secondary target the absolute calibration constant can be determined during the iterative solution of the basic equation. Calculated and experimentally determined relative calibration constants agree to within 5-10% of each other and so do the results obtained from the analysis of an NBS certified alloy using the two sets of constants. (orig.)

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...

Modelling of classical ghost images obtained using scattered light

International Nuclear Information System (INIS)

Crosby, S; Castelletto, S; Aruldoss, C; Scholten, R E; Roberts, A

2007-01-01

The images obtained in ghost imaging with pseudo-thermal light sources are highly dependent on the spatial coherence properties of the incident light. Pseudo-thermal light is often created by reducing the coherence length of a coherent source by passing it through a turbid mixture of scattering spheres. We describe a model for simulating ghost images obtained with such partially coherent light, using a wave-transport model to calculate the influence of the scattering on initially coherent light. The model is able to predict important properties of the pseudo-thermal source, such as the coherence length and the amplitude of the residual unscattered component of the light which influence the resolution and visibility of the final ghost image. We show that the residual ballistic component introduces an additional background in the reconstructed image, and the spatial resolution obtainable depends on the size of the scattering spheres

Modelling of classical ghost images obtained using scattered light

Energy Technology Data Exchange (ETDEWEB)

Crosby, S; Castelletto, S; Aruldoss, C; Scholten, R E; Roberts, A [School of Physics, University of Melbourne, Victoria, 3010 (Australia)

2007-08-15

The images obtained in ghost imaging with pseudo-thermal light sources are highly dependent on the spatial coherence properties of the incident light. Pseudo-thermal light is often created by reducing the coherence length of a coherent source by passing it through a turbid mixture of scattering spheres. We describe a model for simulating ghost images obtained with such partially coherent light, using a wave-transport model to calculate the influence of the scattering on initially coherent light. The model is able to predict important properties of the pseudo-thermal source, such as the coherence length and the amplitude of the residual unscattered component of the light which influence the resolution and visibility of the final ghost image. We show that the residual ballistic component introduces an additional background in the reconstructed image, and the spatial resolution obtainable depends on the size of the scattering spheres.

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

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

Refitting density dependent relativistic model parameters including Center-of-Mass corrections

International Nuclear Information System (INIS)

Avancini, Sidney S.; Marinelli, Jose R.; Carlson, Brett Vern

2011-01-01

Full text: Relativistic mean field models have become a standard approach for precise nuclear structure calculations. After the seminal work of Serot and Walecka, which introduced a model Lagrangian density where the nucleons interact through the exchange of scalar and vector mesons, several models were obtained through its generalization, including other meson degrees of freedom, non-linear meson interactions, meson-meson interactions, etc. More recently density dependent coupling constants were incorporated into the Walecka-like models, which are then extensively used. In particular, for these models a connection with the density functional theory can be established. Due to the inherent difficulties presented by field theoretical models, only the mean field approximation is used for the solution of these models. In order to calculate finite nuclei properties in the mean field approximation, a reference set has to be fixed and therefore the translational symmetry is violated. It is well known that in such case spurious effects due to the center-of-mass (COM) motion are present, which are more pronounced for light nuclei. In a previous work we have proposed a technique based on the Pierls-Yoccoz projection operator applied to the mean-field relativistic solution, in order to project out spurious COM contributions. In this work we obtain a new fitting for the density dependent parameters of a density dependent hadronic model, taking into account the COM corrections. Our fitting is obtained taking into account the charge radii and binding energies for He 4 , O 16 , Ca 40 , Ca 48 , Ni 56 , Ni 68 , Sn 100 , Sn 132 and Pb 208 . We show that the nuclear observables calculated using our fit are of a quality comparable to others that can be found in the literature, with the advantage that now a translational invariant many-body wave function is at our disposal. (author)

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

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

A new preprocessing parameter estimation based on geodesic active contour model for automatic vestibular neuritis diagnosis.

Science.gov (United States)

Ben Slama, Amine; Mouelhi, Aymen; Sahli, Hanene; Manoubi, Sondes; Mbarek, Chiraz; Trabelsi, Hedi; Fnaiech, Farhat; Sayadi, Mounir

2017-07-01

The diagnostic of the vestibular neuritis (VN) presents many difficulties to traditional assessment methods This paper deals with a fully automatic VN diagnostic system based on nystagmus parameter estimation using a pupil detection algorithm. A geodesic active contour model is implemented to find an accurate segmentation region of the pupil. Hence, the novelty of the proposed algorithm is to speed up the standard segmentation by using a specific mask located on the region of interest. This allows a drastically computing time reduction and a great performance and accuracy of the obtained results. After using this fast segmentation algorithm, the obtained estimated parameters are represented in temporal and frequency settings. A useful principal component analysis (PCA) selection procedure is then applied to obtain a reduced number of estimated parameters which are used to train a multi neural network (MNN). Experimental results on 90 eye movement videos show the effectiveness and the accuracy of the proposed estimation algorithm versus previous work. Copyright © 2017 Elsevier B.V. All rights reserved.

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

A novel method for identification of lithium-ion battery equivalent circuit model parameters considering electrochemical properties

Science.gov (United States)

Zhang, Xi; Lu, Jinling; Yuan, Shifei; Yang, Jun; Zhou, Xuan

2017-03-01

This paper proposes a novel parameter identification method for the lithium-ion (Li-ion) battery equivalent circuit model (ECM) considering the electrochemical properties. An improved pseudo two-dimension (P2D) model is established on basis of partial differential equations (PDEs), since the electrolyte potential is simplified from the nonlinear to linear expression while terminal voltage can be divided into the electrolyte potential, open circuit voltage (OCV), overpotential of electrodes, internal resistance drop, and so on. The model order reduction process is implemented by the simplification of the PDEs using the Laplace transform, inverse Laplace transform, Pade approximation, etc. A unified second order transfer function between cell voltage and current is obtained for the comparability with that of ECM. The final objective is to obtain the relationship between the ECM resistances/capacitances and electrochemical parameters such that in various conditions, ECM precision could be improved regarding integration of battery interior properties for further applications, e.g., SOC estimation. Finally simulation and experimental results prove the correctness and validity of the proposed methodology.

Effects of 1% cyclopentolate hydrochloride on anterior segment parameters obtained with Pentacam in young adults

Directory of Open Access Journals (Sweden)

Ceyhun Arici

2014-08-01

Full Text Available Purpose: To investigate the effects of topically applied 1% cyclopentolate hydrochloride on anterior segment parameters obtained with a Pentacam rotating Scheimpflug camera in healthy young adults. Methods: Anterior segment analyses of 25 eyes from 25 young adults (Group 1, before and after 45 min of 1% cyclopentolate hydrochloride application, were performed. For a control group (cycloplegia-free, Group 2, 24 eyes of 24 age- and sex-matched healthy cases were evaluated twice at 45 min intervals. The results obtained from the groups were compared statistically. Results: The mean ages of the groups were 23.04 ± 3.42 (range, 18-29 and 22.4 ± 2.05 (range, 18-27 years for Groups 1 and 2, respectively (p=0.259. In Group 1, measurements between the two analyses were significantly different for the values of anterior chamber depth (ACD, anterior chamber angle (ACA, and anterior chamber volume (ACV (p<0.05, whereas no statistical difference was found for the central corneal thickness (CCT and keratometry (K1, K2 measurements. In Group 2, none of these parameters were statistically different between the two analyses. Conclusions: Topically applied 1% cyclopentolate hydrochloride caused an increase in the ACD and ACV values, and a decrease in the ACA value. However, it had no significant effect on the CCT and keratometry measurements. It is important to consider these effects when using the Pentacam device on young adults with cycloplegia and when applying it for various reasons.

Parameter estimation methods for gene circuit modeling from time-series mRNA data: a comparative study.

Science.gov (United States)

Fan, Ming; Kuwahara, Hiroyuki; Wang, Xiaolei; Wang, Suojin; Gao, Xin

2015-11-01

Parameter estimation is a challenging computational problem in the reverse engineering of biological systems. Because advances in biotechnology have facilitated wide availability of time-series gene expression data, systematic parameter estimation of gene circuit models from such time-series mRNA data has become an important method for quantitatively dissecting the regulation of gene expression. By focusing on the modeling of gene circuits, we examine here the performance of three types of state-of-the-art parameter estimation methods: population-based methods, online methods and model-decomposition-based methods. Our results show that certain population-based methods are able to generate high-quality parameter solutions. The performance of these methods, however, is heavily dependent on the size of the parameter search space, and their computational requirements substantially increase as the size of the search space increases. In comparison, online methods and model decomposition-based methods are computationally faster alternatives and are less dependent on the size of the search space. Among other things, our results show that a hybrid approach that augments computationally fast methods with local search as a subsequent refinement procedure can substantially increase the quality of their parameter estimates to the level on par with the best solution obtained from the population-based methods while maintaining high computational speed. These suggest that such hybrid methods can be a promising alternative to the more commonly used population-based methods for parameter estimation of gene circuit models when limited prior knowledge about the underlying regulatory mechanisms makes the size of the parameter search space vastly large. © The Author 2015. Published by Oxford University Press. For Permissions, please email: journals.permissions@oup.com.

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.

Updating parameters of the chicken processing line model

DEFF Research Database (Denmark)

Kurowicka, Dorota; Nauta, Maarten; Jozwiak, Katarzyna

2010-01-01

A mathematical model of chicken processing that quantitatively describes the transmission of Campylobacter on chicken carcasses from slaughter to chicken meat product has been developed in Nauta et al. (2005). This model was quantified with expert judgment. Recent availability of data allows...... updating parameters of the model to better describe processes observed in slaughterhouses. We propose Bayesian updating as a suitable technique to update expert judgment with microbiological data. Berrang and Dickens’s data are used to demonstrate performance of this method in updating parameters...... of the chicken processing line model....

Determination of the atrazine migration parameters in Vertisol

Science.gov (United States)

Raymundo-Raymundo, E.; Hernandez-Vargas, J.; Nikol'Skii, Yu. N.; Guber, A. K.; Gavi-Reyes, F.; Prado-Pano, B. L.; Figueroa-Sandoval, B.; Mendosa-Hernandez, J. R.

2010-05-01

The parameters of the atrazine migration in columns with undisturbed Vertisol sampled from an irrigated plot in Guanajuato, Mexico were determined. A model of the convection-dispersion transport of the chemical compounds accounting for the decomposition and equilibrium adsorption, which is widely applied for assessing the risk of contamination of natural waters with pesticides, was used. The model parameters were obtained by solving the inverse problem of the transport equation on the basis of laboratory experiments on the transport of the 18O isotope and atrazine in soil columns with an undisturbed structure at three filtration velocities. The model adequately described the experimental data at the individual selection of the parameters for each output curve. Physically unsubstantiated parameters of the atrazine adsorption and degradation were obtained when the parameter of the hydrodynamic dispersion was determined from the data on the 18O migration. The simulation also showed that the use of parameters obtained at water content close to saturation in the calculations for an unsaturated soil resulted in the overestimation of the leaching rate and the maximum concentration of atrazine in the output curve compared to the experimental data.

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

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.

Estimating model parameters for an impact-produced shock-wave simulation: Optimal use of partial data with the extended Kalman filter

International Nuclear Information System (INIS)

Kao, Jim; Flicker, Dawn; Ide, Kayo; Ghil, Michael

2006-01-01

This paper builds upon our recent data assimilation work with the extended Kalman filter (EKF) method [J. Kao, D. Flicker, R. Henninger, S. Frey, M. Ghil, K. Ide, Data assimilation with an extended Kalman filter for an impact-produced shock-wave study, J. Comp. Phys. 196 (2004) 705-723.]. The purpose is to test the capability of EKF in optimizing a model's physical parameters. The problem is to simulate the evolution of a shock produced through a high-speed flyer plate. In the earlier work, we have showed that the EKF allows one to estimate the evolving state of the shock wave from a single pressure measurement, assuming that all model parameters are known. In the present paper, we show that imperfectly known model parameters can also be estimated accordingly, along with the evolving model state, from the same single measurement. The model parameter optimization using the EKF can be achieved through a simple modification of the original EKF formalism by including the model parameters into an augmented state variable vector. While the regular state variables are governed by both deterministic and stochastic forcing mechanisms, the parameters are only subject to the latter. The optimally estimated model parameters are thus obtained through a unified assimilation operation. We show that improving the accuracy of the model parameters also improves the state estimate. The time variation of the optimized model parameters results from blending the data and the corresponding values generated from the model and lies within a small range, of less than 2%, from the parameter values of the original model. The solution computed with the optimized parameters performs considerably better and has a smaller total variance than its counterpart using the original time-constant parameters. These results indicate that the model parameters play a dominant role in the performance of the shock-wave hydrodynamic code at hand

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

Using Mathematical Modeling Methods for Estimating Entrance Flow Heterogeneity Impact on Aviation GTE Parameters and Performances

Directory of Open Access Journals (Sweden)

Yu. A. Ezrokhi

2017-01-01

Full Text Available The paper considers methodological approaches to the mathematical models (MM of various levels, dedicated to estimate an impact of the entrance flow heterogeneity on the main parameters and performances of the aviation GTE and it units. By an example of calculation of a twin-shaft turbofan engine in cruiser mode, demonstrates engineering mathematical model capabilities to define the impact of the total pressure field distortion on engine trust and air flow parameters, and also gas dynamic stability margin of the both compressors.It is shown that the presented first level mathematical model allows us to estimate sufficiently the impact of entrance total pressure heterogeneity on the engine parameters. Here reliability of calculations is proved to be true by their comparison with the results, obtained owing to well fulfilled 2D & 3D mathematical models of the engine, which have been repeatedly identified by the results of experiments.It is shown that received results including those on decreasing values of stability margin of both compressors can be used for tentative estimates when choosing a desirable stability margin, providing steady operation of compressors and engine in an entire range of its operating modes. Carrying out a definitive testing calculation using the specialized engine MM of a higher level will not only confirm the results obtained, but also reduce their expected error with regard to the real values reached as a result of tests.

Setting Parameters for Biological Models With ANIMO

NARCIS (Netherlands)

Schivo, Stefano; Scholma, Jetse; Karperien, Hermanus Bernardus Johannes; Post, Janine Nicole; van de Pol, Jan Cornelis; Langerak, Romanus; André, Étienne; Frehse, Goran

2014-01-01

ANIMO (Analysis of Networks with Interactive MOdeling) is a software for modeling biological networks, such as e.g. signaling, metabolic or gene networks. An ANIMO model is essentially the sum of a network topology and a number of interaction parameters. The topology describes the interactions

Parameter Estimation of Dynamic Multi-zone Models for Livestock Indoor Climate Control

DEFF Research Database (Denmark)

Wu, Zhuang; Stoustrup, Jakob; Heiselberg, Per

2008-01-01

, the livestock, the ventilation system and the building on the dynamic performance of indoor climate. Some significant parameters employed in the climate model as well as the airflow interaction between each conceptual zone are identified with the use of experimental time series data collected during spring......In this paper, a multi-zone modeling concept is proposed based on a simplified energy balance formulation to provide a better prediction of the indoor horizontal temperature variation inside the livestock building. The developed mathematical models reflect the influences from the weather...... and winter at a real scale livestock building in Denmark. The obtained comparative results between the measured data and the simulated output confirm that a very simple multi-zone model can capture the salient dynamical features of the climate dynamics which are needed for control purposes....

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.

Parameters estimation of the single and double diode photovoltaic models using a Gauss–Seidel algorithm and analytical method: A comparative study

International Nuclear Information System (INIS)

Et-torabi, K.; Nassar-eddine, I.; Obbadi, A.; Errami, Y.; Rmaily, R.; Sahnoun, S.; El fajri, A.; Agunaou, M.

2017-01-01

Highlights: • Comparative study of two methods: a Gauss Seidel method and an analytical method. • Five models are implemented to estimate the five parameters for single diode. • Two models are used to estimate the seven parameters for double diode. • The parameters are estimated under changing environmental conditions. • To choose method/model combination more adequate for each PV module technology. - Abstract: In the photovoltaic (PV) panels modeling field, this paper presents a comparative study of two parameter estimation methods: the iterative method called Gauss Seidel, applied on the single diode model, and the analytical method used on the double diode model. These parameter estimation methods are based on the manufacturer's datasheets. They are also tested on three PV modules of different technologies: multicrystalline (kyocera KC200GT), monocrystalline (Shell SQ80), and thin film (Shell ST40). For the iterative method, five existing mathematical models classified from 1 to 5 are used to estimate the parameters of these PV modules under varying environmental conditions. Only two models of them are used for the analytical method. Each model is based on the combination of the photocurrent and the reverse saturation current’s expressions in terms of temperature and irradiance. In addition, the results of the models’ simulation are compared with the experimental data obtained from the PV modules’ datasheets, in order to evaluate the accuracy of the models. The simulation shows that the I-V characteristics obtained are matching to the experimental data. In order to validate the reliability of the two methods, both the Absolute Error (AE) and the Root Mean Square Error (RMSE) were calculated. The results suggest that the analytical method can be very useful for monocrystalline and multicrystalline modules, but for the thin film module, the iterative method is the most suitable.

Impact of ambiguity resolution and application of transformation parameters obtained by regional GNSS network in Precise Point Positioning

Science.gov (United States)

Gandolfi, S.; Poluzzi, L.; Tavasci, L.

2012-12-01

Precise Point Positioning (PPP) is one of the possible approaches for GNSS data processing. As known this technique is faster and more flexible compared to the others which are based on a differenced approach and constitute a reliable methods for accurate positioning of remote GNSS stations, even in some remote area such as Antarctica. Until few years ago one of the major limits of the method was the impossibility to resolve the ambiguity as integer but nowadays many methods are available to resolve this aspect. The first software package permitting a PPP solution was the GIPSY OASIS realized, developed and maintained by JPL (NASA). JPL produce also orbits and files ready to be used with GIPSY. Recently, using these products came possible to resolve ambiguities improving the stability of solutions. PPP permit to estimate position into the reference frame of the orbits (IGS) and when coordinate in others reference frames, such al ITRF, are needed is necessary to apply a transformation. Within his products JPL offer, for each day, a global 7 parameter transformation that permit to locate the survey into the ITRF RF. In some cases it's also possible to create a costumed process and obtain analogous parameters using local/regional reference network of stations which coordinates are available also in the desired reference frame. In this work some tests on accuracy has been carried out comparing different PPP solutions obtained using the same software packages (GIPSY) but considering the ambiguity resolution, the global and regional transformation parameters. In particular two test area have been considered, first one located in Antarctica and the second one in Italy. Aim of the work is the evaluation of the impact of ambiguity resolution and the use of local/regional transformation parameter in the final solutions. Tests shown how the ambiguity resolution improve the precision, especially in the EAST component with a scattering reduction about 8%. And the use of global

Automatic J–A Model Parameter Tuning Algorithm for High Accuracy Inrush Current Simulation

Directory of Open Access Journals (Sweden)

Xishan Wen

2017-04-01

Full Text Available Inrush current simulation plays an important role in many tasks of the power system, such as power transformer protection. However, the accuracy of the inrush current simulation can hardly be ensured. In this paper, a Jiles–Atherton (J–A theory based model is proposed to simulate the inrush current of power transformers. The characteristics of the inrush current curve are analyzed and results show that the entire inrush current curve can be well featured by the crest value of the first two cycles. With comprehensive consideration of both of the features of the inrush current curve and the J–A parameters, an automatic J–A parameter estimation algorithm is proposed. The proposed algorithm can obtain more reasonable J–A parameters, which improve the accuracy of simulation. Experimental results have verified the efficiency of the proposed algorithm.

Investigation on Insar Time Series Deformation Model Considering Rheological Parameters for Soft Clay Subgrade Monitoring

Science.gov (United States)

Xing, X.; Yuan, Z.; Chen, L. F.; Yu, X. Y.; Xiao, L.

2018-04-01

The stability control is one of the major technical difficulties in the field of highway subgrade construction engineering. Building deformation model is a crucial step for InSAR time series deformation monitoring. Most of the InSAR deformation models for deformation monitoring are pure empirical mathematical models, without considering the physical mechanism of the monitored object. In this study, we take rheology into consideration, inducing rheological parameters into traditional InSAR deformation models. To assess the feasibility and accuracy for our new model, both simulation and real deformation data over Lungui highway (a typical highway built on soft clay subgrade in Guangdong province, China) are investigated with TerraSAR-X satellite imagery. In order to solve the unknows of the non-linear rheological model, three algorithms: Gauss-Newton (GN), Levenberg-Marquarat (LM), and Genetic Algorithm (GA), are utilized and compared to estimate the unknown parameters. Considering both the calculation efficiency and accuracy, GA is chosen as the final choice for the new model in our case study. Preliminary real data experiment is conducted with use of 17 TerraSAR-X Stripmap images (with a 3-m resolution). With the new deformation model and GA aforementioned, the unknown rheological parameters over all the high coherence points are obtained and the LOS deformation (the low-pass component) sequences are generated.

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

Intelligent Modeling Combining Adaptive Neuro Fuzzy Inference System and Genetic Algorithm for Optimizing Welding Process Parameters

Science.gov (United States)

Gowtham, K. N.; Vasudevan, M.; Maduraimuthu, V.; Jayakumar, T.

2011-04-01

Modified 9Cr-1Mo ferritic steel is used as a structural material for steam generator components of power plants. Generally, tungsten inert gas (TIG) welding is preferred for welding of these steels in which the depth of penetration achievable during autogenous welding is limited. Therefore, activated flux TIG (A-TIG) welding, a novel welding technique, has been developed in-house to increase the depth of penetration. In modified 9Cr-1Mo steel joints produced by the A-TIG welding process, weld bead width, depth of penetration, and heat-affected zone (HAZ) width play an important role in determining the mechanical properties as well as the performance of the weld joints during service. To obtain the desired weld bead geometry and HAZ width, it becomes important to set the welding process parameters. In this work, adaptative neuro fuzzy inference system is used to develop independent models correlating the welding process parameters like current, voltage, and torch speed with weld bead shape parameters like depth of penetration, bead width, and HAZ width. Then a genetic algorithm is employed to determine the optimum A-TIG welding process parameters to obtain the desired weld bead shape parameters and HAZ width.

Ammonium removal from aqueous solutions by clinoptilolite: determination of isotherm and thermodynamic parameters and comparison of kinetics by the double exponential model and conventional kinetic models.

Science.gov (United States)

Tosun, Ismail

2012-03-01

The adsorption isotherm, the adsorption kinetics, and the thermodynamic parameters of ammonium removal from aqueous solution by using clinoptilolite in aqueous solution was investigated in this study. Experimental data obtained from batch equilibrium tests have been analyzed by four two-parameter (Freundlich, Langmuir, Tempkin and Dubinin-Radushkevich (D-R)) and four three-parameter (Redlich-Peterson (R-P), Sips, Toth and Khan) isotherm models. D-R and R-P isotherms were the models that best fitted to experimental data over the other two- and three-parameter models applied. The adsorption energy (E) from the D-R isotherm was found to be approximately 7 kJ/mol for the ammonium-clinoptilolite system, thereby indicating that ammonium is adsorbed on clinoptilolite by physisorption. Kinetic parameters were determined by analyzing the nth-order kinetic model, the modified second-order model and the double exponential model, and each model resulted in a coefficient of determination (R(2)) of above 0.989 with an average relative error lower than 5%. A Double Exponential Model (DEM) showed that the adsorption process develops in two stages as rapid and slow phase. Changes in standard free energy (∆G°), enthalpy (∆H°) and entropy (∆S°) of ammonium-clinoptilolite system were estimated by using the thermodynamic equilibrium coefficients.

Ammonium Removal from Aqueous Solutions by Clinoptilolite: Determination of Isotherm and Thermodynamic Parameters and Comparison of Kinetics by the Double Exponential Model and Conventional Kinetic Models

Directory of Open Access Journals (Sweden)

İsmail Tosun

2012-03-01

Full Text Available The adsorption isotherm, the adsorption kinetics, and the thermodynamic parameters of ammonium removal from aqueous solution by using clinoptilolite in aqueous solution was investigated in this study. Experimental data obtained from batch equilibrium tests have been analyzed by four two-parameter (Freundlich, Langmuir, Tempkin and Dubinin-Radushkevich (D-R and four three-parameter (Redlich-Peterson (R-P, Sips, Toth and Khan isotherm models. D-R and R-P isotherms were the models that best fitted to experimental data over the other two- and three-parameter models applied. The adsorption energy (E from the D-R isotherm was found to be approximately 7 kJ/mol for the ammonium-clinoptilolite system, thereby indicating that ammonium is adsorbed on clinoptilolite by physisorption. Kinetic parameters were determined by analyzing the nth-order kinetic model, the modified second-order model and the double exponential model, and each model resulted in a coefficient of determination (R2 of above 0.989 with an average relative error lower than 5%. A Double Exponential Model (DEM showed that the adsorption process develops in two stages as rapid and slow phase. Changes in standard free energy (∆G°, enthalpy (∆H° and entropy (∆S° of ammonium-clinoptilolite system were estimated by using the thermodynamic equilibrium coefficients.

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.

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

Adaptive Parameter Optimization of a Grid-based Conceptual Hydrological Model

Science.gov (United States)

Samaniego, L.; Kumar, R.; Attinger, S.

2007-12-01

). The main difference with the standard SA is the parameter search routine which uses adaptive heuristic rules to improve its efficiency. These rules are based on the relative behavior of the efficiency criteria. The efficiency of the model is evaluated with the Nash-Sutcliffe efficiency coefficient (NS) and the RMSE obtained for various short and long term runoff characteristics such as daily flows; semiannual high and low flow characteristics such as total drought duration frequency of high flows; and annual specific discharge at various gauging stations. Additionally, the parameter search was constrained with the 95% confidence bands of the runoff characteristics mentioned above. The proposed method was calibrated in the Upper Neckar River basin covering an area of approximately 4000~km2 during the period from 1961 to 1993. The spatial and temporal resolutions used were a grid size of (1000 × 1000)~m and 12~h intervals respectively. The results of the study indicate significant improvement in model performance (e.g. Nash-Sutcliffe of various runoff characteristics ~ 0.8) and a significant reduction in computational burden of at least 25%.

Identification of biomolecule mass transport and binding rate parameters in living cells by inverse modeling

Directory of Open Access Journals (Sweden)

Shirmohammadi Adel

2006-10-01

Full Text Available Abstract Background Quantification of in-vivo biomolecule mass transport and reaction rate parameters from experimental data obtained by Fluorescence Recovery after Photobleaching (FRAP is becoming more important. Methods and results The Osborne-Moré extended version of the Levenberg-Marquardt optimization algorithm was coupled with the experimental data obtained by the Fluorescence Recovery after Photobleaching (FRAP protocol, and the numerical solution of a set of two partial differential equations governing macromolecule mass transport and reaction in living cells, to inversely estimate optimized values of the molecular diffusion coefficient and binding rate parameters of GFP-tagged glucocorticoid receptor. The results indicate that the FRAP protocol provides enough information to estimate one parameter uniquely using a nonlinear optimization technique. Coupling FRAP experimental data with the inverse modeling strategy, one can also uniquely estimate the individual values of the binding rate coefficients if the molecular diffusion coefficient is known. One can also simultaneously estimate the dissociation rate parameter and molecular diffusion coefficient given the pseudo-association rate parameter is known. However, the protocol provides insufficient information for unique simultaneous estimation of three parameters (diffusion coefficient and binding rate parameters owing to the high intercorrelation between the molecular diffusion coefficient and pseudo-association rate parameter. Attempts to estimate macromolecule mass transport and binding rate parameters simultaneously from FRAP data result in misleading conclusions regarding concentrations of free macromolecule and bound complex inside the cell, average binding time per vacant site, average time for diffusion of macromolecules from one site to the next, and slow or rapid mobility of biomolecules in cells. Conclusion To obtain unique values for molecular diffusion coefficient and

Parameter estimation and statistical test of geographically weighted bivariate Poisson inverse Gaussian regression models

Science.gov (United States)

Amalia, Junita; Purhadi, Otok, Bambang Widjanarko

2017-11-01

Poisson distribution is a discrete distribution with count data as the random variables and it has one parameter defines both mean and variance. Poisson regression assumes mean and variance should be same (equidispersion). Nonetheless, some case of the count data unsatisfied this assumption because variance exceeds mean (over-dispersion). The ignorance of over-dispersion causes underestimates in standard error. Furthermore, it causes incorrect decision in the statistical test. Previously, paired count data has a correlation and it has bivariate Poisson distribution. If there is over-dispersion, modeling paired count data is not sufficient with simple bivariate Poisson regression. Bivariate Poisson Inverse Gaussian Regression (BPIGR) model is mix Poisson regression for modeling paired count data within over-dispersion. BPIGR model produces a global model for all locations. In another hand, each location has different geographic conditions, social, cultural and economic so that Geographically Weighted Regression (GWR) is needed. The weighting function of each location in GWR generates a different local model. Geographically Weighted Bivariate Poisson Inverse Gaussian Regression (GWBPIGR) model is used to solve over-dispersion and to generate local models. Parameter estimation of GWBPIGR model obtained by Maximum Likelihood Estimation (MLE) method. Meanwhile, hypothesis testing of GWBPIGR model acquired by Maximum Likelihood Ratio Test (MLRT) method.

Model of twelve properties of a set of organic solvents with graph-theoretical and/or experimental parameters.

Science.gov (United States)

Pogliani, Lionello

2010-01-30

Twelve properties of a highly heterogeneous class of organic solvents have been modeled with a graph-theoretical molecular connectivity modified (MC) method, which allows to encode the core electrons and the hydrogen atoms. The graph-theoretical method uses the concepts of simple, general, and complete graphs, where these last types of graphs are used to encode the core electrons. The hydrogen atoms have been encoded by the aid of a graph-theoretical perturbation parameter, which contributes to the definition of the valence delta, delta(v), a key parameter in molecular connectivity studies. The model of the twelve properties done with a stepwise search algorithm is always satisfactory, and it allows to check the influence of the hydrogen content of the solvent molecules on the choice of the type of descriptor. A similar argument holds for the influence of the halogen atoms on the type of core electron representation. In some cases the molar mass, and in a minor way, special "ad hoc" parameters have been used to improve the model. A very good model of the surface tension could be obtained by the aid of five experimental parameters. A mixed model method based on experimental parameters plus molecular connectivity indices achieved, instead, to consistently improve the model quality of five properties. To underline is the importance of the boiling point temperatures as descriptors in these last two model methodologies. Copyright 2009 Wiley Periodicals, Inc.

Admissible Estimators in the General Multivariate Linear Model with Respect to Inequality Restricted Parameter Set

Directory of Open Access Journals (Sweden)

Shangli Zhang

2009-01-01

Full Text Available By using the methods of linear algebra and matrix inequality theory, we obtain the characterization of admissible estimators in the general multivariate linear model with respect to inequality restricted parameter set. In the classes of homogeneous and general linear estimators, the necessary and suffcient conditions that the estimators of regression coeffcient function are admissible are established.

Connecting Majorana phases to the geometric parameters of the Majorana unitarity triangle in a neutrino mass matrix model

Science.gov (United States)

Verma, Surender; Bhardwaj, Shankita

2018-05-01

We have investigated a possible connection between the Majorana phases and geometric parameters of Majorana unitarity triangle (MT) in two-texture zero neutrino mass matrix. Such analytical relations can, also, be obtained for other theoretical models viz. hybrid textures, neutrino mass matrix with vanishing minors and have profound implications for geometric description of C P violation. As an example, we have considered the two-texture zero neutrino mass model to obtain a relation between Majorana phases and MT parameters that may be probed in various lepton number violating processes. In particular, we find that Majorana phases depend on only one of the three interior angles of the MT in each class of two-texture zero neutrino mass matrix. We have also constructed the MT for class A , B , and C neutrino mass matrices. Nonvanishing areas and nontrivial orientations of these Majorana unitarity triangles indicate nonzero C P violation as a generic feature of this class of mass models.

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.)

Role of deceleration parameter and interacting dark energy in singularity avoidance

Science.gov (United States)

Abdussattar; Prajapati, S. R.

2011-02-01

A class of non-singular bouncing FRW models are obtained by constraining the deceleration parameter in the presence of an interacting dark energy represented by a time-varying cosmological constant. The models being geometrically closed, initially accelerate for a certain period of time and decelerate thereafter and are also free from the entropy and cosmological constant problems. Taking a constant of integration equal to zero one particular model is discussed in some detail and the variation of different cosmological parameters are shown graphically for specific values of the parameters of the model. For some specific choice of the parameters of the model the ever expanding models of Ozer & Taha and Abdel-Rahman and the decelerating models of Berman and also the Einstein de-Sitter model may be obtained as special cases of this particular model.

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.

Solliton-like order parameter distributions in the critical region

Directory of Open Access Journals (Sweden)

A.V.Babich

2006-01-01

Full Text Available Some exact one-component order parameter distributions for the Michelson thermodynamic potential are obtained. The phase transition of second kind in Ginzburg-Landau type model is investigated. The exact partial distribution of the order parameter in the form of Jakobi elliptic function is obtained. The energy of this distribution is lower at some temperature interval than for the best known models.

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 practical method to assess model sensitivity and parameter uncertainty in C cycle models

Science.gov (United States)

Delahaies, Sylvain; Roulstone, Ian; Nichols, Nancy

2015-04-01

data streams or by considering longer observation windows no systematic analysis has been carried out so far to explain the large differences among results. We consider adjoint based methods to investigate inverse problems using DALEC and various data streams. Using resolution matrices we study the nature of the inverse problems (solution existence, uniqueness and stability) and show how standard regularization techniques affect resolution and stability properties. Instead of using standard prior information as a penalty term in the cost function to regularize the problems we constraint the parameter space using ecological balance conditions and inequality constraints. The efficiency and rapidity of this approach allows us to compute ensembles of solutions to the inverse problems from which we can establish the robustness of the variational method and obtain non Gaussian posterior distributions for the model parameters and initial carbon stocks.

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