BumiWasesa

Energi bumi yang lestari

BumiWasesa - Energi bumi yang lestari

Geothermal heat diffusion simulation powered by Julia and Python

📘 View Theory & Documentation

Simulation Parameters

Spatial Domain

Time Domain

Default: 1000 years for geothermal simulations
Default: 1×10⁻⁶ (volcanic zone rocks)
Typical: 10-20°C depending on location
Default: 30°C/km (geothermal area)
Default: 2700 (granite)
Default: 790 (granite)
📘 Heat Capacity Conversion:

Volumetric heat generation (q) in W/m³ is automatically converted to temperature change rate using the material's heat capacity:
q[K/s] = S[W/m³] / (ρ × cp)
This ensures physically realistic temperature increases. For granite (ρ=2700, cp=790), 1 W/m³ → 4.69×10⁻⁷ K/s.

Leave empty for 0.01 × duration
Note: Time step (dt) is calculated automatically based on CFL condition with 0.5 safety factor
Formula: dt = 0.5 × (dx² · dz²) / (2κ(dx² + dz²))
Number of time steps (nt) is calculated from: nt = (t_max - t_min) / dt
Geological Time Scale Simulations:
The explicit finite difference method is ideal for geological time scales with realistic thermal diffusivity values (κ ~ 1×10⁻⁶ m²/s for rock).
• For κ = 1×10⁻⁶ m²/s: time step dt ≈ 40 years
• 1000 year simulation: only ~25 time steps needed
• Computationally efficient and numerically stable for diffusion-dominated processes

Heat Sources

Optional: Horizontal Layer Cake Model

Add horizontal layers stacked from top to bottom. Each layer has its own thermal properties. The last layer extends to infinity (halfspace). If no layers are added, the simulation uses uniform background parameters.