# 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] Xiaoni Chi, Zhongping Wan, Zijun Hao. A full-modified-Newton step $O(n)$ infeasible interior-point method for the special weighted linear complementarity problem. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021082 [2] Amanda E. Diegel. A C0 interior penalty method for the Cahn-Hilliard equation. Electronic Research Archive, , () : -. doi: 10.3934/era.2021030 [3] Hongsong Feng, Shan Zhao. A multigrid based finite difference method for solving parabolic interface problem. Electronic Research Archive, , () : -. doi: 10.3934/era.2021031 [4] 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 [5] Xiaofei Liu, Yong Wang. Weakening convergence conditions of a potential reduction method for tensor complementarity problems. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021080 [6] Jie-Wen He, Chi-Chon Lei, Chen-Yang Shi, Seak-Weng Vong. An inexact alternating direction method of multipliers for a kind of nonlinear complementarity problems. Numerical Algebra, Control & Optimization, 2021, 11 (3) : 353-362. doi: 10.3934/naco.2020030 [7] 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 [8] Li Chu, Bo Wang, Jie Zhang, Hong-Wei Zhang. Convergence analysis of a smoothing SAA method for a stochastic mathematical program with second-order cone complementarity constraints. Journal of Industrial & Management Optimization, 2021, 17 (4) : 1863-1886. doi: 10.3934/jimo.2020050 [9] Wan-Hua He, Chufang Wu, Jia-Wen Gu, Wai-Ki Ching, Chi-Wing Wong. Pricing vulnerable options under a jump-diffusion model with fast mean-reverting stochastic volatility. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021057 [10] Lunji Song, Wenya Qi, Kaifang Liu, Qingxian Gu. A new over-penalized weak galerkin finite element method. Part Ⅱ: Elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2581-2598. doi: 10.3934/dcdsb.2020196 [11] Yueqiang Shang, Qihui Zhang. A subgrid stabilizing postprocessed mixed finite element method for the time-dependent Navier-Stokes equations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 3119-3142. doi: 10.3934/dcdsb.2020222 [12] Derrick Jones, Xu Zhang. A conforming-nonconforming mixed immersed finite element method for unsteady Stokes equations with moving interfaces. Electronic Research Archive, , () : -. doi: 10.3934/era.2021032 [13] Claudia Lederman, Noemi Wolanski. An optimization problem with volume constraint for an inhomogeneous operator with nonstandard growth. Discrete & Continuous Dynamical Systems, 2021, 41 (6) : 2907-2946. doi: 10.3934/dcds.2020391 [14] Marita Holtmannspötter, Arnd Rösch, Boris Vexler. A priori error estimates for the space-time finite element discretization of an optimal control problem governed by a coupled linear PDE-ODE system. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021014 [15] 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 [16] Mohsen Abdolhosseinzadeh, Mir Mohammad Alipour. Design of experiment for tuning parameters of an ant colony optimization method for the constrained shortest Hamiltonian path problem in the grid networks. Numerical Algebra, Control & Optimization, 2021, 11 (2) : 321-332. doi: 10.3934/naco.2020028 [17] Sergei Avdonin, Julian Edward. An inverse problem for quantum trees with observations at interior vertices. Networks & Heterogeneous Media, 2021, 16 (2) : 317-339. doi: 10.3934/nhm.2021008 [18] Fioralba Cakoni, Pu-Zhao Kow, Jenn-Nan Wang. The interior transmission eigenvalue problem for elastic waves in media with obstacles. Inverse Problems & Imaging, 2021, 15 (3) : 445-474. doi: 10.3934/ipi.2020075 [19] Haodong Chen, Hongchun Sun, Yiju Wang. A complementarity model and algorithm for direct multi-commodity flow supply chain network equilibrium problem. Journal of Industrial & Management Optimization, 2021, 17 (4) : 2217-2242. doi: 10.3934/jimo.2020066 [20] 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

2019 Impact Factor: 1.366