基于半隐欧拉-拉格朗日法盐水入侵数学模型

建立了适合河口复杂边界的二维潮流盐度数学模型.其中,网格模块是通过多元最小二乘重构的无结构三角网格;潮流模块基于消除了稳定性条件限制的半隐的欧拉-拉格朗日法,并用干湿判断法实现动边界的处理;盐度输运模块采用有限体积法进行离散,并通过用周围单元平均浓度值重构界面浓度的方法得到与连续性离散方程相协调的二阶对流扩散离散方程.通过纯对流和纯扩散数值测试对模型进行了验证,结果表明模型能够较好地模拟盐度输运的对流扩散问题,且具有较高的精度.最后,将模型应用于长江口盐水入侵的模拟计算,计算结果表明:模型计算的潮位、流速和盐度过程与实测资料一致.