A Stochastic Framework for Saint-Venant Torsion in Spherical Shells: Monte Carlo Implementation of the Feynman–Kac Approach

Abstract

This research introduces an innovative probabilistic method for examining torsional stress behavior in spherical shell structures through Monte Carlo simulation techniques. The spherical geometry of these components creates distinctive computational difficulties for conventional analytical and deterministic numerical approaches when solving torsion-related problems. The authors develop a comprehensive mesh-free Monte Carlo framework built upon the Feynman–Kac formula, which maintains the geometric symmetry of the domain while offering a probabilistic solution representation via stochastic processes on spherical surfaces. The technique models Brownian motion paths on spherical surfaces using the Euler–Maruyama numerical scheme, converting the Saint-Venant torsion equation into a problem of stochastic integration. The computational implementation utilizes the Fibonacci sphere technique for achieving uniform point placement, employs adaptive time-stepping strategies to address pole singularities, and incorporates efficient algorithms for boundary identification. This symmetry-maintaining approach circumvents the mesh generation complications inherent in finite element and finite difference techniques, which typically compromise the problem’s natural symmetry, while delivering comparable precision. Performance evaluations reveal nearly linear parallel computational scaling across up to eight processing cores with efficiency rates above 70%, making the method well-suited for multi-core computational platforms. The approach demonstrates particular effectiveness in analyzing torsional stress patterns in thin-walled spherical components under both symmetric and asymmetric boundary scenarios, where traditional grid-based methods encounter discretization and convergence difficulties. The findings offer valuable practical recommendations for material specification and structural design enhancement, especially relevant for pressure vessel and dome structure applications experiencing torsional loads. However, the probabilistic characteristics of the method create statistical uncertainty that requires cautious result interpretation, and computational expenses may surpass those of deterministic approaches for less complex geometries. Engineering analysis of the outcomes provides actionable recommendations for optimizing material utilization and maintaining structural reliability under torsional loading conditions.