comsol 多孔介质相对渗透率曲线绘制案例，适合新手入门
打开COMSOL新建模型时，新手容易在茫茫多物理场中迷路。别慌，这次咱们只需要用到两个核心模块：多孔介质流和数学函数库。先选"多孔介质和地下水流"，然后在全局定义里戳开插值函数——这就是绘制曲线的秘密武器。
材料参数设置有个小陷阱：绝对渗透率别手滑输成相对值。记得单位换算，比如把10达西换算成9.87e-12平方米。代码片段里渗透率张量写成对角阵更保险：
model.param.set('k_abs', '9.87e-12[m^2]'); model.component("comp1").material("mat1").propertyGroup("Perm").set("k", {"k_abs","0","0","0","k_abs","0","0","0","k_abs"});这里用全对称设置避免各向异性带来的麻烦，适合新手试水。
饱和度函数定义是重头戏。在全局定义里新建三次插值函数krw(Sw)，把油水两相数据填进去。注意Sw范围强制在0-1之间的小技巧：
if(Sw<0.01, 0, if(Sw>0.99, krw_max, krw_table))这个条件判断能避免计算溢出，之前没加这个的时候我的模型报错报得亲妈都不认识。
后处理阶段别急着导数据，直接在结果表里用参数化扫描。设置Sw从0.2扫到0.8，步长0.05。突然发现曲线出现锯齿？八成是插值方法选成了线性，切到三次样条马上丝滑。最后用全局评估把krw和kro同时拖进表格，右键导出直接扔Excel里画图。
记得保存mph文件时把结果数据打包进去，不然下次打开只有干巴巴的参数。遇到过导出CSV数据错位的惨剧吗？用这个格式字符串保平安：
String.format("%.3f\t%.4e\t%.4e", Sw, krw, kro)三列数据用制表符分隔，小数点位数固定，扔进Origin里都不用再处理。
折腾两小时后突然发现相对渗透率曲线交叉了？别怀疑人生，检查是不是把油相端点饱和度设成了0.2而水相是0.8。这种低级错误我每年都要犯个三四次，正常操作。
