Computational Fluid Dynamics
Thaw
Consolidation
1D Heat Diffusion & Settlement Model — Permafrost Terrain
In Development — Graduate Math Certificate Capstone
This module will simulate thaw front advance through frozen soil using finite difference methods applied to the 1D heat equation, coupled with Terzaghi consolidation to predict ice-rich soil settlement over time. Developed in parallel with a graduate mathematics course. User inputs will include surface temperature, elapsed time, and key soil parameters including P200 fines content as a frost susceptibility proxy.
Planned Input Controls
Parameter 01
Surface Temperature
Mean annual or seasonal temp above freezing (°C)
Parameter 02
Elapsed Time
Simulation duration in days or thaw seasons
Parameter 03
P200 Fines Content
% passing #200 sieve — frost susceptibility index
Parameter 04
Moisture Content
Initial water/ice content — controls latent heat load
Development Phases
01
Walking Skeleton
In Progress
1D finite difference heat equation solver. Animated soil column showing temperature profile and thaw front depth over time. Interactive controls for surface temperature, elapsed time, and soil parameters.
02
Consolidation Coupling
Deferred
Link thaw front depth to Terzaghi 1D consolidation. Display pore pressure dissipation and cumulative settlement alongside thermal output.
03
Full Capstone Model
Deferred
Enthalpy method for phase change, 2D cross-section geometry, unfrozen water content curve for fine-grained soils. Parameterized for Alaska soil types.