Avetik Arakelyan,Rafayel Barkhudaryan,Henrik Shahgholian,Mohammad M. Salehi
Abstract:In this paper we continue to study a non-local free boundary problem arising in financial bubbles. We focus on the parabolic counterpart of the bubble problem and suggest an iterative algorithm which consists of a sequence of parabolic obstacle problems at each step to be solved, that in turn gives the next obstacle function in the iteration. The convergence of the proposed algorithm is proved. Moreover, we consider the finite difference scheme for this algorithm and obtain its convergence. At the end of the paper we present and discuss computational results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the numerical treatment method for a class of non - local parabolic free - boundary problems in the context of financial bubbles. Specifically, the author focuses on how to analyze time - dependent (i.e., parabolic) financial bubble problems from the perspective of partial differential equations (PDEs). An iterative algorithm is proposed in the paper, which gradually approximates the final solution by solving a series of parabolic obstacle problems at each step, and the convergence of the proposed algorithm is proved. In addition, the author also considers the finite - difference scheme of this algorithm and obtains its convergence. Finally, the paper presents the computational results and discusses them.
### Mathematical Model
The model equation in the paper is a free - boundary problem, which can be expressed as a Hamilton - Jacobi equation:
\[
\min\left( Lu, u(t,x) - u(t, -x) - \psi(t,x) \right) = 0, \quad (t,x) \in \mathbb{R}^+ \times \Omega,
\]
where \(\Omega \subset \mathbb{R}\) is a symmetric bounded domain, satisfying that if \(x\in\Omega\) then \(-x\in\Omega\), and \(\psi\in C^2(\mathbb{R}^+ \times \Omega)\).
### Iterative Algorithm
The author proposes an iterative algorithm, which gradually approximates the final solution by solving a series of parabolic obstacle problems at each step. The specific steps are as follows:
1. **Initial Guess**: Define the initial function \(u_0\) as the solution to the following problem:
\[
\begin{cases}
Lu_0 = 0, & (t,x) \in \mathbb{R}^+ \times \Omega, \\
u_0(0,x) = g(0,x), & x\in\Omega, \\
u_0(t,x) = g(t,x), & (t,x) \in \mathbb{R}^+ \times \partial\Omega.
\end{cases}
\]
2. **Iterative Step**: Recursively define the sequence \(\{u_k\}\) as:
\[
\begin{cases}
\min\left( Lu_{k + 1}, u_{k + 1} - \tilde{u}_k - \psi \right) = 0, & (t,x) \in \mathbb{R}^+ \times \Omega, \\
u_{k + 1}(0,x) = g(0,x), & x\in\Omega, \\
u_{k + 1}(t,x) = g(t,x), & (t,x) \in \mathbb{R}^+ \times \partial\Omega.
\end{cases}
\]
where the obstacle problem corresponding to each \(k\) is \(\tilde{u}_k+\psi\).
### Convergence
The paper proves that the sequence \(\{u_k\}\) generated by this iterative algorithm is non - decreasing and converges to the viscosity solution of problem (1).
### Finite - Difference Scheme
To numerically solve each obstacle problem in the above - mentioned iterative algorithm, the author uses a finite - difference scheme. Specifically, for each iteration step, a nonlinear system needs to be solved:
\[
\begin{cases}
\min\left( L_h u_h^{(k)}, u_h^{(k)} - u^{(k - 1)}(t, -x) - \psi(t,x) \right) = 0, & (t,x) \in \Omega_T^h, \\
u_h^{(k)}(0,x) = g(0,x), & x\in\Omega_h, \\
u_h^{(k)}(n\Delta t, \pm a) = g(n\D