Curvature‐compensated convective transport: SMART, A new boundedness‐ preserving transport algorithm

International Journal for Numerical Methods in Fluids - Tập 8 Số 6 - Trang 617-641 - 1988
Philip Gaskell1, A. K. C. Lau1
1Department of Mechanical Engineering, University of Leeds, Leeds LS2 9JT, U.K.

Tóm tắt

Abstract

The paper describes a new approach to approximating the convection term found in typical steady‐state transport equations. A polynomial‐based discretization scheme is constructed around a technique called ‘curvature compensation’; the resultant curvature‐compensated convective transport approximation is essentially third‐order accurate in regions of the solution domain where the concept of order is meaningful. In addition, in linear scalar transport problems it preserves the boundedness of solutions. Sharp changes in gradient in the dependent variable are handled particularly well. But above all, the scheme, when used in conjunction with an ADI pentadiagonal solver, is easy to implement with relatively low computational cost, representing an effective algorithm for the simulation of multi‐dimensional fluid flows. Two linear test problems, for the case of transport by pure convection, are employed in order to assess the merit of the method.

Từ khóa


Tài liệu tham khảo

Roache P. J., 1976, Computational Fluid Dynamics

Richtmyer R. D., 1967, Difference Methods for Initial‐Value Problems

10.1016/0045-7825(79)90034-3

B. P.Leonard ‘A survey of finite differences of opinion on numerical muddling of incomprehensible defective confusion equation’ paper in ASME Applied Mechanics Division Winter Annual Meeting 1979.

R. K.Agarwal ‘A third‐ordcr‐accurate upwind scheme for Navier‐Stokcs solution at high Reynolds numbers’ Paper No. AIAA‐81‐0112 AIAA 19th Aerospace Science Meeting St. Louis 1982.

B. P.Leonard ‘Adjusted quadratic upstream algorithms for transient incompressible convection’ AIAA‐ J. 226–233(1979).

10.1016/0021-9991(85)90188-3

P. H.GaskellandA. K. C.Lau ‘An efficient solution strategy for use with high‐order discretization schemes’ Report No. T40 Department of Mechanical Engineering University of Leeds 1986.

10.1016/0045-7825(76)90058-X

S. A.Syed A. D.GosmanandM.Peric ‘Assessment of discretization schemes to reduce numerical diffusion in the calculation of complex flows’ Paper No. AIAA‐85‐0441 AIAA 23rd Aerospace Science Meeting Reno January1985.

P. H.GaskellandA. K. C.Lau ‘An assessment of direct stress modelling for elliptic turbulent flows with the aid of a non‐diffusive boundedness preserving discretization scheme’ Proc. 5th int. Conf. on Numerical Methods in Laminar and Turbulent Flow Montreal 1987.

P. H.GaskellandN. G.Wright ‘Multigrids applied to an efficient fully coupled solution technique for recirculating fluid flow problems’ Proc. IMA Conf. on the Simulation and Optimisation of Large Systems Reading 1986.

10.2514/3.60623