|
Key words: Hanging Wall Deformation, Finite Difference, Crank-Nicholson, Simple Shear
A computer simulation based upon a finite difference model of hanging wall deformation is elaborated via Mathematica Programming. Vertical simple shear is usually adapted to set up the boundary conditions because it facilitates the calculations. Mathematica allows a brief algorithm to generate boundary conditions consistent with inclined simple shear. To find the derivatives of the displacement directions during inclined simple shear, an essential sub-routine in the simulation, two approaches are used. The first involves a numerical solution based upon the Crank-Nicholson method. The second requires the calculus of the derivatives with the support of Mathematica built-in functions.
A comparative analysis is carried out on the results of these sub-routines and on the outcome of the numerical solution of the partial differential equation that governs hanging wall deformation.
The graphical and analytical capabilities of Mathematica are used to generate the trajectories and velocity fields of the hanging wall particles with respect to a Eulerian system of reference.