无限水深内完全非线性数值波浪水槽建立

基于无旋、不可压的势流理论,利用高阶边界元法建立了一种可应用于无限水深的无粘三维完全非线性数值波浪水槽.研究中,自由水面条件为完全非线性边界条件,进而可以模拟任意强的非线性波浪.将镜像格林函数应用到无限水深的数值波浪水槽中,以至于两水槽侧壁可以排除于计算域外.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置在计算域内,同时在出流边界布置人工阻尼层来吸收出流波浪.线性的和完全非线性数值实验验证了所建立的数值模型可以用来产生无限水深的任何指定波浪,计算结果与解析解吻合得很好,且在出流边界没有反射.     

数值波浪水槽造波技术及工程应用研究

本文在充分了解数值水槽模型、数值造波、消波理论基础上,不仅分别利用设置造波边界法和质量源项法在Fluent软件中模拟出较为有效的线性波和弱非线性的二阶Stokes波,还将建筑物添加到所构造的波浪数值水槽中,利用波浪数值水槽不拘泥于特定波浪的特点,按照有已知实测资料的波浪和建筑物相互作用的波浪水槽进行建模,提出了适合不同建筑物和不同造波方法的网格构建方法,优化了水槽的构建,并验证了水槽的有效性,最后分析了两种方法的优缺点及适用性。本文研究成果为波浪与建筑物相互作用的深入研究奠定了基础。     

二维数值造波水槽不规则波之数值研究

利用边界元素法所建立数值造波水槽探讨不规则波浪之造波与传递问题，以Lagrange法描述水粒子之运动并利用前进差分处理时间项之微分项。水槽一端为活塞推拉式造波板，运动方式以Mitsuyasu-Bretschneider之波谱造波，另一侧设有一段消波区以减少计算结果受波浪反射之影响。水面上选定若干个测点量取水位变化之历时，利用波谱分析水面上各测点之波谱并与原入射波波谱比较来探讨本模式之准确性。     

二维数值波浪水槽在FLUENT中的实现

为在波浪水槽内得到稳定的、重复性良好的波动过程,基于N-S方程和标准k-ε湍流模型,采用VOF方法追踪自由面,利用CFD商业软件FLUENT建立同时具有造波和消波功能的二维数值波浪水槽.基于线性波浪理论和UDF实现动边界造波和定义水质点速度造波,同时,在数值波浪水槽尾端通过在动量方程中添加附加源项的方法实现动量源消波.计算结果表明,基于FLUENT软件实现的二维数值波浪水槽工作特性良好.     

基于OpenFOAM的数值波浪水槽研究

使用设置造波边界法和质量源造波法建立了数值波浪水槽,采用阻尼消波法进行消波.设置造波边界法根据行波的解析解或数值解来给定边界上的流速和波高值.质量源造波法是在水槽中特定位置,将质量源项加入到质量守恒方程中.使用该数值水槽模拟了线性波浪和二阶Stokes波,并与理论解进行了对比,结果显示两种造波方法都能有效模拟行波.设置造波边界法必须指定速度和波面边界,不适用于波流相互作用等边界未知的情况,而源造波法则能适用.但是由于源造波法需要在两端设定造波区域,计算区域更大,其造波效率低于设置造波边界法.     

基于边界造波法的二阶Stokes波的数值生成

针对滑坡涌浪灾害中波浪的产生和传播分析至关重要的问题,利用FLUENT软件平台,提出了边界造波的数值方法,并通过二次开发实现了二阶Stokes波的数值模拟。该边界造波法以不可压缩粘性流体的Navier-Stokes方程和追踪自由面的VOF法为基础,借助FLUENT软件中的UDF功能定义了入射边界处波浪的速度和波形解析表达式,从而模拟了二阶Stokes波的产生和传播。计算结果表明,该数值水槽能有效可靠的模拟弱非线性波浪。这为后续模拟由滑坡涌浪所产生的波浪提供了一定的研究基础。     

基于改进质量源造波方法的非线性波数值模拟

传统的质量源造波方法在强非线性波浪模拟过程中精度不高,不同波况造波源区初始设置较为复杂,针对这些问题本文提出了改进方法.考虑质量源区双向质量输出(横向和垂向),改进质量源源强函数并进行理论推导;基于OpenFOAM开源程序开展网格和时间步长收敛性验证;讨论在不同相对水深以及不同非线性强度条件下,二阶Stokes波幅值和相位与文献中传统方法模拟结果的相对误差,并对孤立波的传播过程进行模拟分析.数值模拟结果表明,改进方法不仅能提高效率,而且能提高质量源造波在较强非线性范围的模拟精度,尤其在非线性强度参数为0.2时,波幅及相位误差相对于传统方法分别降低21.7%和46.5%.     

基于边界造波法的波浪数值生成

一般边界造波法只是用来生成线性波,至于边界造波法对不规则波、孤立波和强非线性波等是否适用,并没有做过多分析与研究。基于线性波、二阶Stokes波以及孤立波的速度和波形解析表达式,采用边界造波法,借助FLUENT软件中的用户自定义(UDF)功能,通过二次开发在FLUENT中实现波浪的模拟,并对它们的数值解与解析解做出比较,以此来验证边界造波法的适用性。计算结果表明,该数值水槽能有效模拟这三种波浪。     

基于二阶规则波模拟的数值波浪水槽优化

基于CFD开源软件包OpenFOAM,建立了基于Navier-Stokes方程的二维数值波浪水槽模型。模型采用体积函数(VOF)法捕捉自由水面,采用 k-ω SST紊流模型封闭控制方程。针对不同波陡( i =0.01~0.06)的Stokes二阶波,分别采用速度边界法或推板造波法进行造波,采用阻尼消波法、辐射边界法或松弛区域法进行消波;通过对不同造波和消波方法组合方式计算结果的比较,分析了数值波浪水槽在不同波陡工况下的消波效果和造波效果。模拟结果表明:对于小波陡( i =0.01~0.03)工况,采用速度边界造波和松弛区域法消波的组合方式具有较高的模拟精度;对于大波陡( i =0.04~0.06)工况,采用推板造波法造波和基于平方形式阻尼函数的阻尼消波法消波的组合方式具有较高的模拟精度。     

基于光滑粒子方法的推板式波浪水槽实验研究

波浪数值水槽是波浪研究方面的一个十分重要的实验装置。目前针对造波水槽所展开的数值或实验研究,多将造波系统的设定参数作为波浪的参数,而未考虑设定参数与实际参数之间的误差。本文利用推板造波的实验装置进行造波,并借助视频拍摄来捕捉造波板的运动,而波浪参数则由波高仪的数据采集来得到。推板造波的过程,则用光滑粒子方法来进行数值模拟。造波过程中,造波板的运动同样采用实测的参数,并计算了不同时刻的波形曲线和水面高度随时间的变化。光滑粒子的模拟结果与实验结果比,吻合的非常好。     

退化子单元方法在构造数值波浪池中的应用

采用高阶边界元方法构造了一个三维数值波浪水池。为了解决奇异积分问题 ,在单元划分中联合采用连续元和不连续元同时引入一种在退化子单元中进行单元细分的处理方法。在数值波浪水池中模拟的线性波和非线性 Stokes波 ,都得到了较好的结果 ,表明所采用的高阶边界元方法以及开边界条件有着很好的计算精确性和稳定性     

用ALE有限元法建立二维数值波浪水池

从ALE描述下的Navier-Stokes方程和连续方程出发,对数值波浪水池进行分析、通过对各种边界条件的处理建立了合适的二维数值波浪水池。通过对二阶Stokes波的探讨,阐明了此数值波浪水池是有效合理的。     

波浪对船剖面作用的时域计算模型

为了求港口中波浪对船的作用力,用边界元数值求解复速度势建立了波浪对船剖面作用的时域计算模型,采用满足运动学和动力学条件的半拉格朗日法跟踪流体的自由表面。推导了适用本数值模型的辐射边界条件,此边界条件和海绵层边界条件组成的出口边界条件,可以使波浪很好地被吸收或传播出去。数值计算结果和实验结果对比表明此模型能精确模拟波浪对船的作用力。     

有限元与边界元耦合法在波浪力计算中的应用

利用有限元与边界元耦合法对三维无界区域中直立圆柱所受的波浪力进行计算,把整个求解区域分成内域和外域两部分,在内域采用有限元法,对外域采用边界元法计算,由加权余量法的理论知这种耦合在理论上是可行的,根据此耦合方法编制了完整的计算机程序,对海洋中直立圆柱的波浪力进行了分析.数值计算的结果与理论解吻合良好,表明该方法有效.     

有限元波动模拟中傍轴近似透射边界

有限元波动模拟中傍轴近似透射边界张剑锋（大连理工大学工程力学系１１６０２４）关键词：地震波；有限元法／傍轴近似；人工边界分类号：Ｏ２４２．２１；Ｐ３１５．３实际地震波是在半无限域中传播，但有限元有限差分等数值方法只能在有限尺寸模型上求解；因此数值模拟...     

基于图像八叉树的三维比例边界有限元多面体网格生成算法

基于图像八叉树方法,提出了平衡八叉树和多面体网格修剪相结合的三维比例边界有限元多面体网格算法,该算法根据结构尺寸建立恰好完全包含整个结构的立方体网格,再建立结构图像像素信息,根据像素信息,按照2∶1的平衡分割原则递归地进行等分,完成平衡八叉树网格生成。结构内单元网格完全保留,结构外单元网格删除,对于结构边界单元网格,提出采用面-立方体相交判断方法进行边界单元与结构边界相交面的筛选,搜寻结构单元与结构边界表面相交点,通过有序连接相交点形成边界单元切割面,再结合边界单元其他几个面,构成裁剪后的多面体单元。数值算例结果表明,基于本文算法生成的比例边界有限元网格计算结果具有较好的精度和边界适应性。     

数值网格生成方程中源项对网络正交性影响的研究

基于对数值网格生成的椭圆型偏微分控制方程中混合偏导数项非线性系数特性的 研究,从理论上推导了用于控制网格正交性的一种新源项给定形式,并讨论了其具体的数值实 现方法,数值计算结果表明,该新源项形式对生成网格的正交性具有很强的调节和控制能力。