Fractional Brownian motion (fBM) was first introduced within a Hilbert space framework by Kolmogorov [Kol40], while studying spiral curves in Hilbert space; the process was further studied and was coined the name ’fractional Brownian motion’ in the 1968 paper by Mandelbrot and Van Ness [MVN68]. It has been widely used in various scientific fields, most notability in hydrology as first suggested in [Man65]. It also plays an important role in communication technology by enriching the queuing theory in terms of simulating real network traffic. In recent years, it has been steadily gaining traction in finance. This is due to the fact that, traditional stochastic volatility model driven by ordinary Brownian motion implies exponential decay of the implied volatility smile, where empirical studies shows the decay is of power order. This can be explained by the long-memoryness found in autocovariance of the time-series, and this long-memory feature in instantaneous volatility process is called the volatility persistence. Even though this phenomenon is commonly observed on market implied volatility surface, it has been largely ignored because of the difficulty to capture it within the ordinary stochastic volatility framework. Within the framework of time-series analysis, the volatility persistence is concluded as the long-range dependence displaying in the instantaneous volatility time-series, or depicted by prominent components at low frequency under spectral density of the autocovariance function. This phenomenon is so commonly found, that Granger [Gra66] considered such phenomenon as the ”typical spectral shape of an economic variable”. The presence of the long-memory phenomenon has important implication in financial economics and financial engineering, especially in the area of the portfolio optimization: the choice of optimal consumption/saving because now the optimal decision might become very sensitive to the investment horizon instead of an asymptotically time-homogeneous problem if the asset return now is long-range dependent, also in the area of derivative pricing. The introduction of long-memory process is inconsistent with the pre-existing continuous-time process framework commonly utilized in derivative pricing (see, [Mah90b], [Mah90a]), and in [Mer87], [LeR89], Merton and Leroy respectively discussed the relationship of the efficient market hypothesis argument and the presence of the long-memory phenomenon. For detail discussion about the statistical test for the existence of long-range memory, see [Lo91] (the author has also proposed a robust extension of the R/S statistics for the purpose of robust statistical inference). Abudant empirical studies have been done with real data, most exmplemary the study done by Greene and Fielitz [GF77]. The studies was done on securities listed on the New York Stock Exchange, and many were found displaying long-range dependence in their daily returns. Several modeling approaches have been suggested capturing this persistence in conditional variance either via a unit-root or long memory process. In order to keep the pricing-framework largely intact, it is more interesting to study the long memory process, and fBM offers a simplistic extension particular for this purpose, owning to its similarity to the ordinary Brownian motion and its Gaussian properties, on the familiar process. In this paper, several approaches to simulate fractional Brownian motion with H > 1=2 are outlined, including the exact methods and approximate methods, where the Hurst Index H is a parameter used in literature to generalize Brownian motion into fractional Brownian motion, first made popular by Benoit Mandelbrot. Brief introduction of the truncated fractional Brownian motion (long-memory model in continuous time) is also included, as proposed in [CR96], [CR98], which is shown to be inadequate to replicate the fractional Brownian motion. Through the full Monte-Carlo simulation scheme, the implied volatility surface is constructed. One of the main result in the research is that, imposing correlation between the fractional Brownian motion driven stochastic volatility and the ordinary Brownian motion driven asset process does not translates into skewness of the implied volatility surface, this is further supported by E.Alos’s paper [ALV07] on the observation of the close-to-maturity implied volatility surface. Unfortunately, since explicit pricing of the European option under the fractionally-driven stochastic volatility is not available in closed-form due to the non-Markovian nature of the fractional Brownian Motion driven stochastic volatility, and market participants have to rely on computational intensive method such as Monte-Carlo simulation, such dilemma serves as our motivation to explore for an robust approximation. The simulation is built on top of the work of Funahashi [Fun12], which provides a closed-form approximation of European option under the stochastic local-volatility model, this serves as a starting point of our robust simulation approach, reducing the brute-force simulation dimension from three to just one, rendering it a computationally inexpensive pricing scheme. As mentioned before, the skewness cannot be modeled by simply imposing correlation between the fractional Brownian motion and ordinary Brownian motion. To have a significant skewness under our proposed framework, it is necessary to add an ordinary Brownian motion driver in the stochastic volatility as well, in order to establish a correlation between the asset process and the stochastic volatility, thus the correlation. Our approximation-basedsimulation scheme is also capable of pricing European option under this rather complicated stochastic environment, which captures the skewness, smile and persistence of the implied volatility surface. The paper is structured as following: In chapter 1 the motivation is provided in financial context, and a brief description of different common approaches pertaining to the particular problem at hand. In chapter 2, the background information and technical definition of the fractional Brownian motion are outlined. Chapter 3 depicts various financial modeling involves fractional Brownian motion, such as the asset process model, various stochastic volatility models, as well as discussion of presence of arbitrage in the presence of fractional Brownian motion. Chapter 4 outlines various approach to simulate the exact fractional Brownian motion, which we also provide the Fast-Fourier-Transform (FFT) based simulation in detail, as it is our choice of simulation tool. We also go into detail of some approximate simulation approach in chapter 5 for the sake of completion, and comparison, providing the reader choices to of simulation scheme, for example, the spectral method is a great substitution for FFT-based simulation, and the correlated random walk draws an analogy with construction of ordinary Brownian motion by independent random walks. Chapter 6 gives a numerical example of the Comte, Renault approach and point out the inadequacy of this widely adapted approach. Chapter 7 gives us the full simulation scheme for the asset process with fractionally driven stochastic volatility, while exploration the effect of various factors such as the correlation, fractional Brownian motion, and the vol-of-vol for both ordinary Brownian motion and fractional Brownian motion, arguing the necessity for a fully extended hybrid fBM model which has both the correlated ordinary Brownian motion (with the asset return), and uncorrelated fractional Brownian motion to fully captures the stylized phenomenon observed on the market. Chapter 8 depicts the Funahashi approximation scheme and its derivation. Chapter 9 outlines the extension of the Funahashi approximation scheme for our robust simulation scheme which is capable of pricing the fully extended mixture fBM model, and various results depicting the stylized features we aim to capture. In chapter 10, calibration scheme against the market data utilizing the previous pricing scheme is provided, as well as exploring the relationship between parameters. Chapter 11 concludes the finding in this research.