ana1Tests.simiSolTest
Similarity solution module
This module contains functions that compute the similarity solution for a gliding avalanche on a inclined plane according to similarity solution from : Hutter, K., Siegel, M., Savage, S.B. et al. Two-dimensional spreading of a granular avalanche down an inclined plane Part I. theory. Acta Mechanica 100, 37–68 (1993). https://doi.org/10.1007/BF01176861
Functions
Calculate right hand side of the differential equation : dx/dt = F(x,t) F is discribed in Hutter 1993. |
|
Compare analytical and com1DFA results |
|
Compute the early solution for 0<t<t_1 to avoid singularity in the Runge-Kutta integration process |
|
Computes the earth pressure coefficients function of sng of f and g i.e depending on if we are in the active or passive case |
|
Compute coefficients eq 3.2 for the function F |
|
get flow thickness from f and g solutions |
|
get flow velocity in x direction from f and g solutions |
|
get flow velocity in y direction from f and g solutions |
|
get center of mass location |
|
Define earth pressure coefficients |
|
Define release thickness for the similarity solution test |
|
get flow thickness, flow velocity and center location of flow mass of similarity solution for required time step |
|
Compute similarity solution :Parameters: simiSolCfg (pathlib path) -- path to simiSol configuration file |
|
Solve the ODE using a Runge-Kutta method |
|
loop on all DFA simulations and compare then to the anlytic solution |
|
get fields and particles dictionaries for given time step, for com1DFA domain origin is set to 0,0 for particles - so info on domain is required |