GA入渗模型的湿润峰数值解与解析解 [1]模型简介：使用Python 和数值模拟软件COMSOL，复现EI（潘永亮,简文星,李林均,等.基于改进Green-Ampt模型的花岗岩残积土边坡降雨入渗规律研究[J].岩土力学,2020,41(08):2685-2692.），以Green-Ampt入渗模型对无限边坡进行降雨入渗的数值模拟，并同时使用Python计算湿润峰深度的解析解结果以此形成数值解与解析解的对比分析 [2]案例内容：COMSOL数值模拟模型（包括边界条件的设置、后处理结果等），Python解析解代码以及讲解稿一份 [3]模型特色：（1）掌握GA入渗模型在COMSOL数值模拟软件中的设置；（2）了解数值模拟过程中不同降雨边界条件设置的差异；（3）掌握数值解与解析解二者的实现方法 COMSOL6.2版本

搞边坡降雨入渗模拟的兄弟肯定绕不开Green-Ampt模型。这玩意儿虽然是个简化模型，但在工程上够用且不折腾。今儿咱们拿COMSOL6.2和Python配合着玩，看看数值解和解析解到底能差出几个汉堡钱。

先说核心公式。GA模型认为湿润峰像把尖刀垂直插进土里，湿润区与非湿润区泾渭分明。湿润峰深度Z(t)满足这个微分方程：

dZ/dt = (K_s(ψΔθ + Z)) / (ZΔθ)

GA入渗模型的湿润峰数值解与解析解 [1]模型简介：使用Python 和数值模拟软件COMSOL，复现EI（潘永亮,简文星,李林均,等.基于改进Green-Ampt模型的花岗岩残积土边坡降雨入渗规律研究[J].岩土力学,2020,41(08):2685-2692.），以Green-Ampt入渗模型对无限边坡进行降雨入渗的数值模拟，并同时使用Python计算湿润峰深度的解析解结果以此形成数值解与解析解的对比分析 [2]案例内容：COMSOL数值模拟模型（包括边界条件的设置、后处理结果等），Python解析解代码以及讲解稿一份 [3]模型特色：（1）掌握GA入渗模型在COMSOL数值模拟软件中的设置；（2）了解数值模拟过程中不同降雨边界条件设置的差异；（3）掌握数值解与解析解二者的实现方法 COMSOL6.2版本

别急着掏笔算，Python解这个有更骚的操作。直接上代码：

import numpy as np
import matplotlib.pyplot as plt

Ks = 5e-6  # 饱和导水率(m/s)
psi = 0.8  # 湿润锋基质吸力(m)
theta_s = 0.45  # 饱和含水率
theta_i = 0.15  # 初始含水率
q_rain = 1.2e-6  # 降雨强度(m/s)

delta_theta = theta_s - theta_i
t = np.linspace(0, 6*3600, 100)  # 模拟6小时

# 解析解计算公式
Z = (Ks*psi*delta_theta)/(q_rain - Ks) * (np.exp((q_rain - Ks)*t/(psi*delta_theta)) - 1)

这段代码的精髓在最后那个指数项。当降雨强度q_rain超过Ks时，分母会变负，这时候模型其实已经不适用了——说明地表要开始积水了。这个临界点在实际工程中得特别注意。

转到COMSOL这边，设置比代码复杂但更有画面感。关键几步：

1. 在物理场选"多孔介质和地下水流"，别手滑选成达西流
2. 初始条件设置孔隙水压力为-psi（吸力为负压）
3. 上边界用通量边界，下边界设自由排水
4. 材料参数别直接填数值，用变量表达式方便后期调试

有个坑爹细节：COMSOL默认压力单位是Pa，但咱们的参数是米水头。记得在变量定义里加个转换系数rho*g，别问我怎么知道的——说多了都是凌晨三点的眼泪。

跑完模拟把结果拉出来遛遛：数值解的湿润峰前期推进比解析解慢，2小时后开始反超。这可不是误差，解析解假设的是无限边坡瞬态入渗，而数值模型考虑了实际土柱的边界效应。用Python画个对比图立马现形：

plt.plot(t/3600, Z, label='解析解')
plt.plot(comsol_time, comsol_z, '--', label='COMSOL数值解')
plt.xlabel('时间(h)')
plt.ylabel('湿润峰深度(m)') 
plt.legend()
plt.show()

两条曲线像极了异地恋的情侣——开始步调不一致，处久了反而趋向同步。这说明在长时间模拟时两种方法可以相互验证，但短时行为差异能到15%，选方法得看具体需求。

最后说点干活心得：解析解适合快速估算和参数敏感性分析，改个Ks值秒出结果；数值解虽然设置麻烦，但能玩各种骚操作——比如模拟间歇性降雨时，改个边界条件类型分分钟切换干湿状态。下次要是老板催着要结果，先甩解析解应付，背后偷偷用COMSOL调参才是打工人正确姿势。