Recent advances in Schramm-Loewner evolution have driven increasing interest in non-standard Loewner flows. In this work, we propose a novel splitting algorithm to simulate random Loewner curves with rigorous convergence analysis in sup-norm and $L^p$. The algorithm is further extended to explore fractional SLE, driven by fractional Brownian motion, and noise-reinforced SLE, incorporating the effect on long-term memory. These exploratory and numerical extensions enable theoretical predictions on fractal dimensions and other statistical phenomena, providing new insights into such dynamics and opening directions for future research.