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.
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.