A Computational Study for Pricing European- and American-Type Options Under Heston’s Stochastic Volatility Model: Application of the SUPG-YZ $$\beta$$ β Formulation

Süleyman Cengizci () and Ömür Uğur ()
Additional contact information
Süleyman Cengizci: Antalya Bilim University
Ömür Uğur: Middle East Technical University

Computational Economics, 2025, vol. 66, issue 1, No 6, 179-206

Abstract: Abstract The interest of this paper is stabilized finite element approximations for pricing European- and American-type options under Heston’s stochastic volatility model, a generalization of the eminent Black–Scholes–Merton (BSM) framework in which volatility is treated as a constant. For spatial discretizations, the streamline-upwind/Petrov–Galerkin (SUPG) stabilized finite element method is used. The stabilized formulation is also supplemented with a shock-capturing mechanism, the so-called YZ $$\beta$$ β technique, in order to resolve localized sharp layers. The semi-discrete problems, i.e., the systems of time-dependent ordinary differential equations, are discretized in time with the Crank–Nicolson (CN) time-integration scheme. The resulting nonlinear algebraic equation systems are solved with the Newton–Raphson (NR) iterative process. The stabilized bi-conjugate gradient method, preconditioned with the incomplete lower–upper factorization technique, is employed for solving linearized systems. The linear complementarity problems arising in simulating American-type options are handled with an efficient and practical penalty approach, which comes at the cost of introducing a nonlinear source term to the fully discretized formulation. The in-house-developed solvers are verified first for the Heston model with a manufactured solution. Following that, the performances of the proposed method and techniques are evaluated on various test problems, including the digital options, through comparisons with other reported results. In addition to those studied previously, we also introduce new “challenging” parameter sets through which Heston’s model becomes much more convection-dominated and demonstrate the robustness of the proposed formulation and techniques for such cases. Furthermore, for each test case, the results obtained with the classical Galerkin finite element method and SUPG alone without shock-capturing are also presented, revealing that the SUPG-YZ $$\beta$$ β does not degrade the accuracy by introducing excessive numerical dissipation.

Keywords: Heston’s model; Option pricing; Stochastic volatility; SUPG; Shock-capturing (search for similar items in EconPapers)
Date: 2025
References: Add references at CitEc
Citations:

Downloads: (external link)
http://link.springer.com/10.1007/s10614-024-10704-3 Abstract (text/html)
Access to the full text of the articles in this series is restricted.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:kap:compec:v:66:y:2025:i:1:d:10.1007_s10614-024-10704-3

Ordering information: This journal article can be ordered from
http://www.springer. ... ry/journal/10614/PS2

DOI: 10.1007/s10614-024-10704-3

Access Statistics for this article

Computational Economics is currently edited by Hans Amman

More articles in Computational Economics from Springer, Society for Computational Economics Contact information at EDIRC.
Bibliographic data for series maintained by Sonal Shukla () and Springer Nature Abstracting and Indexing ().