水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障碍物复杂边界和层合液体介质声传播问题进行了计算.计算结果表明,与传统时域连续方法如N ew mark方法计算结果对比,所发展方法能较好地消除高梯度和强间断声压力波传播过程中虚假的数值振荡,具有较高的计算精度.问题的求解为进一步流固声耦合问题的研究奠定了基础.     

紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用

利用紧致插值曲线(constrained interpolation profile method in Zhejiang University, CIP-ZJU) 数学模型, 对低科勒冈-卡朋特(Keulegan–Carpenter) 数KC 静止流体中振荡圆柱以及雷诺数Re = 200 时流向强迫振荡圆柱绕流进行了数值模拟. 模型在直角坐标系统下建立, 采用紧致插值曲线方法作为流场的基本求解器离散了纳维-斯托克斯方程, 基于多相流的理论实现流固耦合同步求解, 利用浸入边界方法处理固体边界. 模拟结果与现有文献结果进行比较, 二者吻合情况较好, 验证了此方法对于计算复杂流动问题的可靠性.      

顺序输送管道水力瞬变模拟的有限体积法

针对顺序输送水力瞬变中输送介质不连续的情况,提出采用有限体积法取代传统的特征线法进行求解,以减少插值误差。首先得到水击基本方程基于压力-流速的非守恒有限体积离散格式,通过Rieman求解器得到控制体界面的Rieman解。采用MUSCL-Hancock Primitive方法进行界面数值重构与时间推进,构建时空二阶精度的TVD格式。边界的处理采用Rieman不变量与构造虚拟边界节点相结合的方法,使其与整体精度保持一致。数值计算与实验对比,证明了本文算法能有效抑制非物理振荡,具有精度高、鲁棒性好的特点。     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

粘弹性问题的插值型无单元伽辽金比例边界法

作为一种最近发展起来的半解析数值方法,插值型无单元伽辽金比例边界法不仅无需基本解,且在处理应力奇异性问题和无限域问题时十分有效．为了更有效地求解粘弹性问题,对插值型无单元伽辽金比例边界法应用于此类问题进行了研究,并发展了相应的算法. 通过时域分段展开,将时空耦合的初边值问题转化为一系列递推形式的边值问题,然后采用插值型无单元伽辽金比例边界法进行自适应计算．在径向保持解析特性的基础上,环向采用无单元伽辽金法离散可简化前处理和后处理工作量．此外,改进的插值型移动最小二乘法形函数具有插值性,有效地解决了本质边界条件不能直接施加的困难．最后给出了数值算例,并验证了所提方法的有效性和正确性.     

线性定常对流占优对流扩散问题的无网格解法

应用无网格Galerkin方法求解对流占优对流扩散问题时会出现非物理现象的数值伪振荡,本文将SUPG方法、GLS方法、SGS方法与无网格Galerkin方法相耦合,成功解决了对流扩散方程中对流项占优时的数值伪振荡问题。运用本文构造的方法,采用线性基和具有C2连续的权函数,应用移动最小二乘法可容易地构造高阶导数连续的形函数,从而避免了有限元方法中当采用线性元插值时,因忽略稳定项中二阶导数项而降低计算精度和稳定性的问题。数值实验表明:本文构造的方法具有计算精度高、稳定性好、计算算法实施简单、前后处理方便的优点,这些方法不仅能适用于对流项占优问题,而且也能很好地消除反应项占优时的数值伪振荡问题。     

不可压缩Stokes流动的PSPG无网格法

将应用于有限元法的Pressure-Stabilizing/Petrov-Galerkin(PSPG)稳定化机制引入到无网格法中,有效消除了由于速度和压力的插值模式违反LBB条件而导致的压力场的虚假振荡。采用与有限元法耦合的连续掺混法(Continuous Blending Method)施加本质边界条件,使得边界条件不仅在边界节点上而且在整条边界上都得到严格满足。给出了两个典型算例的数值模拟结果,表明了所建议无网格法模拟不可压缩Stokes流动的有效性。     

三维高频声波的矩阵压缩边界节点法模拟

本文采用边界节点法(Boundary Knot Method, BKM)求解三维高频声场．由于高频赫姆霍兹方程的解是振荡的，极大影响了数值求解的精确度，需要在计算区域增加离散点，这会增加计算量．同时对于大规模声学问题，依靠边界节点法形成的插值矩阵为满秩，导致计算量过高和存储量过大．所以本文采用矩阵压缩技术(Matrix Compression, MC)，在有效继承边界节点法高精确度的基础上减少计算内存需求和时间，从而提高计算效率．数值实验表明，MC-BKM 求解精度高、收敛速度快、计算时间少，在高频大规模声波问题中应用前景广泛．     

基于优化模型的结构声学耦合界面载荷传递算法

针对结构声学耦合系统的界面载荷传递问题,提出了一种基于约束最优化模型的局部参数插值算法。耦合界面处结构和流体单元表面坐标通过各自的形状函数进行插值,把界面载荷在互不匹配的网格节点间的传递问题转化为一个点到用自然坐标表示的有限边界曲面的最小距离问题,以便利用成熟稳定的优化算法对其进行高效求解。与已有方法相比,该算法在耦合界面单元为曲面的情况下仍能保持较高的计算精度。本文给出的数值算例验证了本算法的有效性和可靠性。     

比例边界有限元法求解裂纹面接触问题

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.