Effect of the Water-Air Coupling on the Stability of Rainfall-Induced Landslides Using a Coupled Infiltration and Hydromechanical Model

Algorithm 1

The iteration procedure for solving the displacement-saturation-pressure-porosity system.

1. For the time marching (t_n, t_n+1], let {n_n, , , u_n} be available.

2. Set k =0 and initialize the guess and .

3. Find {, } by solving water and air flow equations (10) and calculate

  ; ; ;

4. Find by solving Eq. (11).

5. Find by solving Eq. (12).

6. Check the convergence criteria:

  IF , THEN collect {, , , } and {, , }, GO TO 7;

  ELSE k = k +1, GO TO 3.

7. Check the time-marching procedure:

  IF , THEN GO TO 8; ELSE GO TO 9.

8. n = n +1, GO TO 1.

9. Output the results.