American Institute of Mathematical Sciences

May  2020, 16(3): 1349-1368. doi: 10.3934/jimo.2019006

Pricing options on investment project contraction and ownership transfer using a finite volume scheme and an interior penalty method

 1 Department of Mathematics & Statistics, Curtin University, GPO Box U1987, WA 6845, Australia 2 School of Mathematical & Software Sciences, Sichuan Normal University, Sichuan, China, 610000 3 Coordinated Innovation Center for Computable Modeling in Management Science, Tianjin University of Finance and Economics Tianjin 300222, China

Received  February 2018 Revised  September 2018 Published  March 2019

In this work we develop partial differential equation (PDE) based computational models for pricing real options to contract the production or to transfer part/all of the ownership of a project when the underlying asset price of the project satisfies a geometric Brownian motion. The developed models are similar to the Black-Scholes equation for valuing conventional European put options or the partial differential linear complementarity problem (LCP) for pricing American put options. A finite volume method is used for the discretization of the PDE models and a penalty approach is applied to the discretized LCP. We show that the coefficient matrix of the discretized systems is a positive-definite $M$-matrix which guarantees that the solution from the penalty equation converges to that of the discretized LCP. Numerical experiments, performed to demonstrate the usefulness of our methods, show that our models and numerical methods are able to produce financially meaningful numerical results for the two non-trivial test problems.

Citation: Nan Li, Song Wang, Shuhua Zhang. Pricing options on investment project contraction and ownership transfer using a finite volume scheme and an interior penalty method. Journal of Industrial & Management Optimization, 2020, 16 (3) : 1349-1368. doi: 10.3934/jimo.2019006
References:

show all references

References:
The values of European and American contracting options when $\kappa = 0.5$ for Test 1
The difference between the values of the European and American contracting options for Test 1
The value of the option to abandon for Test 2
Computed $W - W^*$ and Greeks of $W$ when $\lambda = 0.5$ for Test 2
The American option value and its optimal exercise curve for $\lambda = 0.5$ for Test 2
Project, production and market data and functions used in Test 1
 $Q = 10^4$ million tons $B = 30\%$ per annum $c_0 = {\rm{US}}$ ＄25 $c(t) = c_0\times e^{0.005t}$ $R = 5\%$ per annum $r = 0.06$ per annum $C = {\rm{US}}$ ＄ $5 \times 10^2$ million $T = 1$ year $\sigma = 30\%$ $\delta = 0.02$ {{q}_{0}}=\left\{ \begin{array}{*{35}{l}} \begin{align} &0.01Q\times {{e}^{0.007t}},\ \ t\le T_{1}^{*}, \\ &0,\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ t\in (T_{1}^{*},{{T}^{*}}) \\ \end{align} \\ \end{array} \right. {{q}_{1}}=\left\{ \begin{align} &0.01Q\times {{e}^{0.007t}},\ \ \ \ \ t Q = 10^4$million tons$B = 30\%$per annum$c_0 = {\rm{US}} $＄25$c(t) = c_0\times e^{0.005t}  R = 5\%$per annum$r = 0.06 $per annum$C = {\rm{US}}$＄$5 \times 10^2$million$T = 1$year$\sigma = 30\% \delta = 0.02{{q}_{0}}=\left\{ \begin{array}{*{35}{l}} \begin{align} &0.01Q\times {{e}^{0.007t}},\ \ t\le T_{1}^{*}, \\ &0,\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ t\in (T_{1}^{*},{{T}^{*}}) \\ \end{align} \\ \end{array} \right.{{q}_{1}}=\left\{ \begin{align} &0.01Q\times {{e}^{0.007t}},\ \ \ \ \ t
 [1] Armin Lechleiter, Tobias Rienmüller. Factorization method for the inverse Stokes problem. Inverse Problems & Imaging, 2013, 7 (4) : 1271-1293. doi: 10.3934/ipi.2013.7.1271 [2] Xingchun Wang, Yongjin Wang. Variance-optimal hedging for target volatility options. Journal of Industrial & Management Optimization, 2014, 10 (1) : 207-218. doi: 10.3934/jimo.2014.10.207 [3] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [4] Qiang Guo, Dong Liang. An adaptive wavelet method and its analysis for parabolic equations. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 327-345. doi: 10.3934/naco.2013.3.327 [5] Tao Wu, Yu Lei, Jiao Shi, Maoguo Gong. An evolutionary multiobjective method for low-rank and sparse matrix decomposition. Big Data & Information Analytics, 2017, 2 (1) : 23-37. doi: 10.3934/bdia.2017006 [6] Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems & Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017 [7] Boris Kramer, John R. Singler. A POD projection method for large-scale algebraic Riccati equations. Numerical Algebra, Control & Optimization, 2016, 6 (4) : 413-435. doi: 10.3934/naco.2016018 [8] Petra Csomós, Hermann Mena. Fourier-splitting method for solving hyperbolic LQR problems. Numerical Algebra, Control & Optimization, 2018, 8 (1) : 17-46. doi: 10.3934/naco.2018002 [9] Christina Surulescu, Nicolae Surulescu. Modeling and simulation of some cell dispersion problems by a nonparametric method. Mathematical Biosciences & Engineering, 2011, 8 (2) : 263-277. doi: 10.3934/mbe.2011.8.263 [10] Min Li. A three term Polak-Ribière-Polyak conjugate gradient method close to the memoryless BFGS quasi-Newton method. Journal of Industrial & Management Optimization, 2020, 16 (1) : 245-260. doi: 10.3934/jimo.2018149 [11] Frank Sottile. The special Schubert calculus is real. Electronic Research Announcements, 1999, 5: 35-39. [12] Manfred Einsiedler, Elon Lindenstrauss. On measures invariant under diagonalizable actions: the Rank-One case and the general Low-Entropy method. Journal of Modern Dynamics, 2008, 2 (1) : 83-128. doi: 10.3934/jmd.2008.2.83 [13] Xiaomao Deng, Xiao-Chuan Cai, Jun Zou. A parallel space-time domain decomposition method for unsteady source inversion problems. Inverse Problems & Imaging, 2015, 9 (4) : 1069-1091. doi: 10.3934/ipi.2015.9.1069 [14] Arseny Egorov. Morse coding for a Fuchsian group of finite covolume. Journal of Modern Dynamics, 2009, 3 (4) : 637-646. doi: 10.3934/jmd.2009.3.637 [15] Olena Naboka. On synchronization of oscillations of two coupled Berger plates with nonlinear interior damping. Communications on Pure & Applied Analysis, 2009, 8 (6) : 1933-1956. doi: 10.3934/cpaa.2009.8.1933 [16] Alexandr Mikhaylov, Victor Mikhaylov. Dynamic inverse problem for Jacobi matrices. Inverse Problems & Imaging, 2019, 13 (3) : 431-447. doi: 10.3934/ipi.2019021 [17] Hakan Özadam, Ferruh Özbudak. A note on negacyclic and cyclic codes of length $p^s$ over a finite field of characteristic $p$. Advances in Mathematics of Communications, 2009, 3 (3) : 265-271. doi: 10.3934/amc.2009.3.265 [18] Hyeong-Ohk Bae, Hyoungsuk So, Yeonghun Youn. Interior regularity to the steady incompressible shear thinning fluids with non-Standard growth. Networks & Heterogeneous Media, 2018, 13 (3) : 479-491. doi: 10.3934/nhm.2018021 [19] Diana Keller. Optimal control of a linear stochastic Schrödinger equation. Conference Publications, 2013, 2013 (special) : 437-446. doi: 10.3934/proc.2013.2013.437 [20] Guillaume Bal, Wenjia Jing. Homogenization and corrector theory for linear transport in random media. Discrete & Continuous Dynamical Systems - A, 2010, 28 (4) : 1311-1343. doi: 10.3934/dcds.2010.28.1311

2019 Impact Factor: 1.366