首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 34 毫秒
1.
迎风有限元法在三维潮流数值模拟中的应用   总被引:4,自引:0,他引:4  
董文军  陈虹 《海洋与湖沼》1997,28(3):320-327
为研究和开发迎风有限元法在河口海岸三维波数值模拟中的应用,并为建立三维物质扩散模型,最终为解决实际工程问题打下基础,以非线性浅水波方程为基础,采用一咱平面迎风有奶元与垂向工有限差分相结合地数值计算方法建立三维潮流数学模型,基地实际的物理过程,在计算中将三维流动分成外重力波和内重力波,耦联求解潮位和流速的空间分布。应用本模式计算了天津新港附近渤海海域的三维潮流运动,以较高的分辩率揭示了潮流空间结构特  相似文献   

2.
边界元方法(BEM)是近代计算力学中的一种高效数值方法。由于采用相应的边界积分方程,使处理的问题减少了一维,从而给出的方程组就小得多,要求输入的数据也就大大地减少了,而所得结果的精度却高于有限元法。因此边界元法对于所谓的“区域法”,例如有限差分法(FDM)和有限元法(FEM)具有很强的竞争能力。 大量应用表明:边界元法与有限元法相比,不仅经济、易于使用,而且在许多工程领域,包括某些海洋参数的数值分析和预报方面都是很有发展前途的。 本文首先介绍了边界元方法的基本思想、数学原理和实施步骤,然后分别说明如何将这种方法用于更复杂的、非线性的、依时的问题。最后将讨论边界元法在一般粘性流体流动中的应用。  相似文献   

3.
张二骏  张东生  高飞 《海洋学报》1986,8(5):636-643
本文探讨一种求解工程和数学物理方面具有初边值微分方程的新的数值途径--准分析法。本方法经应用于二维浅水非恒定水流,验算后与实测结果相吻合,并较之按有限差分法和有限单元法所得结果更接近实际。本方法由于使用了各劈开算子的解析解,需要的计算工作量大大减小。又按其基本原理应可以用于求解其他具有初边值的数理方程问题,因而有可能成为一种新的数值求解途径。  相似文献   

4.
为建立高效的Boussinesq类水波数值模型,提出了一种新型的、基于有限差分和有限体积方法的混合数值格式。针对守恒形式的一维控制方程,在等间距矩形控制体内对其进行积分并离散,采用有限体积方法计算界面数值通量,剩余源项采用有限差分方法计算。其中,采用MUSTA格式并结合高精度状态插值方法计算控制体界面数值通量。时间积分则采用具有TVD性质的三阶龙格-库塔多步积分法进行。除验证模型外,重点对MUSTA格式和广泛使用的HLL格式进行了比较。结果表明,MUSTA格式可用于Boussinesq类水波方程数值求解,综合考虑数值精度、计算效率、程序编制和实际应用这几个方面,其较HLL格式更具有优势。  相似文献   

5.
应用比例边界有限元法求解狭缝对双箱水动力的影响   总被引:1,自引:1,他引:0  
比例边界有限元法(SBFEM)是一种半解析数值分析的新方法,既融合了有限元法和边界元法的优点,又有其特有的优点。用该方法可求解有限水深下狭缝对双箱水动力作用的影响,为波浪与多浮体超大型结构的相互作用探索一些规律。整个计算域划分成2个无限子域和4个有限子域,并利用加权余量法在各个子域上推导了SBFEM的积分方程;计算了4个数值算例并与边界元等其它数值方法进行了比较,验证了该方法是一种用很少单元便能得到精确结果的高效方法。应用SBFEM对不同箱体宽度、不同狭缝宽度、不同吃水深度条件的双箱作了计算,得出了狭缝对双箱水动力干涉影响的一些规律,对超大型浮体水动力分析和结构设计具有一定的参考价值。  相似文献   

6.
椭圆型缓坡方程的一个有效的有限元解   总被引:7,自引:0,他引:7  
赵明  滕斌 《海洋学报》2002,24(1):117-123
将海绵层消波的方法用于有限元方法中,提出了缓坡方程的一种有效的有限元求解方法.在应用有限元法求解椭圆型缓坡方程时,通过在方程中加摩阻项,并在入射边界(波浪由此边界进入计算域)处使用不连续单元,将绕射势从总势中分离出来,在边界上利用海绵层进行消波处理,有效地消除了由于引用放射边界条件引起的误差和数值反射现象.  相似文献   

7.
文章以具有自由表面的三维非线性Navier-Stokes方程组为基本控制方程,垂向采用σ坐标变换后得到三维浅海湍流数学模型。数值模拟求解方法采用具有自动迎风、精度高及稳定性好的混合有限分析法求解,最后,以渤海湾潮流模拟算例验证了文中模型及其混合有限分析解法模拟潮流的正确性和应用于工程实际问题的适用性。  相似文献   

8.
海上溢油粒子追踪预测模型中的两种数值方法比较   总被引:2,自引:0,他引:2  
在海上溢油粒子追踪预测模型中,关键的是对拉格朗日微分方程的求解。本文首先通过数值实验比较了欧拉法和龙格-库塔法求解拉格朗日溢油轨迹微分方程的优劣,然后将其应用到2005年4月3日发生在大连附近的“ARTEAGA”油轮溢油事故的油膜粒子追踪模型中。数值实验和应用结果表明,在近岸不均匀流场下,用龙格-库塔方法解拉格朗日油粒子微分方程比用欧拉法求解精度高,用龙格-库塔方法模拟“ARTEAGA”油轮轨迹及其扩散范围与实际观测更为接近,而用欧拉法模拟溢油扩散的面积偏大。  相似文献   

9.
建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。  相似文献   

10.
文章基于Delaunay三角化的思想实现了非结构化数值网格生成技术,并LOP(Local Optimization Procedure)算法对网格进行了优化。用有限体积法在非结构三角形网格上对平面二维潮流方程进行了数值求解。实例应用分析表明,文中的数值网格生成技术和有限体积方法相结合用于求解具有复杂几何边界区域内的潮流问题能取得很好的结果。  相似文献   

11.
众所周知,由于存在插值误差和污染误差的原因,高波数声学问题的有限元解是不可靠的。为了提高有限元法的求解精度,本文提出了一种新型的α有限元法来求解水下声散射问题。在α有限元法中,首先运用基于点的梯度光滑技术得到点光滑的梯度场,然后对点光滑有限元法进行推导。因此,α有限元模型既包含来自点光滑有限元模型的梯度成分又包含来自标准光滑有限元模型的梯度成分,充分利用了点光滑有限元模型的"过软"特性和标准有限元模型的"过刚"特性。为了处理外声场问题,本文运用了Dt N无反射边界条件。数值计算结果表明:与标准有限元法相比,α有限元法在水下声散射计算中具备更高的计算精度和计算效率。  相似文献   

12.
结合抛物型缓坡方程计算波浪辐射应力   总被引:17,自引:6,他引:11  
将波浪辐射应力与抛物型缓坡方程中的待求变量联系起来,提出了一种计算辐射应力的新方法,并用有限差分法对控制方程进行了数值求解。数值结果表明这种方法精度高、编程简单、求解快速,可用于实际大区域波浪辐射应力的计算。  相似文献   

13.
非结构化网格下椭圆型缓坡方程的数值求解   总被引:4,自引:4,他引:0  
魏美芳  唐军  沈永明 《海洋学报》2009,31(2):159-164
椭圆型缓坡方程是一种用线性波浪理论研究近岸波浪传播变形的有效波浪数学模型。非结构化网格下的有限容积法不仅对复杂边界的适应性好,还能保证迭代求解过程的守恒性。建立了非结构化网格下的椭圆型缓坡方程数值模型。在模型中采用非结构化网格下的有限容积法对椭圆型缓坡方程进行了数值离散,结合GPBiCG(m,n)算法求解离散方程。数值计算结果表明,该数值模型可有效地用于模拟近岸缓坡区域复杂边界下波浪的传播。  相似文献   

14.
数值模拟作为海啸预报的主要研究方法在海啸预警中起着关键作用。本文采用Godunov格式的有限体积方法,使用MUSCL-Hancock格式,并利用HLLC Riemann近似求解器计算单元界面上的流体通量,建立了球坐标系下二阶精度的海啸数值模型。模型所基于的全和谐型浅水方程保证了数值的稳定性,而地形重构方法实现了干湿边界的精准模拟。本文模拟了2015年9月16日智利Mw8.3级地震海啸,通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较,验证了模型对实际越洋海啸模拟预报的能力。  相似文献   

15.
曹凤帅  滕斌 《海洋工程》2008,26(1):102-108
比例边界有限元法(SBFEM)是一种半解析的数值方法,比完全数值方法具有更高的精度,该方法结合了有限元和边界元的优点,采用相对少的剖分单元就可以得到较高精度的模拟结果。通过改变有限子域内部比例中心的位置,使这种方法可以应用到多种形式浮体在波浪作用下的水动力特性的计算中。同时还给出了各种形式浮体的波浪力及反射、透射系数的数值结果,并与边界元方法(BEM)计算结果和特征函数展开方法得到解析解进行了比较,均吻合良好。研究表明比例边界有限元不仅可以计算矩形的浮体结构,而且对于多种结构形式的浮体都可以计算,这为多种结构形式浮体的水动力分析提供了一个可行的方法。  相似文献   

16.
为了水下探测和声对抗,舰艇艏部舱室通常安装主动声呐设备,声呐发射声波经由透声窗向外传播,因此,开展透声窗声学特性的相关研究具有重大现实意义。为解决实际工程中透声窗透声性能计算较困难这一问题,一种用于快速预报透声窗的遮挡损失的解析–数值混合方法被提出。该方法用球谐波展开法和叠加原理求解阵列活塞声源声场,得到声激励下透声窗模型表面的声压和法向振速,进而利用COMSOL的Kirchhoff-Helmholtz模块通过再辐射的方法求解模型的遮挡损失。最后,建立了一个二维透声窗模型,分别采用解析–数值混合方法和纯有限元法计算了其全方位的遮挡损失,验证了该方法的准确性、有效性。  相似文献   

17.
初始化方案对有限区域海面风场数值预报模式的影响   总被引:1,自引:1,他引:1  
在一个已应用于海洋要素预报的有限区域海面风场预报模式的基础上,以渤海海域为实验海区,研究初始化方案对有限区域海面风场数值预报模式的影响。对“未采用初始化”与“采用初始化”方法进行了预报实验比较,结果表明,初始化方法对有限区域海面风场模式作用明显,尤其是海面风场的开始阶段,主要体现在海风与陆风预报能力的差异上。  相似文献   

18.
论文研究了多体有限段模型中缆索弹性的两种处理方法:段间弹簧模拟弹性和采用弹性缆段。处理弹性问题,弹性缆段是一个更好的选择,它将多体运动力学和弹性力学结合起来,建立了由弹性缆段组成的多体有限段模型,理论上,这种方法对弹性的处理效果类似采用杆单元考虑几何非线性的有限元法,用于求解运动响应。水下缆索的附加质量作用性质与缆索自身惯性相同,文章将其归入了广义惯性力,并将其它流场力等效为多体方法需要的作用于质心的力系。  相似文献   

19.
长江口水域三维水动力数值模拟研究   总被引:1,自引:0,他引:1  
本文在前人对河口水域三维水动力研究的基础上,建立了长江口水域三维非线性斜压浅海与陆架模式,并首次应用到长江口水域尾水排放水动力数值模拟研究中。考虑到实际河口形状和海底地形,采用平面和垂向三维坐标变换,在空间完全交错的网格系统上离散化控制方程组,并引入ADI计算方法数值求解。数值模拟结果基本上反演了长江口水域的三维水动力状况。  相似文献   

20.
白玉川  于天一 《海洋学报》1998,20(5):126-135
以Navier-Stokes方程为基础,经代数坐标变换导出模式方程,然后在该坐标系下将三维流场进行垂直分层,并在每一层内积分模式方程。在数值模拟求解方面,采用分步、特征差分近似与Galerkin集中质量有限元相结合的方法进行求解。在实际应用方面,对广西廉州湾潮流进行了三维模拟计算。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号