共查询到17条相似文献,搜索用时 93 毫秒
1.
2.
根据研究工程、工业废水及生活污水的水环境影响的需要,提出一个平面二维水流-水质有限体积法及黎曼近似解模型。在无结构网格上对偏微分方程进行有限体积的积分离散,把二维问题转化为一系列局部的一维问题进行求解,模型具备有限单元法及有限差分法两者的优点。由于模型采用黎曼近似解计算水量、动量及污染物输运等通量,不仅提高了数值模拟的精度,而且能模拟包括恒定、非恒定或急流、缓流等水流-水质状态。应用若干理想条件下的精确解对模型精度进行了检验,并结合汉江中下游的水质问题验证了模型的实际应用能力。模型计算结果与精确解拟合极好,与汉江中下游实测的污染分布带也相当吻合,为中线南水北调对汉江中下游水质的影响评价提供了依据。 相似文献
3.
4.
采用平面水流、水质有限体积法及黎曼近似解模型,在无结构网格上对偏微分方程进行有限体积法的积分分离散,模型具备有限单元法及有限体积法两者的优点。采用同步遥感、水量、水质实测资料对模型进行了率定验证,在对长江水源地、沿江污染源调查以及岸线稳定段分析结果基础上,进行污染源概化以及确定主要研究江段。在对概化的污染源计算出在设计条件下的污染带分布图,并计算出污染带长、宽随排污量变化曲线;在上述研究基础上,分析出长江各水源地及岸线稳定段的水质变化状况。对长江江苏段各主要排污口的污染带进行了计算,并对长江江苏段区域供水水源地的水质可达性进行了研究。 相似文献
5.
6.
7.
8.
采用一维-二维耦合数学模型对溃堤洪水进行模拟,可以发挥两者优势,提高计算效率。将基于黎曼近似解Godunov格式的一维、二维水流数学模型通过堰流公式进行耦合,实现交界面的水力要素交互;一维、二维模型均采用中心形式的有限体积法求解,分别采用HLL格式和Roe格式计算界面通量,具有一致的空间同步性。通过经典算例验证了模型的可靠性和稳定性,应用耦合模型对1998年松花江干流胖头泡溃堤洪水进行了模拟计算,较好地模拟了溃堤洪水的演进过程与淹没范围,对200年一遇洪水的分洪效果进行了分析,为此区域的分洪利用及松花江中上游防洪调度提供参考依据。 相似文献
9.
基于有限体积方法和结构化网格,建立了海底滑坡引起的波浪传播数值模型。模型控制方程为考虑了海床随时间变化的二维浅水方程。采用中心迎风格式计算控制体界面数值通量,采用线性重构技术、局部海床高程处理技术和全隐式离散底摩阻项,保证了格式的和谐性、守恒性和水深非负性,有效处理了海岸动边界问题。时间积分采用具有强稳定性质的二阶龙格-库塔方法(ssp RK)。针对滑坡体兴波经典算例开展数值模拟,将计算结果与解析解、实验结果及其他模型计算结果进行比较和分析。结果表明,对于所考虑的计算工况,模型能较合理地模拟滑坡兴波的产生、传播和爬高过程。 相似文献
10.
11.
12.
为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。 相似文献
13.
针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。 相似文献
14.
天然水体具有形状复杂的计算区域和水下地形。本文采用无结构的网格以适合这一情况,并能方便地根据工程应用的要求局部地和适应性地加密网格。相应地,对二维水流计算问题采用了有限体积法的数学表示。跨单元边界的法向数值通量通过用Osher格式求解黎曼问题得到。文中给出了适合于二维浅水方程组的有关单元界面和各种物理边界的法向数值通量公式。分析了这一格式所具有的优点,包括:普适性、守恒性、逆风性、单调保持性、高效性、对间断的高分辨率、边界处理和内部格式相容,以及不引入数值边界条件等。最后,通过长江口南支水流计算的实例,阐明其良好性能。 相似文献
15.
边界层参数化影响“梅花”台风的敏感性试验 总被引:3,自引:0,他引:3
以GRAPES-TCM为试验模式,对1109台风“梅花”进行了36次72 h的预报试验,通过试验分析了2种边界层参数化方案——MRF方案与YSU方案在不同情况下对台风预报的影响.结果显示:“梅花”路径与强度对边界层方案的变化都表现出一定的敏感性,敏感性大小与对流参数化方案、台风的初始强度等因素有关,强度的敏感性比路径更明显;对弱台风的路径与强度,YSU方案的总体预报效果优于MRF方案,对于强台风,2种边界层方案中MRF方案的路径预报效果更好,哪种方案的强度预报效果更好与对流参数化方案有关;无论何种情况,YSU方案预报的“梅花”强度都明显强于MRF方案,YSU方案预报的降水及感热通量与潜热通量总体上大于MRF方案;YSU方案时更多的感热通量和潜热通量与该方案时边界层更强的湍流混合有关,更多的潜热通量导致更多的降水,从而释放更多的潜热,更多的潜热释放以及更多的感热通量导致台风强度更强. 相似文献
16.
青藏公路沿线多年冻土区土壤热通量参数化方案的优化和检验 总被引:1,自引:1,他引:0
土壤热通量是地表能量平衡的重要分量,其估算方案在研究地表能量平衡研究中必不可少。利用青藏公路沿线5个站点0~20 cm的实测土壤层温、湿度及5 cm土壤热通量资料,以翁笃鸣气候学计算方案为基础建立了优化的5 cm土壤热通量计算方案。通过唐古拉和西大滩两个独立站点的检验结果表明,优化方案的结果相对于原方案有较大的改善,唐古拉和西大滩5 cm土壤热通量均方根误差值分别减小了3.2 W·m-2和4.8 W·m-2,而相对误差分别减小了61.9%和36.1%,即新方案能够较好地估算出青藏公路沿线多年冻土区5 cm土壤热通量。使用优化方案模拟了青藏公路沿线11个站点5 cm土壤热通量变化,结果显示,近十年青藏公路沿线土壤热通量呈现出增大的趋势,其中,5 cm土壤热通量增大了近1.0 W·m-2,而且各观测场的年平均土壤热通量值均大于0.0 W·m-2,表明就年尺度而言,热量有盈余,盈余热量用于加热下层土壤,引起活动层厚度增加,平均状况下土壤热通量每增大1.0 W·m-2,活动层厚度增大约21.0 cm。 相似文献
17.
在交错网格上,基于高精度的Gauss积分公式,针对浅水波方程设计了对模拟涌波具有高分辨率的完全二阶精度的数值计算格式。由于采用了交错网格,差分格式不需要解Riemann问题,因此本文格式具有计算简单、工作量少、编程简便等特点。另外,在一维单个方程时,本文格式在CFL(Courant Friedrich Lewy)条件限制下为TVD(Total Variation Diminishing)格式,在二维和三维情况下格式具有MmB(Maximum and Minimum Bounds Preserving)性质。利用国家高性能计算中心(合肥)的曙光1000型分布存储大规模并行机,对在交错网格下所构造的求解浅水方程的高分辨差分格式进行了并行实现,几个算例的计算结果令人满意。 相似文献