Strike Point Control on EAST Using an Isoflux Control Method
|
|
- Gerald McKenzie
- 6 years ago
- Views:
Transcription
1 Plasma Science and Technology, Vol.17, No.9, Sep Strike Point Control on EAST Using an Isoflux Control Method XING Zhe ( ) 1, XIAO Bingjia ( ) 1,2, LUO Zhengping ( ) 1, M. L. WALKER 3, D. A. HUMPHREYS 3 1 Institute of Plasma Physics, Chinese Academy of Sciences, Hefei , China 2 School of Nuclear Science and Technology, University of Science and Technology of China, Hefei , China 3 General Atomics, General Atomics Ct., San Diego, CA 92186, USA Abstract For the advanced tokamak, the particle deposition and thermal load on the divertor is a big challenge. By moving the strike points on divertor target plates, the position of particle deposition and thermal load can be shifted. We could adjust the Poloidal Field (PF) coil current to achieve the strike point position feedback control. Using isoflux control method, the strike point position can be controlled by controlling the X point position. On the basis of experimental data, we establish relational expressions between X point position and strike point position. Benchmark experiments are carried out to validate the correctness and robustness of the control methods. The strike point position is successfully controlled following our command in the EAST operation. Keywords: EAST, plasma control system, strike point control PACS: Dy, Fa DOI: / /17/9/09 (Some figures may appear in colour only in the online journal) 1 Introduction In the advanced tokamak, the plasma is confined by a magnetic field [1] with elongated cross-section [2]. The plasma is directed to the divertor region, which is designed to absorb the heat of the energetic byproduct particles, and then impacts on special plates at an appropriate distance far from the hot plasma core [3]. These plates are referred to as the divertor target plates. Fig. 1 illustrates the EAST geometry and one typical plasma shape. To better adjust the particle deposition and thermal load, the magnetic field should be varied in order to change the strike point position on the divertor target plates. For advanced tokamak operation, a reliable control of the strike point position is necessary for the study of plasma facing materials on the divertor target plates. Furthermore, the strike point control should be integrated into the plasma control system as a part of the plasma shape feedback controller [4]. Moving the strike point position to reduce the local average thermal loading has been demonstrated on JET [5], NSTX [6], and DIII-D [7]. Considering that the Plasma Control System (PCS) on EAST is adopted from DIII-D [8,9], we apply a similar realtime strike point control algorithm using rtefit-isoflux method [10]. At first, we establish relational expressions between X point position and strike point position for the strike point control based on historical experimental data. Benchmark experiments were then carried out to validate the correctness and robustness of the control methods. Finally, the strike point position was successfully controlled following our command in the EAST operation. Our results are discussed in the summary. Fig.1 The plasma shape (black line), the vacuum vessel (green line), the four divertor target plates (blue line) on the first wall (magenta line), the control segments (gray line), and the region for searching X points (gray grids) in rtefit supported by the National Magnetic Confinement Fusion Science Program of China (Nos. 2012GB and 2014GB103000) 774
2 XING Zhe et al.: Strike Point Control on EAST Using an Isoflux Control Method 2 rtefit-isoflux algorithm on EAST To achieve shape control in a tokamak, we need to gather information about the plasma shape during discharge. The plasma shape parameters cannot be measured directly. We may reconstruct the plasma boundary, X point locations and some other parameters from the available diagnostic data, such as external magnetic measurements. By now, many efficient methods and numerical codes for equilibrium reconstruction have been developed. In EAST, the real time equilibrium reconstruction code (rtefit) [10] is applied to obtain plasma shape parameters. The rtefit is the real time version of the equilibrium reconstruction code EFIT [11]. EFIT code provides the approximate equilibrium solution to the Grad-Shafranov equation which best fits the external diagnostic data. The rtefit algorithm is integrated into the isoflux shape control algorithms for PCS [8,12]. In the isoflux control scheme [13,14], we set a series of the control segments (the gray bars in Fig. 1) which intersect the plasma last closed flux surface at the control points. These control points define the top and bottom of the plasma shape, and also the outer and inner gap. Particularly, in the case of divertor shape discharges, the control points include the X points. In the shape controller, one of the control points is chosen as the reference point. For the divertor discharge, the X point is typically used as the reference point. If the X point is not unique, such as double null shape, then the reference point is the one with the higher poloidal flux. The poloidal flux differences between the reference point and each of the other control points are used as error signals to control the plasma boundary. The X point position differences between the computed and the desired location are used as error signals to control the X points. These error signals will be translated to the control command for current adjustments in the poloidal field coils. In the isoflux control method, a preliminary method for achieving strike point control has been included [14]. Taking the inner strike point control method as an example, the desired position on the inner divertor plate and the radial position of X point are selected as control points. The strike point position error is generated as an control error signal, which is used to control the vertical position of X point instead of the X point vertical position error. Consequently, we can control the strike point by controlling the X point position. To achieve the strike point control, we need to fit the relational expressions between the X point position and the strike point position by using statistical tools based on historical experimental data. 3 Strike point position control using isoflux control method There are three kinds of strike point control mode during discharge [14] : 775 (I) Inner strike point control; (II) Outer strike point control; and, (III) Both inner and outer strike point control. Considering the three strike point control modes, we need to fit three groups of expressions to determine the X point locations with the desired strike point positions: (I) Programming the inner strike point positions and the radial position of the X point to compute vertical position of the X point: zxtop = zxtop (rxtop, z strike in, r strike in); zxbot = zxbot (rxbot, z strike in, r strike in). (II) Programming the outer strike point positions and the vertical position of the X point to compute radial position of the X point: rxtop = rxtop (zxtop, z strike out, rxbot = rxbot (zxbot, z strike out, r strike out). (III) Programming both inner strike point and outer strike point positions to compute the X point positions: rxtop = rxtop (z strike in, r strike in, z strike out, zxtop = zxtop (z strike in, r strike in, z strike out, rxbot = rxbot (z strike in, r strike in, z strike out, zxbot = zxbot (z strike in, r strike in, z strike out, where the top X point depends on the upper strike points, and the bottom X point depends on the lower strike points. As noted above, the plasma shape parameters, except for the strike point position, can be evaluated by rtefit. To fit the relational expressions between the X point position and the strike point position, we need to determine the strike point position based on the equilibrium data provided by rtefit. We can then collect a large enough sample from the historical experimental data. 3.1 Determination of strike point position based on rtefit data We retrieve necessary data, such as the poloidal flux at the plasma boundary, the poloidal flux on the rectangular grid points (33 33 rectangular grids shown in Fig. 2), and the corresponding X point positions. Based on the flux values of the rectangular grid points, we adopt an interpolation technique to compute the flux values on the divertor target plates (the area of target plate is divided to rectangular grids). The grid point on the divertor target plate with the flux value approximately equal to the plasma boundary flux is considered as the strike point. For this work, we choose cubic spline interpolation for the two-dimensional interpolation. Spline interpolation uses low-degree polynomials in each of the intervals, and chooses the polynomial pieces such that they fit smoothly together. Like polynomial interpolation, spline interpolation incurs a smaller error than linear interpolation, and the interpolation function is
3 Plasma Science and Technology, Vol.17, No.9, Sep smoother. To test the algorithm, we take the results of shot at s as an example. The estimated result of the strike point positions is plotted (the red hexagrams in Fig. 3) with the poloidal flux contour of the plasma boundary as a reference. The details of the four strike points are shown in Fig. 4. Fig.3 The estimated results of the strike points (four red hexagrams on the divertor target plates). The solid line is the poloidal flux contour of the plasma boundary Fig.2 The plasma shape (black line), the vacuum vessel (green line), the four divertor target plates (black line) on the first wall (blue line). The grids (gray line) are the equilibrium reconstruction area of rtefit 3.2 Fitting function We collect 450 sets of (strike point positions, X point positions) as the sample data to fit the relational expressions between X point position and strike point position. The sample distribution is shown in Fig. 5. (a) Upper strike point on inner plate (b) Upper strike point on inner plate (c) Lower strike point on inner plate (d) Lower strike point on outer plate Fig.4 The details of the four strike points (red hexagram). The distance between the grid point (red circle) of the inner divertor target plate is 0.56 mm, while the outer one is 0.49 mm 776
4 XING Zhe et al.: Strike Point Control on EAST Using an Isoflux Control Method +d [(R strike out R strike in ) 2 +(Z strike out Z strike in ) 2 ], (5) Z x = a + b (Z strike out + Z strike in ) + c (R strike out R strike in ) 2 + (Z strike out Z strike in ) 2 +d [(R strike out R strike in ) 2 +(Z strike out Z strike in ) 2 ]. (6) Taking the expression form Eq. (1) as an example, if we choose z strike in, (rxbot-r strike in) and SQUARE (rxbot-r strike in) as three independent variables, we can do multiple linear fitting instead of nonlinear fitting. Based on the 450 sets of sample data, the coefficients of relational expressions for the three strike point control modes are obtained (shown in Table 1). (a) Sample data in the upper divertor region Table 1. The parameters of the fitting functions for the three strike point control modes Control a b c d mode I Zx top Zx bottom II Rx top Rx bottom III Rx top Zx top Rx bottom Zx bottom (b) Sample data in the lower divertor region Fig.5 The distribution of sample data for fitting For the strike point control mode (I) as noted above: Z xtop = a + b Z strike + c (R xtop R strike ) +d (R xtop R strike ) 2, (1) Z xbot = a + b Z strike + c (R xbot R strike ) +d (R xbot R strike ) 2. (2) For the strike point control mode (II) as noted above: R xtop = a + b R strike + c (Z strike Z xtop ) +d (Z strike Z xtop ) 2, (3) R xbot = a + b R strike + c (Z strike Z xbot ) +d (Z strike Z xbot ) 2. (4) For the strike point control mode (III) as noted above: R x = a + b (R strike out + R strike in ) + c (R strike out R strike in ) 2 + (Z strike out Z strike in ) Benchmark test To benchmark the relational expressions, we selected ten time slices from the experimental data at random, and calculated the relative errors of X point positions between the rtefit results and the estimated results by using the relational expressions. As shown in Fig. 6, the relative errors are less than 1%. We connect the strike point and the estimated X point with a straight line, which is easy for us to examine whether or not the calculation result fits with the poloidal flux contour of the plasma boundary (as shown in Fig. 7). The results of the benchmark proved the validity of the relational expressions for the three strike point control modes. We then implemented the strike point control experiments based on these relational expressions using an isoflux control method. 4 Strike point control experiments We modified the strike point control module in the isoflux control algorithm on EAST. The flexibility of the strike point control was proven in the EAST experiment operation. Our experimental results are analyzed below. 777
5 Plasma Science and Technology, Vol.17, No.9, Sep (a) Top X point (b) Bottom X point Fig.6 The relative error between the rtefit results and the calculation results using fitting function of top (a) and bottom (b) X points position Fig.7 The estimated results of the X points (red star) by using fitting expressions in three strike point control mode: (a) Mode I, (b) Mode II, (c) Mode III. The red hexagrams on the target plates are strike points. The plasma last closed flux surface is illustrated in blue solid line 4.1 Reference experiment With the constant target value of the X point position, we can estimate the desired strike point position. From the experimental results shown in Fig. 9, we found that the strike points had oscillated around the desired value, and the margin of difference was 1%5%. At first, we took shot without controlling the strike point position as the reference experiment. In the reference shot, we turned on the X point control in the discharge process. The control effect of X point is shown in Fig. 8. Fig.8 The experimental X point position (black line) and the target value (red dashed line) during discharge process of shot Fig.9 The strike point position (black line) and the desired value (red dashed line) during discharge process of shot (without strike point control) 778
6 XING Zhe et al.: Strike Point Control on EAST Using an Isoflux Control Method 4.2 Inner strike point control experiment In shot 36768, we tested control mode I as the representation, controlling both upper and lower inner strike points. With the same discharge setting of shot 36744, we turned on the strike point control since 3.2 s (target Z of the inner lower strike point is 94 cm and target Z of the inner upper strike point is 93 cm). Referring to the experimental results of shot (shown in Fig. 10), we observed that both upper and lower strike points had been controlled to oscillate around the target position. Furthermore, in comparison to the reference experiment, the oscillation amplitude of shot is less than that of shot Consequently, the inner strike point control works well by using the modified algorithm. We achieved the purpose of controlling the strike points to stay around our target position on EAST. referring to the sample space shown in Fig. 5. On the other hand, we carry out the strike point control by controlling X point positions using the isoflux control method. The control precision of the strike points depends on the control effect of the X points. To find the cause of this problem, we analyzed the control effect of X points (shown in Fig. 12). In control mode I, the active control variables are the radial positions of X point and inner strike point positions. The desired vertical positions of X point are estimated as the red dashed line in Fig. 12. Fig.11 The strike point position (black line) and the target value (red dashed line) during discharge process of shot (with inner strike point control) Fig.10 The strike point position (black line) and the target value (red dashed line) during discharge process of shot (with inner strike point control) Considering the primary aim of strike point control, we should control the strike points to move over an appropriate area of the divertor target plates during discharge. In shot and shot 36770, we set timevarying target positions of inner strike points in the discharge process. Taking shot as an example, the target strike control points varied from 3.2 s to 4.5 s, with Z of the lower inner strike control point shifting from 94 cm to 89 cm and Z of the upper inner strike control point shifting from 93 cm to 88 cm. The experimental results of shot are shown in Fig. 11. From 3.2 s to 4.5 s, the inner strike points shifted on the divertor target plate when the target value had changed. After 4.5 s, the strike points oscillated slightly around a constant position, the amplitude was less than 1%. We also noticed that there was an obvious gap between the actual strike point position and the target position. The active control failed to move the strike point to the target position accurately after 4.5 s. On one hand, the relational expressions between X point position and strike point position are valid in this shot Fig.12 The experimental X point position (black line) and the desired value (red dashed line) during discharge process of shot (with inner strike point control) From 3.2 s to 4.5 s, the target R of X points remains constant (shown in Fig. 12) while the inner strike points are expected to shift on the divertor target plates. However, before 3.2 s, the actual R of X point was far from the target value. During 3.2 s to 4.5 s, it was controlled to approach the target position instead of staying around the target position as expected. Even when the target position had been reached after 4 s, the oscillation of actual R of X point was a little big. This 779
7 explains why the strike points could not reach the target position till 4.5 s, even though they were successfully controlled to shift in the right direction. As a result, the actual Z of X point could not agree with the desired value very well while the strike points shift. Similar inaccurate control effect of X points was also found in the reference shot (shown in Fig. 8). We need to improve the control effect of X point position by adjusting control parameters in the future. While the strike points move as programmed, we should control the X points to reach the target position. 5 Summary A reliable control of the strike point position is very important for advanced tokamak operation. We applied the strike point control algorithm using rtefit-isoflux control method on EAST. We fitted the relational expressions for the strike point control by a statistical technique that was based on historical experimental data, and then incorporated them into the plasma control system as a part of the plasma shape controller. Benchmark experiments on EAST were carried out to validate the correctness and robustness of the control methods. We achieved the aim of moving the strike points on the divertor target plates under control. The control precision of the strike point positions was limited to the control effect of the X point position. We will adjust the control parameters to improve the control effect of X point position in the future. Furthermore, we can in this way implement the study of plasma facing materials on divertor target plates more effectively. It is of great importance to better adjust the particle deposition and thermal load for future experiments on EAST. Plasma Science and Technology, Vol.17, No.9, Sep References 1 Wesson J. 1997, Tokamaks. Oxford University Press, Oxford, UK 2 Pironti A and Walker M. 2005, IEEE Control Systems Magazine, 25: 30 3 Beghi A and Cenedese A. 2005, IEEE Control Systems Magazine, 25: 44 4 Mirnov S, Wesley J, Fujisawa N, et al. 1999, Nuclear Fusion, 39: Ariola M and Pironti A. 2005, IEEE Control Systems Magazine, 25: 65 6 Kolemen E, Gates D A, Rowley C W, et al. 2010, Nuclear Fusion, 50: Petrie T W, Fenstermacher M E, Lasnier C J. 2001, Fusion Technology, 39: Ferron J R, Penaflor B, Walker M L, et al. 1996, Fusion Engineering, 2: Xiao Bingjia, Humphreys D A, Walker M L, et al. 2008, Fusion Engineering and Design, 83: Ferron J R, Walker M L, Lao L L, et al. 1998, Nuclear Fusion, 38: Lao L L, John H St, Stambaugh R D, et al. 1985, Nuclear Fusion, 25: Penaflor B G, Ferron J R, Walker M L, et al. 2008, Fusion Engineering and Design, 83: Walker M L, Humphreys D A, Ferron J R. 1997, Proceedings of the 36th IEEE Conference, 4: Ferron J R, Walker M L, Penaflor B, et al. 2003, Application Programmer s Guide to the Plasma Control System, DIII-D Documents. General Atomics, Gneral Atomics Ct., San Diego (Manuscript received 16 October 2014) (Manuscript accepted 26 December 2014) address of XING Zhe: xzpb@ipp.ac.cn 780
FX Smile Modelling. 9 September September 9, 2008
FX Smile Modelling 9 September 008 September 9, 008 Contents 1 FX Implied Volatility 1 Interpolation.1 Parametrisation............................. Pure Interpolation.......................... Abstract
More informationConvergence of statistical moments of particle density time series in scrape-off layer plasmas
Convergence of statistical moments of particle density time series in scrape-off layer plasmas R. Kube and O. E. Garcia Particle density fluctuations in the scrape-off layer of magnetically confined plasmas,
More informationSegmentation and Scattering of Fatigue Time Series Data by Kurtosis and Root Mean Square
Segmentation and Scattering of Fatigue Time Series Data by Kurtosis and Root Mean Square Z. M. NOPIAH 1, M. I. KHAIRIR AND S. ABDULLAH Department of Mechanical and Materials Engineering Universiti Kebangsaan
More informationEfficient Utilization Condition of MACD on Stock Market and Nontrend Status Detecting Indicatior
Indian Journal of Science and Technology, Vol 9(24), DOI: 1017485/ijst/2016/v9i24/96033, June 2016 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 Efficient Utilization Condition of MACD on Stock Market
More informationEstimating A Smooth Term Structure of Interest Rates
E STIMATING A SMOOTH LTA 2/98 TERM STRUCTURE P. 159 177 OF INTEREST RATES JARI KÄPPI 1 Estimating A Smooth Term Structure of Interest Rates ABSTRACT This paper extends the literature of the term structure
More informationPROBABILITY CONTENT OF ERROR ELLIPSE AND ERROR CONTOUR (navell-08.mcd)
PROBABILITY CONTENT OF ERROR ELLIPSE AND ERROR CONTOUR (navell-8.mcd) 6.. Conditions of Use, Disclaimer This document contains scientific work. I cannot exclude, that the algorithm or the calculations
More informationStatistical and Machine Learning Approach in Forex Prediction Based on Empirical Data
Statistical and Machine Learning Approach in Forex Prediction Based on Empirical Data Sitti Wetenriajeng Sidehabi Department of Electrical Engineering Politeknik ATI Makassar Makassar, Indonesia tenri616@gmail.com
More informationPrentice Hall Connected Mathematics 2, 7th Grade Units 2009 Correlated to: Minnesota K-12 Academic Standards in Mathematics, 9/2008 (Grade 7)
7.1.1.1 Know that every rational number can be written as the ratio of two integers or as a terminating or repeating decimal. Recognize that π is not rational, but that it can be approximated by rational
More informationEmpirical Study on Short-Term Prediction of Shanghai Composite Index Based on ARMA Model
Empirical Study on Short-Term Prediction of Shanghai Composite Index Based on ARMA Model Cai-xia Xiang 1, Ping Xiao 2* 1 (School of Hunan University of Humanities, Science and Technology, Hunan417000,
More informationANALYSIS OF THE BINOMIAL METHOD
ANALYSIS OF THE BINOMIAL METHOD School of Mathematics 2013 OUTLINE 1 CONVERGENCE AND ERRORS OUTLINE 1 CONVERGENCE AND ERRORS 2 EXOTIC OPTIONS American Options Computational Effort OUTLINE 1 CONVERGENCE
More informationComputational Finance Binomial Trees Analysis
Computational Finance Binomial Trees Analysis School of Mathematics 2018 Review - Binomial Trees Developed a multistep binomial lattice which will approximate the value of a European option Extended the
More informationStock Trading System Based on Formalized Technical Analysis and Ranking Technique
Stock Trading System Based on Formalized Technical Analysis and Ranking Technique Saulius Masteika and Rimvydas Simutis Faculty of Humanities, Vilnius University, Muitines 8, 4428 Kaunas, Lithuania saulius.masteika@vukhf.lt,
More information38050 Povo Trento (Italy), Via Sommarive 14
UNIVERSITY OF TRENTO DEPARTMENT OF INFORMATION AND COMMUNICATION TECHNOLOGY 38050 Povo Trento (Italy), Via Sommarive 14 http://www.dit.unitn.it STOCHASTIC AND REACTIVE METHODS FOR THE DETERMINATION OF
More informationSpline Methods for Extracting Interest Rate Curves from Coupon Bond Prices
Spline Methods for Extracting Interest Rate Curves from Coupon Bond Prices Daniel F. Waggoner Federal Reserve Bank of Atlanta Working Paper 97-0 November 997 Abstract: Cubic splines have long been used
More informationResale Price and Cost-Plus Methods: The Expected Arm s Length Space of Coefficients
International Alessio Rombolotti and Pietro Schipani* Resale Price and Cost-Plus Methods: The Expected Arm s Length Space of Coefficients In this article, the resale price and cost-plus methods are considered
More informationA MATHEMATICAL PROGRAMMING APPROACH TO ANALYZE THE ACTIVITY-BASED COSTING PRODUCT-MIX DECISION WITH CAPACITY EXPANSIONS
A MATHEMATICAL PROGRAMMING APPROACH TO ANALYZE THE ACTIVITY-BASED COSTING PRODUCT-MIX DECISION WITH CAPACITY EXPANSIONS Wen-Hsien Tsai and Thomas W. Lin ABSTRACT In recent years, Activity-Based Costing
More informationPARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS
PARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS Melfi Alrasheedi School of Business, King Faisal University, Saudi
More informationSAMPLE. HSC formula sheet. Sphere V = 4 πr. Volume. A area of base
Area of an annulus A = π(r 2 r 2 ) R radius of the outer circle r radius of the inner circle HSC formula sheet Area of an ellipse A = πab a length of the semi-major axis b length of the semi-minor axis
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationSmooth estimation of yield curves by Laguerre functions
Smooth estimation of yield curves by Laguerre functions A.S. Hurn 1, K.A. Lindsay 2 and V. Pavlov 1 1 School of Economics and Finance, Queensland University of Technology 2 Department of Mathematics, University
More informationOn the use of time step prediction
On the use of time step prediction CODE_BRIGHT TEAM Sebastià Olivella Contents 1 Introduction... 3 Convergence failure or large variations of unknowns... 3 Other aspects... 3 Model to use as test case...
More informationTraditional Optimization is Not Optimal for Leverage-Averse Investors
Posted SSRN 10/1/2013 Traditional Optimization is Not Optimal for Leverage-Averse Investors Bruce I. Jacobs and Kenneth N. Levy forthcoming The Journal of Portfolio Management, Winter 2014 Bruce I. Jacobs
More informationF A S C I C U L I M A T H E M A T I C I
F A S C I C U L I M A T H E M A T I C I Nr 38 27 Piotr P luciennik A MODIFIED CORRADO-MILLER IMPLIED VOLATILITY ESTIMATOR Abstract. The implied volatility, i.e. volatility calculated on the basis of option
More informationAvailable online at ScienceDirect. Procedia Economics and Finance 32 ( 2015 ) Andreea Ro oiu a, *
Available online at www.sciencedirect.com ScienceDirect Procedia Economics and Finance 32 ( 2015 ) 496 502 Emerging Markets Queries in Finance and Business Monetary policy and time varying parameter vector
More informationsinc functions with application to finance Ali Parsa 1*, J. Rashidinia 2
sinc functions with application to finance Ali Parsa 1*, J. Rashidinia 1 School of Mathematics, Iran University of Science and Technology, Narmak, Tehran 1684613114, Iran *Corresponding author: aliparsa@iust.ac.ir
More informationBacktesting Performance with a Simple Trading Strategy using Market Orders
Backtesting Performance with a Simple Trading Strategy using Market Orders Yuanda Chen Dec, 2016 Abstract In this article we show the backtesting result using LOB data for INTC and MSFT traded on NASDAQ
More informationOption Pricing Using Bayesian Neural Networks
Option Pricing Using Bayesian Neural Networks Michael Maio Pires, Tshilidzi Marwala School of Electrical and Information Engineering, University of the Witwatersrand, 2050, South Africa m.pires@ee.wits.ac.za,
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationApplication of Innovations Feedback Neural Networks in the Prediction of Ups and Downs Value of Stock Market *
Proceedings of the 6th World Congress on Intelligent Control and Automation, June - 3, 006, Dalian, China Application of Innovations Feedback Neural Networks in the Prediction of Ups and Downs Value of
More informationDynamic Interaction Network to Model the Interactive Patterns of International Stock Markets
World Academy of Science, Engineering and Technology 59 29 Dynamic Interaction Network to Model the Interactive Patterns of International Stock Markets Laura Lukmanto, Harya Widiputra, Lukas Abstract Studies
More informationAn Adjusted Trinomial Lattice for Pricing Arithmetic Average Based Asian Option
American Journal of Applied Mathematics 2018; 6(2): 28-33 http://www.sciencepublishinggroup.com/j/ajam doi: 10.11648/j.ajam.20180602.11 ISSN: 2330-0043 (Print); ISSN: 2330-006X (Online) An Adjusted Trinomial
More informationBinary Options Trading Strategies How to Become a Successful Trader?
Binary Options Trading Strategies or How to Become a Successful Trader? Brought to You by: 1. Successful Binary Options Trading Strategy Successful binary options traders approach the market with three
More informationThese notes essentially correspond to chapter 13 of the text.
These notes essentially correspond to chapter 13 of the text. 1 Oligopoly The key feature of the oligopoly (and to some extent, the monopolistically competitive market) market structure is that one rm
More informationDeep RL and Controls Homework 1 Spring 2017
10-703 Deep RL and Controls Homework 1 Spring 2017 February 1, 2017 Due February 17, 2017 Instructions You have 15 days from the release of the assignment until it is due. Refer to gradescope for the exact
More informationDiploma in Financial Management with Public Finance
Diploma in Financial Management with Public Finance Cohort: DFM/09/FT Jan Intake Examinations for 2009 Semester II MODULE: STATISTICS FOR FINANCE MODULE CODE: QUAN 1103 Duration: 2 Hours Reading time:
More informationQuadrant marked mesh patterns in 123-avoiding permutations
Quadrant marked mesh patterns in 23-avoiding permutations Dun Qiu Department of Mathematics University of California, San Diego La Jolla, CA 92093-02. USA duqiu@math.ucsd.edu Jeffrey Remmel Department
More informationThe Missing Knob on Your Random Vibration Controller Philip Van Baren, Vibration Research Corporation, Jenison, Michigan
The Missing Knob on Your Random Vibration Controller Philip Van Baren, Vibration Research Corporation, Jenison, Michigan Random vibration testing is the industry workhorse for simulating the environment
More informationIncome inequality and the growth of redistributive spending in the U.S. states: Is there a link?
Draft Version: May 27, 2017 Word Count: 3128 words. SUPPLEMENTARY ONLINE MATERIAL: Income inequality and the growth of redistributive spending in the U.S. states: Is there a link? Appendix 1 Bayesian posterior
More informationRisk-Adjusted Futures and Intermeeting Moves
issn 1936-5330 Risk-Adjusted Futures and Intermeeting Moves Brent Bundick Federal Reserve Bank of Kansas City First Version: October 2007 This Version: June 2008 RWP 07-08 Abstract Piazzesi and Swanson
More informationReinforcement Learning Analysis, Grid World Applications
Reinforcement Learning Analysis, Grid World Applications Kunal Sharma GTID: ksharma74, CS 4641 Machine Learning Abstract This paper explores two Markov decision process problems with varying state sizes.
More informationAn Improvement of Vegetation Height Estimation Using Multi-baseline Polarimetric Interferometric SAR Data
PIERS ONLINE, VOL. 5, NO. 1, 29 6 An Improvement of Vegetation Height Estimation Using Multi-baseline Polarimetric Interferometric SAR Data Y. S. Zhou 1,2,3, W. Hong 1,2, and F. Cao 1,2 1 National Key
More informationOn Stochastic Evaluation of S N Models. Based on Lifetime Distribution
Applied Mathematical Sciences, Vol. 8, 2014, no. 27, 1323-1331 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.412 On Stochastic Evaluation of S N Models Based on Lifetime Distribution
More information9. Logit and Probit Models For Dichotomous Data
Sociology 740 John Fox Lecture Notes 9. Logit and Probit Models For Dichotomous Data Copyright 2014 by John Fox Logit and Probit Models for Dichotomous Responses 1 1. Goals: I To show how models similar
More informationTechnical Summary. Contract for Cooling Water Plant Installation Works in Building 61
Contract for Cooling Water Plant Installation Works in Building 61 1 Purpose The purpose of this Call for Nomination is to establish a contract for early cooling water plant installation in the Site Services
More informationA Skewed Truncated Cauchy Uniform Distribution and Its Moments
Modern Applied Science; Vol. 0, No. 7; 206 ISSN 93-844 E-ISSN 93-852 Published by Canadian Center of Science and Education A Skewed Truncated Cauchy Uniform Distribution and Its Moments Zahra Nazemi Ashani,
More informationDynamic Linkages between Newly Developed Islamic Equity Style Indices
ISBN 978-93-86878-06-9 9th International Conference on Business, Management, Law and Education (BMLE-17) Kuala Lumpur (Malaysia) Dec. 14-15, 2017 Dynamic Linkages between Newly Developed Islamic Equity
More informationOPENING RANGE BREAKOUT STOCK TRADING ALGORITHMIC MODEL
OPENING RANGE BREAKOUT STOCK TRADING ALGORITHMIC MODEL Mrs.S.Mahalakshmi 1 and Mr.Vignesh P 2 1 Assistant Professor, Department of ISE, BMSIT&M, Bengaluru, India 2 Student,Department of ISE, BMSIT&M, Bengaluru,
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationSeasonal Pathloss Modeling at 900MHz for OMAN
2011 International Conference on Telecommunication Technology and Applications Proc.of CSIT vol.5 (2011) (2011) IACSIT Press, Singapore Seasonal Pathloss Modeling at 900MHz for OMAN Zia Nadir + Electrical
More informationSetting JET on Track. Prof. Donato Palumbo Honorary General Director Research EURATOM
Setting JET on Track Prof. Donato Palumbo Honorary General Director Research EURATOM 25th JET Anniversary, Culham 20th May 2004 Fusion Reaction: D + T 4 He + n T D Alpha particle : 20 % of reaction energy
More informationstarting on 5/1/1953 up until 2/1/2017.
An Actuary s Guide to Financial Applications: Examples with EViews By William Bourgeois An actuary is a business professional who uses statistics to determine and analyze risks for companies. In this guide,
More informationKERNEL PROBABILITY DENSITY ESTIMATION METHODS
5.- KERNEL PROBABILITY DENSITY ESTIMATION METHODS S. Towers State University of New York at Stony Brook Abstract Kernel Probability Density Estimation techniques are fast growing in popularity in the particle
More informationChapter DIFFERENTIAL EQUATIONS: PHASE SPACE, NUMERICAL SOLUTIONS
Chapter 10 10. DIFFERENTIAL EQUATIONS: PHASE SPACE, NUMERICAL SOLUTIONS Abstract Solving differential equations analytically is not always the easiest strategy or even possible. In these cases one may
More informationSolving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function?
DOI 0.007/s064-006-9073-z ORIGINAL PAPER Solving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function? Jules H. van Binsbergen Michael W. Brandt Received:
More informationState Switching in US Equity Index Returns based on SETAR Model with Kalman Filter Tracking
State Switching in US Equity Index Returns based on SETAR Model with Kalman Filter Tracking Timothy Little, Xiao-Ping Zhang Dept. of Electrical and Computer Engineering Ryerson University 350 Victoria
More informationOptimizing Modular Expansions in an Industrial Setting Using Real Options
Optimizing Modular Expansions in an Industrial Setting Using Real Options Abstract Matt Davison Yuri Lawryshyn Biyun Zhang The optimization of a modular expansion strategy, while extremely relevant in
More informationApril, 2006 Vol. 5, No. 4
April, 2006 Vol. 5, No. 4 Trading Seasonality: Tracking Market Tendencies There s more to seasonality than droughts and harvests. Find out how to make seasonality work in your technical toolbox. Issue:
More informationImpact of Weekdays on the Return Rate of Stock Price Index: Evidence from the Stock Exchange of Thailand
Journal of Finance and Accounting 2018; 6(1): 35-41 http://www.sciencepublishinggroup.com/j/jfa doi: 10.11648/j.jfa.20180601.15 ISSN: 2330-7331 (Print); ISSN: 2330-7323 (Online) Impact of Weekdays on the
More informationANALYSIS OF MACROECONOMIC FACTORS AFFECTING SHARE PRICE OF PT. BANK MANDIRI Tbk
ANALYSIS OF MACROECONOMIC FACTORS AFFECTING SHARE PRICE OF PT. BANK MANDIRI Tbk Camalia Zahra 1 Management Study Program, Faculty of Business, President University, Indonesia Camalia.zahra@gmail.com Purwanto
More informationAn IMEX-method for pricing options under Bates model using adaptive finite differences Rapport i Teknisk-vetenskapliga datorberäkningar
PROJEKTRAPPORT An IMEX-method for pricing options under Bates model using adaptive finite differences Arvid Westlund Rapport i Teknisk-vetenskapliga datorberäkningar Jan 2014 INSTITUTIONEN FÖR INFORMATIONSTEKNOLOGI
More informationApplied Mathematics Letters. On local regularization for an inverse problem of option pricing
Applied Mathematics Letters 24 (211) 1481 1485 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml On local regularization for an inverse
More informationThe mathematical model of portfolio optimal size (Tehran exchange market)
WALIA journal 3(S2): 58-62, 205 Available online at www.waliaj.com ISSN 026-386 205 WALIA The mathematical model of portfolio optimal size (Tehran exchange market) Farhad Savabi * Assistant Professor of
More informationImpact of Devaluation on Trade Balance in Pakistan
Page 16 Oeconomics of Knowledge, Volume 3, Issue 3, 3Q, Summer 2011 Impact of Devaluation on Trade Balance in Pakistan Muhammad ASIF, Lecturer Management Sciences Department CIIT, Abbottabad, Pakistan
More information1 SE = Student Edition - TG = Teacher s Guide
Mathematics State Goal 6: Number Sense Standard 6A Representations and Ordering Read, Write, and Represent Numbers 6.8.01 Read, write, and recognize equivalent representations of integer powers of 10.
More informationModeling, Analysis, and Characterization of Dubai Financial Market as a Social Network
Modeling, Analysis, and Characterization of Dubai Financial Market as a Social Network Ahmed El Toukhy 1, Maytham Safar 1, Khaled Mahdi 2 1 Computer Engineering Department, Kuwait University 2 Chemical
More informationInflation Regimes and Monetary Policy Surprises in the EU
Inflation Regimes and Monetary Policy Surprises in the EU Tatjana Dahlhaus Danilo Leiva-Leon November 7, VERY PRELIMINARY AND INCOMPLETE Abstract This paper assesses the effect of monetary policy during
More informationTHE ECONOMICS OF BANK ROBBERIES IN NEW ENGLAND 1. Kimberly A. Leonard, Diane L. Marley & Charlotte A. Senno
THE ECONOMICS OF BANK ROBBERIES IN NEW ENGLAND 1 The Economics of Bank Robberies in New England Kimberly A. Leonard, Diane L. Marley & Charlotte A. Senno The University of Rhode Island, STA308 Comment
More informationPower law in market capitalization Title and Shanghai bubble periods. Mizuno, Takayuki; Ohnishi, Takaaki; Author(s) Tsutomu
Power law in market capitalization Title and Shanghai bubble periods Mizuno, Takayuki; Ohnishi, Takaaki; Author(s) Tsutomu Citation Issue 2016-07 Date Type Technical Report Text Version publisher URL http://hdl.handle.net/10086/27965
More informationSTOCK PRICE PREDICTION: KOHONEN VERSUS BACKPROPAGATION
STOCK PRICE PREDICTION: KOHONEN VERSUS BACKPROPAGATION Alexey Zorin Technical University of Riga Decision Support Systems Group 1 Kalkyu Street, Riga LV-1658, phone: 371-7089530, LATVIA E-mail: alex@rulv
More informationCategorical. A general name for non-numerical data; the data is separated into categories of some kind.
Chapter 5 Categorical A general name for non-numerical data; the data is separated into categories of some kind. Nominal data Categorical data with no implied order. Eg. Eye colours, favourite TV show,
More informationModelling of selected S&P 500 share prices
MPRA Munich Personal RePEc Archive Modelling of selected S&P 5 share prices Ivan Kitov and Oleg Kitov IDG RAS 22. June 29 Online at http://mpra.ub.uni-muenchen.de/15862/ MPRA Paper No. 15862, posted 22.
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
More informationECONOMIC GROWTH AND UNEMPLOYMENT RATE OF THE TRANSITION COUNTRY THE CASE OF THE CZECH REPUBLIC
ECONOMIC GROWTH AND UNEMPLOMENT RATE OF THE TRANSITION COUNTR THE CASE OF THE CZECH REPUBLIC 1996-2009 EKONOMIE Elena Mielcová Introduction In early 1960 s, the economist Arthur Okun documented the negative
More informationA Fast Procedure for Transmission Loss Allocation of Power Systems by Using DC-OPF
Bangladesh J. Sci. Ind. Res. 42(3), 249-256, 2007 Introduction A Fast Procedure for Transmission Loss Allocation of Power Systems by Using DC-OPF M. H. Kabir Department of Computer Science and Telecommunication
More informationApproximate Variance-Stabilizing Transformations for Gene-Expression Microarray Data
Approximate Variance-Stabilizing Transformations for Gene-Expression Microarray Data David M. Rocke Department of Applied Science University of California, Davis Davis, CA 95616 dmrocke@ucdavis.edu Blythe
More informationGovernment spending in a model where debt effects output gap
MPRA Munich Personal RePEc Archive Government spending in a model where debt effects output gap Peter N Bell University of Victoria 12. April 2012 Online at http://mpra.ub.uni-muenchen.de/38347/ MPRA Paper
More informationA Comparative Study of Ensemble-based Forecasting Models for Stock Index Prediction
Association for Information Systems AIS Electronic Library (AISeL) MWAIS 206 Proceedings Midwest (MWAIS) Spring 5-9-206 A Comparative Study of Ensemble-based Forecasting Models for Stock Index Prediction
More informationPinpointing Entry & Exit Points
This tutorial was originally titled The Four Steps to 80% Day Trading Success and was recorded at the Online Trading Expo Pinpointing Entry & Exit Points Dr. John F. Clayburg Email - clayburg@pionet.net
More informationOPTIMIZATION STUDY OF RSI EXPERT SYSTEM BASED ON SHANGHAI SECURITIES MARKET
0 th February 013. Vol. 48 No. 005-013 JATIT & LLS. All rights reserved. ISSN: 199-8645 www.jatit.org E-ISSN: 1817-3195 OPTIMIZATION STUDY OF RSI EXPERT SYSTEM BASED ON SHANGHAI SECURITIES MARKET HUANG
More informationGT RETROFIT AND RE-DESIGN OF A FOUR-STAGE CENTRIFUGAL COMPRESSOR
Proceedings of ASME Turbo Expo 2015: Turbine Technical Conference and Exposition GT2015 June 15-19, 2015, Montreal GT2015-43760 RETROFIT AND RE-DESIGN OF A FOUR-STAGE CENTRIFUGAL COMPRESSOR Chao Chen SoftInWay,
More informationF.2 Factoring Trinomials
1 F.2 Factoring Trinomials In this section, we discuss factoring trinomials. We start with factoring quadratic trinomials of the form 2 + bbbb + cc, then quadratic trinomials of the form aa 2 + bbbb +
More informationTaxation and Efficiency : (a) : The Expenditure Function
Taxation and Efficiency : (a) : The Expenditure Function The expenditure function is a mathematical tool used to analyze the cost of living of a consumer. This function indicates how much it costs in dollars
More informationNotes on a California Perspective of the Dairy Margin Protection Program (DMPP)
Notes on a California Perspective of the Dairy Margin Protection Program (DMPP) Leslie J. Butler Department of Agricultural & Resource Economics University of California-Davis If I were a California dairy
More informationCOMMUNICATION FROM THE COMMISSION TO THE EUROPEAN PARLIAMENT AND THE COUNCIL EU CONTRIBUTION TO A REFORMED ITER PROJECT. {SWD(2017) 232 final}
EUROPEAN COMMISSION Brussels, 14.6.2017 COM(2017) 319 final COMMUNICATION FROM THE COMMISSION TO THE EUROPEAN PARLIAMENT AND THE COUNCIL EU CONTRIBUTION TO A REFORMED ITER PROJECT {SWD(2017) 232 final}
More informationMulti-factor Stock Selection Model Based on Kernel Support Vector Machine
Journal of Mathematics Research; Vol. 10, No. 5; October 2018 ISSN 1916-9795 E-ISSN 1916-9809 Published by Canadian Center of Science and Education Multi-factor Stock Selection Model Based on Kernel Support
More informationDevelopment on Methods for Evaluating Structure Reliability of Piping Components
Transactions of the 17 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 17) Prague, Czech Republic, August 17, 00 Paper # M0- Development on Methods for Evaluating Structure
More informationFlight Ascend Values Analyzer User Guide. Introduction. Values Home Page. Navigation around the product
Flight Ascend Values Analyzer User Guide Introduction Flight Ascend Values Analyzer (FAVA) is a new platform designed to replace ASO Values and become the unified platform for accessing all FlightGlobal
More informationNon-linear logit models for high frequency currency exchange data
Non-linear logit models for high frequency currency exchange data N. Sazuka 1 & T. Ohira 2 1 Department of Physics, Tokyo Institute of Technology, Japan 2 Sony Computer Science Laboratories, Japan Abstract
More informationDynamic Optimization of Open-loop Input Signals for Ramp-up Current Profiles in Tokamak Plasmas 1
25 Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4. International license http://creativecommons.org/licenses/by-nc-nd/4./ Dynamic Optimization of Open-loop Input
More informationCopula-Based Pairs Trading Strategy
Copula-Based Pairs Trading Strategy Wenjun Xie and Yuan Wu Division of Banking and Finance, Nanyang Business School, Nanyang Technological University, Singapore ABSTRACT Pairs trading is a technique that
More informationThe Cost of Poor Reliability for Future Particle Accelerator Driven Systems in the UK Energy Market
The Cost of Poor Reliability for Future Particle Accelerator Driven Systems in the UK Energy Market Steven J Steer, William J Nuttall, Geoff T Parks & Leonardo V N Gonçalves Accelerator Driven Nuclear
More informationAn Empirical Study about Catering Theory of Dividends: The Proof from Chinese Stock Market
Journal of Industrial Engineering and Management JIEM, 2014 7(2): 506-517 Online ISSN: 2013-0953 Print ISSN: 2013-8423 http://dx.doi.org/10.3926/jiem.1013 An Empirical Study about Catering Theory of Dividends:
More informationTHE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE
THE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE GÜNTER ROTE Abstract. A salesperson wants to visit each of n objects that move on a line at given constant speeds in the shortest possible time,
More informationLong-Term Monitoring of Low-Volume Road Performance in Ontario
Long-Term Monitoring of Low-Volume Road Performance in Ontario Li Ningyuan, P. Eng. Tom Kazmierowski, P.Eng. Becca Lane, P. Eng. Ministry of Transportation of Ontario 121 Wilson Avenue Downsview, Ontario
More informationEffect of Mechanical Mismatch on the Stress Intensity Factors of Inclined Cracks under Mode I Tension Loading
Effect of Mechanical Mismatch on the Stress Intensity Factors of Inclined Cracks under Mode I Tension Loading M.A Azlan 1,a and A.E Ismail 1,b 1 Faculty of Mechanical and Manufacturing Engineering, Universiti
More informationChapter 14 : Statistical Inference 1. Note : Here the 4-th and 5-th editions of the text have different chapters, but the material is the same.
Chapter 14 : Statistical Inference 1 Chapter 14 : Introduction to Statistical Inference Note : Here the 4-th and 5-th editions of the text have different chapters, but the material is the same. Data x
More informationIs a Binomial Process Bayesian?
Is a Binomial Process Bayesian? Robert L. Andrews, Virginia Commonwealth University Department of Management, Richmond, VA. 23284-4000 804-828-7101, rlandrew@vcu.edu Jonathan A. Andrews, United States
More informationAn Empirical Analysis of Effect on Copper Futures Yield. Based on GARCH
An Empirical Analysis of Effect on Copper Futures Yield Based on GARCH Feng Li 1, Ping Xiao 2 * 1 (School of Hunan University of Humanities, Science and Technology, Hunan 417000, China) 2 (School of Hunan
More informationA Novel Iron Loss Reduction Technique for Distribution Transformers Based on a Combined Genetic Algorithm Neural Network Approach
16 IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS PART C: APPLICATIONS AND REVIEWS, VOL. 31, NO. 1, FEBRUARY 2001 A Novel Iron Loss Reduction Technique for Distribution Transformers Based on a Combined
More informationScholars Journal of Arts, Humanities and Social Sciences
DOI: 10.21276/sjahss.2016.4.6.11 Scholars Journal of Arts, Humanities and Social Sciences Sch. J. Arts Humanit. Soc. Sci. 2016; 4(6B):686-699 Scholars Academic and Scientific Publishers (SAS Publishers)
More information