Multi-dimensional finite volume scheme for the vorticity transport equations

被引:9
作者
Foti, Daniel [1 ]
Duraisamy, Karthik [1 ]
机构
[1] Univ Michigan, Dept Aerosp Engn, Ann Arbor, MI 48109 USA
关键词
Vorticity transport equation; Multi-dimensional upwind scheme; Wave propagation approach; ADAPTIVE MESH REFINEMENT; NAVIER-STOKES EQUATIONS; WAVE-PROPAGATION METHOD; CONSERVATION-LAWS; EULER EQUATIONS; SIMULATIONS; COMPUTATION; ALGORITHMS; ADVECTION; LIBRARY;
D O I
10.1016/j.compfluid.2018.02.024
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A finite-volume scheme is developed for the three-dimensional, incompressible vorticity transport equations (VTE) using multi-dimensional upwinding, with the goal of efficient computations of vortex-dominated flows. By modifying the VTE with a term proportional to the divergence of the vorticity, a stable hyperbolic PDE system with a simple structure is revealed. The structure of the resulting eigen-system, including the vortex stretching term, makes the formulation and solution of the generalized Riemann problem and multi-dimensional upwinding more natural, when compared to the Euler equations. To reduce the computational costs of determining the transverse fluxes, a flux-based wave propagation approach is employed. In this approach, the transverse fluxes are computed via direct manipulation of the one-dimensional generalized Riemann problem with no additional Riemann problem solutions needed. The numerical scheme is implemented within an adaptive mesh refinement framework and evaluated on a series of canonical vortex-dominated flows. A translating vortex flow reveals that the multi-dimensional upwinding substantially outperforms one-dimensional schemes in terms of accuracy and computational time, especially when the vortex propagates oblique to cell surfaces. By including transverse fluxes in simulations of propagating vortex rings, key physical attributes including propagation velocity and impulse can be captured with better accuracy compared to one-dimensional schemes. Further vortex ring simulations demonstrate that the proposed multi-dimensional scheme can preserve vorticity with relatively coarse grids compared to simulations employing the incompressible Euler equations. Integrated quantities from leapfrogging vortex ring problems and energy spectra from decaying turbulent flows confirm that the multi-dimensional scheme can accurately reproduce flows with complex vortex stretching and multiple length scales. (C) 2018 Elsevier Ltd. All rights reserved.
引用
收藏
页码:17 / 32
页数:16
相关论文
共 53 条
[11]   MULTIDIMENSIONAL UPWIND METHODS FOR HYPERBOLIC CONSERVATION-LAWS [J].
COLELLA, P .
JOURNAL OF COMPUTATIONAL PHYSICS, 1990, 87 (01) :171-200
[12]   A novel velocity-vorticity formulation of the Navier-Stokes equations with applications to boundary layer disturbance evolution [J].
Davies, C ;
Carpenter, PW .
JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 172 (01) :119-165
[13]   High resolution wake capturing methodology for hovering rotors [J].
Duraisamy, Karthikeyan ;
Baeder, James D. .
JOURNAL OF THE AMERICAN HELICOPTER SOCIETY, 2007, 52 (02) :110-122
[14]   LOCAL REFINEMENT TECHNIQUES FOR ELLIPTIC PROBLEMS ON CELL-CENTERED GRIDS .1. ERROR ANALYSIS [J].
EWING, RE ;
LAZAROV, RD ;
VASSILEVSKI, PS .
MATHEMATICS OF COMPUTATION, 1991, 56 (194) :437-461
[15]  
Falgout RD, 2002, LECT NOTES COMPUT SC, V2331, P632
[16]   Multidimensional upwinding. Part II. Decomposition of the Euler equations into advection equations [J].
Fey, M .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 143 (01) :181-199
[17]   Wake meandering statistics of a model wind turbine: Insights gained by large eddy simulations [J].
Foti, Daniel ;
Yang, Xiaolei ;
Guala, Michele ;
Sotiropoulos, Fotis .
PHYSICAL REVIEW FLUIDS, 2016, 1 (04)
[18]  
GODUNOV S. K., 1972, Numer. Methods Mech. Continuum Medium, V1, P26
[19]  
Greenshields CJ., 2015, OpenFOAM Foundation Ltd
[20]  
Hariharan N, 2014, P 52 AER SCI M, P2014