POLLUTEv10 Example 16: Monte Carlo Simulation of Leachate Collection System Failure Timing

Monte Carlo simulation diagram showing variability in landfill leachate system failure timing
Share the knowledge

Introduction

POLLUTEv10 Example 16 introduces probabilistic modeling into landfill contaminant transport analysis using Monte Carlo simulation. Building on Example 15, this case evaluates how uncertainty in the timing of primary leachate collection system (LCS) failure affects contaminant migration.

Rather than assuming a single failure time, this example models a range of possible failure scenarios, providing a more realistic assessment of long-term environmental risk.


Conceptual Model Overview

This model is identical to Example 15 except for one key enhancement:

👉 The onset of primary LCS failure is uncertain

Key Features

  • Finite contaminant source
  • Passive Sink for secondary leachate collection
  • Variable Properties for time-dependent flow
  • Monte Carlo simulation for uncertainty analysis

Modeling Uncertainty in System Failure

Instead of a fixed failure time, the model uses a triangular probability distribution:

ParameterValue
Minimum failure time20 years
Most likely (mode)25 years
Maximum failure time50 years

Interpretation

  • Some simulations show early failure (20 years)
  • Most cluster around 25 years
  • Others extend to late failure (up to 50 years)

This approach captures the natural variability and uncertainty in landfill system performance.


Monte Carlo Simulation Approach

Monte Carlo simulation works by:

  1. Randomly sampling a failure time from the triangular distribution
  2. Running the POLLUTE model for that scenario
  3. Repeating the process many times
  4. Analyzing the distribution of results

Output

Instead of a single deterministic result, you obtain:

  • Range of contaminant concentrations
  • Probability distributions
  • Risk-based insights

Hydraulic Behavior with Variable Failure Timing

As in Example 15:

  • Initial Darcy velocity: 0.01 m/a
  • After failure: 0.1 m/a

However, the timing of this transition varies per simulation.

Implications

  • Earlier failure → greater contaminant migration
  • Later failure → reduced long-term impact
  • Results reflect uncertainty envelope, not a single prediction

Passive Sink and Variable Properties Integration

The modeling framework remains consistent:

Passive Sink

  • Represents secondary leachate collection system
  • Provides lateral drainage

Variable Properties

  • Controls time-dependent Darcy velocity
  • Now linked to random failure timing

Important Note

  • Darcy velocities from both features are multiplied internally
  • Best practice:
    • Set one feature to 1.0
    • Input actual values in the other

Model Parameters

All parameters are identical to Example 15, with the addition of Monte Carlo inputs:

PropertyValueUnits
Diffusion Coefficient0.02m²/a
Dispersivity0.4m
Porosity (soil)0.4
Porosity (granular)0.3
Source Concentration1000mg/L
Reference Height7.5m
Landfill Length200m
Aquifer Velocity4m/a

Monte Carlo Parameters

ParameterValue
Minimum failure time20 years
Mode25 years
Maximum50 years

Graphical Output: Probability vs Concentration

PDF Report

Loader Loading…
EAD Logo Taking too long?

Reload Reload document
| Open Open in new tab

Key Insights

  • Uncertainty significantly affects predicted contaminant migration
  • Deterministic models may underestimate or overestimate risk
  • Monte Carlo simulation provides:
    • Probability-based outcomes
    • Better support for decision-making
  • Early failure scenarios dominate worst-case risk

Practical Applications

This type of analysis is valuable for:

  • Risk assessments
  • Landfill design optimization
  • Regulatory submissions
  • Long-term monitoring strategies

It helps answer questions like:

👉 “What is the probability that contamination reaches the aquifer within 50 years?”


Numerical and Modeling Considerations

  • Results depend on number of Monte Carlo realizations
  • More simulations → smoother probability distributions
  • Computational demand increases with:
    • Number of realizations
    • Model complexity

Important Disclaimer

⚠️ This example is for demonstration purposes only.

  • Not a design guideline
  • Not suitable for direct application without expert review

Use of Monte Carlo simulation with Variable Properties requires:

  • Advanced understanding of POLLUTEv10
  • Strong background in hydrogeology and risk modeling
  • Consultation with program developers for critical applications

Conclusion

POLLUTEv10 Example 16 elevates contaminant transport modeling by incorporating uncertainty into system performance. Through Monte Carlo simulation, it provides a more realistic and robust framework for evaluating landfill risks over time.

This approach reflects real-world conditions—where system failure is not a fixed event, but a probabilistic process—making it a powerful tool for modern environmental engineering.


Learn more about our Contaminant Transport Modeling Solutions

POLLUTE Examples

Comparison between POLLUTE and MIGRATE

1 / ?