Full Article - Open Access.

Idioma principal


Kazemi-Kamyab*, V.; Zuijlen, A.H. van; Bijl, H.;

Full Article:

Thermal interaction of fluids and solids, or conjugate heat transfer (CHT), is en- countered in many engineering applications. Since time-accurate computations of such cou- pled problems can be computationally expensive, we consider loosely-coupled and strongly- coupled solution algorithms in which higher order multi-stage Runge-Kutta schemes are em- ployed for time integration. The higher order time integration schemes have the potential to improve the computational efficiency at arriving at a certain accuracy relative to the tradition- ally used 1st and 2nd order implicit schemes. The spatial coupling between the subdomains is realized using Dirichlet-Neumann interface conditions and the coupled domains are solved in a sequential manner at each stage (Block Gauss-Seidel). In this paper, the stability of two par- titioned algorithms is analyzed by considering a one dimensional model problem. The model problem consists of two thermally coupled domains where the governing equation within each subdomain is unsteady linear heat conduction. In the loosely-coupled approach, a family of multi-stage IMEX schemes is used for time integration. By observing similarities between the second stage of the IMEX schemes and the θ scheme with θ = 0.5 (Crank-Nicolson), the stability of the partitioned algorithm in which the Crank-Nicolson scheme is used for time integration is first analyzed by applying the stability theory of Godunov-Ryabenkii. The stability of the IMEX schemes is next investigated by numerically solving the model problem and comparing the results to the conclusions of the stability analysis for the Crank-Nicolson scheme. Due to partly explicit nature of the IMEX schemes, the loosely-coupled algorithm becomes unstable for sufficiently large Fourier numbers (similar to the Crank-Nicolson scheme). When the ratio of the thermal effusivities of the coupled domains is much smaller than unity, time step restriction due to stability is sufficiently weak that computations can be performed with reasonably large Fourier numbers. Furthermore, the results show better stability properties of the IMEX schemes compared to the Crank-Nicolson scheme. In the strongly-coupled approach, the stability and rate of convergence of performing (Gauss-Seidel) subiterations at each stage of the higher order implicit ESDIRK time integra- tion schemes are analyzed. From the stability analysis, an expression for the rate of conver- gence of the iterations (σ) is obtained. For cases where σ ≪ 1, subiterations will convergence rapidly. However, when σ ≈ 1, the convergence rate of the iterations is slow. The results ob- tained by solving the model problem numerically are in line with the performed analytical stability analysis.

Full Article:

Palavras-chave: High order implicit time integration, Conjugate heat, Partitioned algorithm, Sta- bility.,


DOI: 10.5151/meceng-wccm2012-19084

Referências bibliográficas
  • [1] Henshaw W. D., Chand K. K., “A composite grid solver for conjugate heat transfer in fluid-structure systems”. Journal of Computational Physics 228, 3708-3741, 2009.
  • [2] van Zuijlen A. H., Bijl H., “Implicit and explicit higher order time integration schemes for structural dynamics and fluid-structure interaction computations”. Computers and Structures 83, 93-105, 2005.
  • [3] Giles M. B., “Stability analysis of numerical interface conditions in fluid-structure thermal analysis”. Int. J. Numer. Meth. Fluids 25, 421-436, 1997.
  • [4] Kennedy C. A., Carpenter M. H., “Additive Runge-Kutta schemes for convectiondiffusion- reaction equations”. Applied Numerical Mathematics 44, 139-181, 2003.
  • [5] Kazemi-Kamyb V., van Zuijlen A. H., Bijl H., “Accuracy and stability analysis of a second order time-accurate loosely-coupled partitioned algorithm for transient conjugate heat transfer problems”. Submitted to Comput. Methods Appl. Mech. Engrg, 2011.
  • [6] Roe B., Jaiman R., Haselbacher A., Geubelle P. H., “Combined interface boundary condition method for coupled thermal simulations”. Int. J. Numer. Meth. Fluids 57, 329-354, 2008.
Como citar:

Kazemi-Kamyab*, V.; Zuijlen, A.H. van; Bijl, H.; "ON THE STABILITY OF LOOSE AND STRONG PARTITIONED ALGORITHMS FOR THERMAL COUPLING OF DOMAINS USING HIGHER ORDER IMPLICIT TIME INTEGRATION SCHEMES", p. 2954-2968 . In: In Proceedings of the 10th World Congress on Computational Mechanics [= Blucher Mechanical Engineering Proceedings, v. 1, n. 1]. São Paulo: Blucher, 2014.
ISSN 2358-0828, DOI 10.5151/meceng-wccm2012-19084

últimos 30 dias | último ano | desde a publicação