我正在寻找解决理查兹方程。它模拟多孔介质中的水流，是一个非线性的、可能是退化的抛物线微分方程，其形式为

$\partial_t \Theta(\psi) - \nabla \cdot (K(\Theta)\nabla(\psi+z))=0$

其中，是压头，是水含量，是电导率，是垂直高度。压力是主要的未知数。$\psi$$\Theta$$K$$z$

这是一个“著名的”方程，更多信息可以在学者 [1,2] 上找到。

我想计算一个 3D 解决方案，我正在寻找可能的解决方案。

1. Clawpack 适用于这种类型的方程吗？

2. 是否可以将 Clawpack 与 3D 非结构化网格一起使用，例如从 CGAL 生成的网格（目前我有一个使用 CGAL 定义的表面网格）。

非常感谢任何见解。

[1] http://www.sciencedirect.com/science/article/pii/S037704270301001X

[2] http://onlinelibrary.wiley.com/doi/10.1029/WR026i007p01483/abstract