用于三维非均质流场计算的改进流线迎风Petrov-Galerkin(SUPG)方法

苏波 1; 韩向科 2

doi:10.7652/xjtuxb201403023

您当前的位置：

首页 >

文章列表页 >

用于三维非均质流场计算的改进流线迎风Petrov-Galerkin(SUPG)方法

更新时间：2025-04-29

- 用于三维非均质流场计算的改进流线迎风Petrov-Galerkin(SUPG)方法
- Modified Streamline Upwind Petrov-Galerkin(SUPG)Strategy for Calculation of Three-Dimensional Heterogeneous Flow Problems
- 西安交通大学学报 2014年48卷第3期页码：128-134
- 作者机构：
  
  江苏大学土木工程与力学学院,江苏,镇江,212013
- 作者简介：
- 基金信息：
  
  : 2013-08-05。作者简介: 苏波(1977-),男,博士,讲师。基金项目: 国家自然科学基金青年基金资助项目(51108210);江苏省博士后基金资助项目(1301048C);国家自然科学基金专项数学天元基金资助项目(11226308)
- DOI：10.7652/xjtuxb201403023
  中图分类号： O357.1
- 网络首发：2014-03-10，
  
  纸质出版：2014
- 稿件说明：
移动端阅览
苏波 1, 韩向科 2. 用于三维非均质流场计算的改进流线迎风Petrov-Galerkin(SUPG)方法[J]. 西安交通大学学报, 2014,48(3):128-134.

Modified Streamline Upwind Petrov-Galerkin(SUPG)Strategy for Calculation of Three-Dimensional Heterogeneous Flow Problems[J]. 2014, 48(3): 128-134.
苏波 1, 韩向科 2. 用于三维非均质流场计算的改进流线迎风Petrov-Galerkin(SUPG)方法[J]. 西安交通大学学报, 2014,48(3):128-134. DOI： 10.7652/xjtuxb201403023.

Modified Streamline Upwind Petrov-Galerkin(SUPG)Strategy for Calculation of Three-Dimensional Heterogeneous Flow Problems[J]. 2014, 48(3): 128-134. DOI： 10.7652/xjtuxb201403023.

摘要

在流线迎风Petrov-Galerkin(SUPG)稳定有限元方法基础上

通过引入双时间步法和变量分裂算法

发展了一种可用于三维非均质流场计算的改进SUPG方法。该方法摒弃了传统不可压缩流动问题中密度为常数的假定

采用包含密度输运方程的流体运动控制方程; 基于变量分裂算法

速度、压强场采用同阶插值函数进行空间离散

使改进SUPG方法具有简明的有限元格式

同时对速度场、压强场进行迭代求解

以降低线性代数方程组的阶数。双时间步法的引入

有利于提高SUPG方法对复杂非定常问题的求解效率。采用该方法对非均质、非定常三维矩形管道重力作用下的自由流动问题进行了分析

研究了重力作用下两种不同密度液体的相对运动过程。计算分析表明:在采用较大时间步的情况下

速度场和压强场在整个流动过程中随时间平稳过渡且分布光滑

没有出现数值波动现象; 旋涡位置及其随时间变化的规律与经典文献结果相符

没有出现跳跃和不连续现象。算例分析表明

改进SUPG方法具有良好的计算精度及数值稳定性

可用于三维非均质流动类似问题的研究。

Abstract

By introducing the dual-time stepping method and the variable splitting algorithm

a modified SUPG strategy is developed based on the stable streamline upwind Petrov-Galerkin(SUPG)finite element method

where the traditional assumption of constant density for incompressible flow problem is abandoned and the density transportation equation is introduced into the governing equations. The velocity and pressure fields are discretized with interpolating function of the same order

thus finite element scheme of the modified SUPG strategy gets simple and clear

and the order of algebraic equations is reduced. The dual-time stepping method is also introduced to enhance calculation stability of the SUPG strategy for complex unsteady problems. A free flow with heterogeneous

unsteady field of three-dimensional rectangular pipe under gravity and the whole relative motion between two types of liquids with different density are analyzed. The calculation results indicate that the velocity and pressure fields distribute and transmit smoothly with time and no numerical wave occurs in the case of greater time steps; the vortex position and its varying regulation coincide well with the results in the classic literatures without jumps and discontinuities. Example proves the numerical stability and accuracy of the modified SUPG method.

关键词

Keywords

references

HUANG Cheng, ZHOU Dai, BAO Yan, et al. A stabilized finite element technique and its application for turbulent flow with high Reynolds number[J]. Wind and Structures, 2011, 14(5): 465-480.

ZIENKIEWICZ O C, TAYLOR R L. The finite element method for fluid dynamics[M]. 5th ed. Amsterdam, The Netherlands: Elsevier, 2000.

钱若军, 董石麟, 袁行飞. 流固耦合理论研究进展[J]. 空间结构, 2008, 4(1): 3-15.

QIAN Ruo-jun, DONG Shi-lin, YUAN Xing-fei. Advances in research on fluid-structure interaction theory[J]. Spacial Structures, 2008, 4(1): 3-15.

韩向科, 苏波. 黏性不可压流体的一种FCBIS三角形单元[J]. 力学与实践, 2010, 32(3): 22-25.

HAN Xiangke, SU Bo. A FCBIS triangular element for incompressible fluid flows[J]. Mechanics in Engineering, 2010, 32(3): 22-25.

KOHNO H, BATHE K J. A nine-node quadrilateral FCBI element for incompressible fluid flows[J]. Communications in Numerical Methods in Engineering, 2006, 22(8): 917-931.

BROOKS A N, HUGHES T J R. Streamline upwind/Petrov-Galerkin formulation for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations[J]. Computer Methods in Applied Mechanics and Engineering, 1982, 32(1/2/3): 199-259.

JAMESON A. Time dependent calculation using multigrid with applications to unsteady flows past airfoils and wings[C]∥Proceedings of AIAA 10th Computational Fluid Dynamics Conference. Reston, VA, USA: AIAA, 1991: 1-8.

NITHIARASU P, CODINA R, ZIENKIEWICZ O C. The characteristic-based split(CBS)scheme: a unified approach to fluid dynamics[J]. Numerical Methods in Engineering, 2006, 66(10): 1514-1546.

ADINA RD Inc. ADINA CFDFSI: theory and modeling guide[M]. Watertown, MA, USA: ADINA RD Inc., 2005.

DONEA J, HUERTA A. Finite element methods for flow problems[M]. London, UK: John Wiley Sons, 2003: 92-93.

GRESHO P M, SANI R L. Incompressible flow and the finite element method[M]. New York, USA: John Wiley Sons, 2000: 302.

TRYGGVASON G. Numerical simulation of the Rayleigh-Taylor instability[J]. Journal of Computational Physics, 1988, 75(2): 253-282.

CALGARO C, CREUSE E, GOUDON T. A hybrid finite volume-finite element method for variable density incompressible flows[J]. Journal of Computational Physics, 2008, 227(9): 4671-4696.

GUERMOND J L, QUARTAPELLE L. A projection FEM for variable density incompressible flows[J]. Journal of Computational Physics, 2000, 165(1): 167-188.

GUERMOND J L, SALGADO A. A splitting method for incompressible flows with variable density based on a pressure Poisson equation[J]. Journal of Computational Physics, 2009, 228(8): 2834-2846.

赵立波,徐龙起,热合曼艾比布力,等.矩形微悬臂梁的流固耦合谐振频率分析.2013,47(11):60-64.[doi:10.7652/xjtuxb201311011]

何联格,左正兴,向建华.气缸盖中两相流沸腾换热热机耦合仿真分析.2013,47(7):23-28.[doi:10.7652/xjtuxb201307 005]

李连忠,牛文全,魏正英,等.微灌压力调节器调压特性流固耦合数值计算及分析.2013,47(5):131-136.[doi:10.7652/xjtuxb201305024]

何联格,左正兴,向建华.气缸盖冷却水腔内两相流动沸腾传热仿真研究.2013,47(1):21-26.[doi:10.7652/xjtuxb 201301005]

张颖莉,种道彤,刘继平,等.方管内混合对流与管壁导热耦合换热的数值模拟.2012,46(5):19-24.[doi:10.7652/xjtuxb201205004]

刘小民,王星,许运宾.运动罐体内液体晃动的双向流固耦合数值分析.2012,46(5):120-126.[doi:10.7652/xjtuxb201205 021]

张慧贤,寇子明,吴娟,等.液压激波作用下管道流固耦合的动力学建模.2012,46(3):94-99.[doi:10.7652/xjtuxb201203 017]

康伟,张家忠,李凯伦.利用本征正交分解的非线性Galerkin降维方法.2011,45(11):58-62.[doi:10.7652/xjtuxb201111 011]

梅冠华,张家忠.时滞惯性流形在三维壁板颤振数值分析中的应用.2011,45(9):40-46.[doi:10.7652/xjtuxb201109008]

苏波,钱若军,韩向科.一种用于流固耦合分析的有限元网格简捷更新方法.2011,45(3):16-24.[doi:10.7652/xjtuxb 201103003]

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

用于多导体传输参数提取的等效电流密度法

叶轮机械通道内压力的修正方程和数值模拟

改进的基于特征线理论的流体力学有限元法

膨胀管膨胀参数优化和膨胀模拟

基于微动滑移模型的叶片阻尼结构参数研究