Considering the fractional Brownian motion \(\textbf{B}_t\) as a generalized stochastic process in \(\textbf{S}'(\mathbb {R})\otimes (\textbf{S})_{-1}\) , the fractional white noise \(\textbf{W}^H_t\) is defined as the distributional derivative of \(\textbf{dB}^H_t\) . Using the framework of white noise analysis and the integral representation of \(\textbf{B}^H_t\) , we explicitly calculate the coefficients of their chaos expansion and provide a recurrence formula for their effective calculation. As a novel stochastic model, we introduce the notion of fractional Brownian motion and fractional white noise with a distributed-order Hurst parameter \(\textbf{H}\in (0,1)\) . Building on this construction, we derive numerical simulations of \(\textbf{B}^H_t\) and \(\textbf{W}^H_t\) by truncating the obtained chaos expansions, which allows us to approximate their sample paths and estimate the truncation error. We illustrate the results with two examples from finance and life insurance: a stock price model and a mortality model, both driven by a fractional white noise process with distributed-order Hurst parameter.