论文部分内容阅读
本文利用数值方法研究了海上建筑物周围的泥沙冲刷和海床演变,建立了用于模拟建筑物周围海床演变的数值模型。由于直接利用三维模型计算量非常大,为了提高计算速度和节省计算机存储量,本文利用深度平均的二维模型计算波浪场和水流场,并利用在海底边界层内对动量方程进行垂向积分的数值模型计算海底剪切应力。本文中所有的数值方法利用有限元方法实现。 本文提出的数值模型主要用于模拟大尺度建筑物周围的泥沙冲刷情况,所以在计算时必须考虑波浪对建筑物的绕射和反射。在本文中首先分别利用有限元方对缓坡方程和Boussinesq方程进行了数值求解,并在此基础上分别建立了基于缓坡方程和Boussinesq方程的波浪作用下建筑物周围的海床演变数值计算模型。在求解波浪方程时应用了边界处利用解析解拟合和内部边界造波、边界阻尼层消波的方法,避免了利用数值方法近似处理边界条件引起的误差。在求解波浪引起的水流场时利用二维浅水方程,将辐射应力梯度项加入浅水方程中以考虑波浪对水流的影响。计算海底剪切应力时应用波流共同作用模型,而不是纯波浪和纯水流两种情况下的剪切应力的叠加。计算了波浪周期内各个时刻瞬时的悬移质和推移质泥沙输沙率,并对输沙率进行时间平均求得时间平均的输沙率。计算输沙率时考虑了Lagrangian漂移速度和由于波浪边界层不均匀产生的竖向流速的影响,在计算海底地形演变时通过求解用输沙率表示的泥沙质量守恒方程来实现。推导出了利用波面随时间的分布来表示的辐射应力计算公式。对泥沙在斜坡上运动时的临界Shields数进行了修正,并将其应用于本文提出的数值模型中。利用本文数值模型的计算结果与他人的数值结果和试验结果进行了比较。利用数值方法研究了KC数、泥沙粒径和模型比尺对圆柱周围的泥沙冲刷情况的影响。 本文利用物理模型试验对波流共同作用下大尺度圆柱周围的泥沙冲刷进行了研究,研究了在波流作用下不同直径的圆柱周围的泥沙冲刷情况,并建立了用于模拟波流共同作用下建筑物周围的泥沙冲刷的数值模型,将此数值模型的计算结果与物理模型试验中的实测值进行了比较。由于实际工程中水流速度相对波浪水质点速度很小,所以在数值模型中忽略了水流对波浪的影响,而考虑波浪辐射应力对水流的影响。