加权最小二乘无网格法的随机稳态温度场分析

对加权最小二乘无网格法在随机稳态温度场中的应用进行了研究.在移动最小二乘近似的基础上,采用罚函数法满足边界条件,通过变分原理详细推导了求解稳态温度场问题的加权最小二乘无网格公式,与无网格伽辽金法相比,该方法无须进行高斯积分,具有计算量小、处理方便等优点.同时考虑结构物理参数和边界条件随机性的影响,利用Neumann展开蒙特卡罗法对含有随机参数温度场的加权最小二乘无网格方程进行求解,得到了温度场响应量的统计特征值并考察了各随机参数对节点温度的影响.通过数值算例分析结果与有限元方法所得结果进行比较,验证了本方法的正确性和有效性.     

柱体绕流的自适应数值模拟方法

发展了基于四叉树数据结构的网格生成和流动的Navier-Stokes方程数值求解器。采用压力梯度或者密度梯度的绝对值作为网格自适应的控制参量,同时采用基于最小二乘法的无网格方法处理对于一般Cartesian网格难于处理的物面边界条件。采取了绕方柱流动和绕圆柱流动的经典二维定常和非定常层流算例对所发展的方法进行了验证。计算的结果验证了所发展的方法在处理绕流流动时的合理性和有效性。从而也为将来数值模拟具有较复杂几何外形的流动提供了一种网格布局合理、高效,边界处理简单易行的新思路。     

Waspaloy高温合金涡轮盘锻造无网格法数值模拟

针对锻造工艺中的轴对称问题进行无网格伽辽金法数值模拟研究,基于Mises屈服准则和刚塑性流动理论给出轴对称问题的应力应变关系.通过移动最小二乘近似函数建立无网格伽辽金法的基本公式.在刚塑性本构关系和最小二乘法的基础上,基于马尔可夫变分原理建立无网格伽辽金法的求解列式.采用C++语言编制了轴对称无网格伽辽金法数值模拟程序.基于编制的程序进行了Waspaloy合金涡轮盘锻造过程的数值模拟,预测了锻件的最终形状和锻造过程中材料流动规律,获得了锻造过程中的节点和等效塑性应变分布规律.     

方柱绕流的大涡模拟

采用大涡模拟(LES)方法对雷诺数为2.2×104的方柱绕流流场进行了数值模拟.使用非交错网格的有限差分法,分别对准三维物理模型和真实三维物理模型求解不可压N-S方程,将沿流向方向方柱水平中心线上的时均速度的计算结果与实验数据进行了比较,结果表明,三维模型的模拟结果优于准三维模型的模拟结果.比较升力系数和阻力系数发现,与二维模拟(RNG)方法相比, 三维模拟的结果更加接近实验测试数据.     

用加权最小二乘无网格法求解抛物方程

加权最小二乘无网格法是一种新的高效无网格法,鉴于传统数值方法求解动态问题网格限制的缺陷,将传统差分法和加权最小二乘无网格法结合构造差分一加权最小二乘无网格方法,应用于求解一维与时间相关的线性抛物方程;该方法在空间域上的离散彻底摆脱了网格的束缚,算例表明:该方法计算量较小,并能够保证较高的精度.     

断裂力学的复变量无网格方法

在移动最小二乘法的基础上, 讨论了复变量移动最小二乘法. 复变量移动最小二乘法的优点是不形成病态方程组、精度高, 所形成的无网格方法计算量小. 利用裂纹尖端解析解将复变量移动最小二乘法的基函数进行扩展, 推导了相应的逼近函数; 从最小势能原理出发提出了断裂力学的复变量无网格方法, 推导了相应的复变量无网格方法的求解方程. 与传统的无网格方法相比, 断裂力学的复变量无网格方法具有计算量小、精度高的优点. 最后给出了数值算例.     

用加权最小二乘无网格法求解抛物方程

加权最小二乘无网格法是一种新的高效无网格法,鉴于传统数值方法求解动态问题网格限制的缺陷,将传统差分法和加权最小二乘无网格法结合构造差分一加权最小二乘无网格方法,应用于求解一维与时间相关的线性抛物方程;该方法在空间域上的离散彻底摆脱了网格的束缚,算例表明：该方法计算量较小,并能够保证较高的精度．     

三维最小二乘无网格方法精度与误差分析

为了得到最小二乘无网格方法在求解三维Euler方程时的计算精度及其应用于三维流场求解的准确性,通过数学方法对该无网格方法在正交等间距离散点分布情况下的计算精度进行了论证,给出了基于正交等间距和非正交均匀布点的4个三维验证算例.并将该文计算结果与其他网格方法得到的计算结果和理论结果进行了比较,两种方法得到的计算结果吻合较好,与理论结果的误差在合理范围以内.结果表明将最小二乘无网格方法应用于三维流场的求解是完全可行的.     

弹性动力学的边界无单元法

讨论了Hilbert空间上的改进的移动最小二乘法, 并将弹性动力学的边界积分方程方法与改进的移动最小二乘法结合, 提出了弹性动力学的边界无单元法. 改进的移动最小二乘法避免了求解病态方程组, 既提高了效率, 又提高了精度. 边界无单元法是边界积分方程的无网格方法的直接列式法, 容易引入边界条件, 且具有更高的精度. 最后给出了弹性动力学的边界无单元法的数值实现方法和具体的算例.     

自动网格体系在柱体绕流大涡模拟中的适用性评估

为了评估基于snappyHexMesh方法生成的自动网格体系在二维柱体绕流大涡模拟中的适用性,比较了该自动网格体系与人工网格体系对于Re为3 900圆柱绕流和Re为22 000方柱绕流的数值模拟结果。通过设置合理的计算域以及数值格式,采用snappyHexMesh自动网格以及人工网格的算例都表现出良好的数值稳定性。将不同网格体系的数值模拟结果与物理试验结果进行对比,结果表明,采用snappyHexMesh网格可以提高数值求解效率;圆柱绕流对网格体系的变化比较敏感,不同密度的snappyHexMesh网格会显著影响圆柱气动力特征以及尾流区域的流场结果;snappyHexMesh网格体系可以准确预测方柱绕流,在方柱绕流大涡模拟中具有相较于圆柱绕流更好的适用性。