1D consolidation

1-D Consolidation — Finite Difference Solver

Explicit FD scheme with β = 1/6. Supports open (double drainage) and half-closed (impermeable bottom) layers. Interpolates coarse initial profile to a finer FD grid.

Coarse number of layers (im)

Total layer height h (m)

Coefficient of consolidation c_v (m²/s)

Time to compute t (s)

Refinement factor (default 10 → m = 10 × im)

Layer case Open (double drainage) Half-closed (impermeable bottom)

Beta (β) — stability parameter

Use formulaic n? (shows n from textbook formulas) YesNo (use dt from β)

Initial pore pressure input at coarse depths (one value per coarse node). Enter numbers separated by commas or newlines. Top (z=0) first, bottom (z=h) last. 100 60 0

If number of values doesn't match im+1, the vector will be linearly interpolated between first and last values.

Actions

Computed parameters

Refined intervals m

Nodes (m+1)

dz (m)

dt (s) from β (= β·dz²/cv)

Time steps (n) computed

Textbook formulas (for β ≈ 1/6): open → n = 1.5·m²·Tv (Tv = cv·t/(h/2)²), half-closed → n = 6·m²·Tv. If Use formulaic n = Yes, those n are shown here.

Interpolated initial profile & FD grid

Build to see interpolated grid & initial u(z).

Results / Plot

Not run yet.

Snapshot time (s) — use slider to inspect

Show grid nodes AllCoarse only