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

Ffunction

Calculate right hand side of the differential equation : dx/dt = F(x,t) F is discribed in Hutter 1993.

analyzeResults

Compare analytical and com1DFA results

calcEarlySol

Compute the early solution for 0<t<t_1 to avoid singularity in the Runge-Kutta integration process

computeEarthPressCoeff

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

computeFCoeff

Compute coefficients eq 3.2 for the function F

computeH

get flow thickness from f and g solutions

computeU

get flow velocity in x direction from f and g solutions

computeV

get flow velocity in y direction from f and g solutions

computeXC

get center of mass location

defineEarthPressCoeff

Define earth pressure coefficients

getReleaseThickness

Define release thickness for the similarity solution test

getSimiSolParameters

get flow thickness, flow velocity and center location of flow mass of similarity solution for required time step

mainSimilaritySol

Compute similarity solution :Parameters: simiSolCfg (pathlib path) -- path to simiSol configuration file

odeSolver

Solve the ODE using a Runge-Kutta method

postProcessSimiSol

loop on all DFA simulations and compare then to the anlytic solution

prepareParticlesFieldscom1DFA

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