【摘要】:边界条件:式中:C0为土壤初始含盐量,本文直接以土壤含盐率百分数表示,%。
模型的初始条件:由于棉花4月中旬播种后到6月中旬没有进行灌溉,所以第一次灌水前的含水率和含盐率为初始值,降雨、蒸发、地下水位均以实测值赋值。典型单元两侧以膜间(未覆膜区)为左右边界。该典型剖面二维模型生育期膜间为大气边界,土面蒸发和降水入渗为实测值;滴头处(2cm直径半圆)为变流量边界1,用实际灌水量表示;滴头以外的覆膜区域为变流量边界2,生育期是为零通量;左右边界为无流量边界;底部为变水头边界,据实测地下水水位赋值。溶质边界上下边界均为第三类边界,左右为无流量边界。选取如图12-15所示模拟区域。
图12-15 模拟剖面示意图
对土壤剖面进行数值模拟计算,即
式中:X、Z分别为模拟计算区域水平和垂直方向的最大距离,cm;(x,z)为计算点的位置坐标;θ0(x,z)为初始含水量,cm3/cm3。
上边界条件:软件模拟的是膜下滴灌棉田,主要分为滴头处、覆膜处和膜间处。膜间处表达式如下:
滴头处表达式如下:
覆膜处表达式如下:
式中:E(t)为蒸发及降雨,cm/d;D(θ)为土壤扩散率,cm2/d;Q(t)为滴头处流量,cm/d;S(t)为蒸腾,cm/d。(www.daowen.com)
下边界条件:模型底部设置为变水头边界,据实测地下水水位赋值。
左边界条件:
右边界条件:
(二)土壤溶质运移初始及边界条件
初始条件可以表示为
式中:C0为土壤初始含盐量,本文直接以土壤含盐率百分数表示,%。
边界条件:
式中:c为已知盐分浓度;q为边界处流量;X、Z为土壤横纵剖面长度。
免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。