Numerical meshfree path integration method for non-linear dynamic systems

We present a new numerical meshfree path integration (MPI) method for non-linear dynamic systems. The obtained MPI method can be performed in the irregular computational domain and the probability density values of the random nodes in the domain can be calculated via the MPI method and the ordinary differential equations for the first and second-order moments on the basis of Gaussian closure method. The piecewise linear interpolation based on adaptive least squares is utilized as a post processor to approximate the probability density values on arbitrary positions. The good performance of the resulting method is finally shown in the numerical examples by using three specific non-linear dynamic systems: Duffing oscillator subjected to both harmonic and stochastic excitations, Duffing-Rayleigh oscillator subjected to both harmonic and stochastic excitations, and CHEN system driven by three different Gaussian white noises.