An algorithm for the simulation of thermally coupled low speed flow problems

Abstract

SUMMARY In this paper, we propose a computational algorithm for the solution of thermally coupled flows in subsonic regime. The formulation is based upon the compressible Navier–Stokes equations, written in nonconservation form. An efficient modular implementation is obtained by solving the energy equation separately and then using the computed temperature as a known value in the momentum‐continuity system. If an explicit single‐step time integration scheme for the energy equation is used, the decoupling results to be natural. Integration of the momentum‐continuity system is carried out using a semi‐explicit method, combining Runge–Kutta and Backward Euler schemes for the momentum and continuity equations, respectively. Implicit treatment of pressure leads to favorable time step estimates even in the low Mach number (Ma ≪ 1) regimes. The numerical dissipation introduced by the Backward Euler scheme ensures absence of the spurious high frequencies in the numerical solution. The key point of the method is the assumption of linear variation of the temperature within a time step. Combined with a fractional splitting of the momentum‐continuity system, it allows to solve the continuity only once per time step. Omitting the necessity of solving for the pressure at every intermediate step of the Runge–Kutta scheme minimizes the computational cost associated to the implicit step and leads to an efficiency close to that of a purely explicit scheme. The method is tested using two benchmark examples.Copyright © 2011 John Wiley & Sons, Ltd.

Lagrangian finite element model for the 3D simulation of glass forming processes

“…(1)) for the mechanical model. More details on the thermo-mechanical coupling can be found in [26], [27].…”

Finite Element Modeling of Free Surface Flow in Variable Porosity Media

“…where r i with i = 1, 2, 3, 4 are the residuals of the stationary form of Eq.36 evaluated at [48], a linear variation of pressure is assumed within the time step. It should be remarked that this assumption leads the velocity eld to be divergence-free only at the end of the step.…”

A portable OpenCL-based unstructured edge-based finite element Navier–Stokes solver on graphics hardware

“…In this scheme we solve explicitly the momentum equation (using a 4-step Runge-Kutta scheme) and implicitly the pressure correction step. A detailed discussion of the algorithm is presented in [7] for the incompressible case and in [8] for the solution of Low-Mach compressible problems.…”