PhysSandbox
Classical MechanicsWaves & SoundElectricity & MagnetismOptics & LightGravity & OrbitsLabs
🌙Astronomy & The Sky🌡️Thermodynamics🌍Biophysics, Fluids & Geoscience📐Math Visualization🔧Engineering🧪Chemistry

Related simulators

Continue with similar topics in this category — or all 44 in Biophysics, Fluids & Geoscience.

View category →
NewSchool

Darcy Flow & Penetration Depth

Launch Simulator

1D head h(z,t): ∂h/∂t = K ∂²h/∂z²; recharge at surface, fixed head at depth — relaxation time and diffusion length.

NewUniversity / research

Infinite Slope Stability

Launch Simulator

Limit-equilibrium factor of safety for a shallow planar slide: slope angle, cohesion, friction, depth, and rainfall-driven pore pressure set the landslide threshold.

NewUniversity / research

Ekman Spiral in the Ocean

Launch Simulator

Wind stress, Coriolis f, eddy viscosity A_z: classical Ekman spiral of horizontal velocity with depth; deflection angle.

NewUniversity / research

Kelvin–Helmholtz Shear Instability

Launch Simulator

Two parallel flows with slip U: Fourier-mode growth rate sketch and animated interface billows (linear toy model).

NewSchool

0-D Energy Balance & Ice–Albedo

Launch Simulator

Budyko-style (1−α)S/4 = εσT⁴ with icy vs clear albedo, freeze/melt hysteresis when S is swept, and optional smooth α(T) relaxation.

NewSchool

Layered Medium: P/S Fronts & Hodograph

Launch Simulator

P and S in a two-layer straight-ray toy: reflected legs and optional head branches when the lower layer is faster; hodograph plus minimum-time first arrivals (no conversions or amplitudes).

PhysSandbox

Interactive physics, chemistry, and engineering simulators for students, teachers, and curious minds.

Physics

  • Classical Mechanics
  • Waves & Sound
  • Electricity & Magnetism

Science

  • Optics & Light
  • Gravity & Orbits
  • Astronomy & The Sky

More

  • Thermodynamics
  • Biophysics, Fluids & Geoscience
  • Math Visualization
  • Engineering
  • Chemistry

© 2026 PhysSandbox. Free interactive science simulators.

PrivacyTermsContact
Home/Biophysics, Fluids & Geoscience/Glacier Flow 1D (Glen’s Law)

Glacier Flow 1D (Glen’s Law)

Shallow-ice thickness H(x,t) on a tilted bed: nonlinear flux q ∝ H^{n+2}|∂s/∂x|^{n−1}∂s/∂x plus accumulation; τ_b ≈ ρgH sinα readout.

SIA + Glen (toy)

8°
3
0.035
0.012

1D shallow-ice approximation: ∂H/∂t = −∂q/∂x + M with q ∝ H^{n+2}|∂s/∂x|^{n−1}∂s/∂x, s = b+H, x∈[0,1]. Left: zero flux (divide). Right: H=0 (calving). τ_b = ρgH sinα uses thickness at the dashed probe ×320 m for a kPa scale — schematic, not an inversion.

Measured values

Model time (arb.)0.00
Max thickness (model)0.000
Mean H (model)0.000
τ_b readout (probe)0.00 kPa

Real glaciers need 2D/3D thermomechanics and sliding; this page isolates how a nonlinear diffusive flux redistributes thickness once surface slopes drive q in Glen’s n-law form.

Live graphs

About this model

The shallow-ice approximation (SIA) treats ice as slowly sliding material where vertically integrated flux follows Glen’s (power) law: effective diffusivity grows strongly with thickness and surface slope. In 1D, conservation ∂H/∂t + ∂q/∂x = M with q ∝ −H^{n+2}|∂s/∂x|^{n−1}∂s/∂x and surface s = b + H produces characteristic thickening near divides and drawdown toward a calving front. This page uses explicit finite volumes on x ∈ [0,1] with a planar bed tilt α, zero flux at the divide, H = 0 at the terminus, and uniform accumulation M for intuition only — no thermomechanical coupling, sliding, or hydrostatic bed.

Who it's for: Cryosphere introductions after ODEs; companion reading to SIA derivations in glaciometry texts.

Key terms

  • Shallow-ice approximation
  • Glen’s law
  • Ice flux
  • Basal shear stress
  • Accumulation

How it works

Cross-section evolution of ice thickness on a planar bed: nonlinear diffusion from Glen’s flow law in the shallow-ice approximation, with a textbook basal shear τ_b ≈ ρgH sinα readout at a fixed probe.

Key equations

∂H/∂t + ∂q/∂x = M; q ∝ −H^{n+2}|∂s/∂x|^{n−1}∂s/∂x; s = b + H; τ_b ≈ ρ g H sin α

Frequently asked questions

Are ρ, g, and the τ_b number field-calibrated?
No. ρ and g are standard constants for the τ_b ≈ ρgH sinα panel; model thickness is rescaled by a fixed factor for display in kPa. The evolution uses a lumped mobility control, not a measured rate factor A(T).
Why can the profile look noisy at high mobility?
Explicit integration of a strongly nonlinear diffusion needs small effective Courant numbers; raise mobility slowly or pause if short-wavelength ripples appear.