A posteriori error estimation and adaptive mesh refinement for a multiscale operator decomposition approach to fluid–solid heat transfer

We analyze a multiscale operator decomposition finite element method for a conjugate heat transfer problem consisting of a fluid and a solid coupled through a common boundary. We derive accurate a posteriori error estimates that account for all sources of error, and in particular the transfer of error between fluid and solid domains. We use these estimates to guide adaptive mesh refinement. In addition, we provide compelling numerical evidence that the order of convergence of the operator decomposition method is limited by the accuracy of the transferred gradient information, and adapt a so-called boundary flux recovery method developed for elliptic problems in order to regain the optimal order of accuracy in an efficient manner. In an appendix, we provide an argument that explains the numerical results provided sufficient smoothness is assumed.     

Mixed Convection in Viscoelastic Boundary Layer Flow and Heat Transfer over a Stretching Sheet

An investigation is carried out on mixed convection boundary layer flow of an incompressible and electrically conducting viscoelastic fluid over a linearly stretching surface in which the heat transfer includes the effects of viscous dissipation, elastic deformation, thermal radiation, and non-uniform heat source/sink for two general types of non-isothermal boundary conditions. The governing partial differential equations for the fluid flow and temperature are reduced to a nonlinear system of ordinary differential equations which are solved analytically using the homotopy analysis method (HAM). Graphical and numerical demonstrations of the convergence of the HAM solutions are provided, and the effects of various parameters on the skin friction coefficient and wall heat transfer are tabulated. In addition, it is demonstrated that previously reported solutions of the thermal energy equation given in [1] do not converge at the boundary, and therefore, the boundary derivatives reported are not correct.     

Melting phenomenon in magneto hydro-dynamics steady flow and heat transfer over a moving surfacein the presence of thermal radiation

The Lie group method is applied to present an analysis of the magneto hydro-dynamics(MHD) steady laminar flow and the heat transfer from a warm laminar liquid flow to a melting moving surface in the presence of thermal radiation.By using the Lie group method,we have presented the transformation groups for the problem apart from the scaling group.The application of this method reduces the partial differential equations(PDEs) with their boundary conditions governing the flow and heat transfer to a system of nonlinear ordinary differential equations(ODEs) with appropriate boundary conditions.The resulting nonlinear system of ODEs is solved numerically using the implicit finite difference method(FDM).The local skin-friction coefficients and the local Nusselt numbers for different physical parameters are presented in a table.     

DNS of buoyancy-dominated turbulent flows on a bluff body using the immersed boundary method

A novel immersed boundary (IB) method has been developed for simulating multi-material heat transfer problem – a cylinder in a channel heated from below with mixed convection. The method is based on a second-order velocity/scalar reconstruction near the IB. A novel algorithm has been developed for the IB method to handle conjugate heat transfer. The fluid–solid interface is constructed as a collection of disjoint faces of control volumes associated to different material zones. Coupling conditions for the material zones have been developed such that continuity and conservation of the scalar flux are satisfied by a second-order interpolation. Predictions of the local Nusselt number on the cylinder surface show good agreement with the experimental data. The effect of the Boussinesq approximation on this problem was also investigated. Comparison with the variable density formulation suggests that, in spite of a small thermal expansion coefficient of water, the variable density formulation in a transitional flow with mixed convection is preferable.     

A Cartesian grid embedded boundary method for solving the Poisson and heat equations with discontinuous coefficients in three dimensions

We present a method for solving Poisson and heat equations with discontinuous coefficients in two- and three-dimensions. It uses a Cartesian cut-cell/embedded boundary method to represent the interface between materials, as described in Johansen and Colella (1998). Matching conditions across the interface are enforced using an approximation to fluxes at the boundary. Overall second order accuracy is achieved, as indicated by an array of tests using non-trivial interface geometries. Both the elliptic and heat solvers are shown to remain stable and efficient for material coefficient contrasts up to 10⁶, thanks in part to the use of geometric multigrid. A test of accuracy when adaptive mesh refinement capabilities are utilized is also performed. An example problem relevant to nuclear reactor core simulation is presented, demonstrating the ability of the method to solve problems with realistic physical parameters.     

Heat transfer and fluid flow of benard-cell convection in rectangular container with free surface sensed by infrared thermography

The natural convection flow phenomena that occur inside an enclosed space are very interesting examples of complex fluid systems that may yield to analytical, empirical and numerical solutions, and many reports have looked into this basic problem. In the present study, heat transfer and fluid flow for natural convection in a horizontal rectangular container with a free surface are investigated using infrared thermography. The temperature field was measured and visualized at a gas-liquid (air — silicon oil) interface using infrared thermography. The heat transfer phenomena were also investigated by statistically analyzing the temperature data. The applicability of the infrared thermography to quantitative heat transfer measurement at the gas-liquid interface was evaluated. It is revealed that infrared thermography is effective not only in visualization of a gas-liquid interface but also in heat transfer measurement. A new heat transfer correlation is proposed for the gas-liquid interface of this flow system. The coefficient of heat transfer can be summarized by a specific heat transfer correlation formula regardless of several conditions, including container aspect ratio, fluid viscosity and fluid layer depth.     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

Application of the Green functions to the problem of the thermally induced vibration of a beam

This paper considers the problem of the transverse vibrations of a beam induced by a mobile heat source. The formulation of the problem is based on the differential equations of heat conduction and transverse vibrations of the beam, which are complemented by suitable initial and boundary conditions. The effect of internal and external damping on the vibrations of the beam is considered. The solution to the problem in analytical form is obtained by using the properties of the Green functions. A time partitioning method has been used to avoid the difficulties associated with the slow convergence of the series occurring in the solution to the heat conduction problem. The numerical results of the thermally induced vibration of the beam are presented.     

A matched interface and boundary method for solving multi-flow Navier–Stokes equations with applications to geodynamics

We have developed a second-order numerical method, based on the matched interface and boundary (MIB) approach, to solve the Navier–Stokes equations with discontinuous viscosity and density on non-staggered Cartesian grids. We have derived for the first time the interface conditions for the intermediate velocity field and the pressure potential function that are introduced in the projection method. Differentiation of the velocity components on stencils across the interface is aided by the coupled fictitious velocity values, whose representations are solved by using the coupled velocity interface conditions. These fictitious values and the non-staggered grid allow a convenient and accurate approximation of the pressure and potential jump conditions. A compact finite difference method was adopted to explicitly compute the pressure derivatives at regular nodes to avoid the pressure–velocity decoupling. Numerical experiments verified the desired accuracy of the numerical method. Applications to geophysical problems demonstrated that the sharp pressure jumps on the clast-Newtonian matrix are accurately captured for various shear conditions, moderate viscosity contrasts and a wide range of density contrasts. We showed that large transfer errors will be introduced to the jumps of the pressure and the potential function in case of a large absolute difference of the viscosity across the interface; these errors will cause simulations to become unstable.     

The Linear KdV Equation with an Interface

The interface problem for the linear Korteweg–de Vries (KdV) equation in one-dimensional piecewise homogeneous domains is examined by constructing an explicit solution in each domain. The location of the interface is known and a number of compatibility conditions at the boundary are imposed. We provide an explicit characterization of sufficient interface conditions for the construction of a solution using Fokas’s Unified Transform Method. The problem and the method considered here extend that of earlier papers to problems with more than two spatial derivatives.     

The influence of an electric field on the heat transfer rate during film boiling of stagnant fluids

This paper describes the influence of an electric field on the heat transfer rate during film boiling of stagnant fluids (pool boiling) on a horizontal surface. Modelling the influence of an electric field on the heat transfer during film boiling requires the knowledge of the equilibrium shape of the liquid-vapour interface. This equilibrium shape is calculated using an iterative solution technique. In each iteration step a fourth order Runge-Kutta technique is used to calculate the shape of the interface and a boundary element method is used to calculate the electric field. It is shown that the vapour bubbles formed on this interface become elongated in the presence of an electric field. This effect is accounted for in the existing heat transfer models in order to describe the influence of the electric field on the heat transfer coefficient during film boiling.     

蒙特卡罗方法在解微分方程边值问题中的应用

介绍了蒙特卡罗方法的基本原理以及随机数的产生方法。基于蒙特卡罗方法的思想,结合有限差分方法,建立了求解微分方程边值问题的随机概率模型,并以第一类边界条件的拉普拉斯方程和一个给定初值及边界条件的非稳态热传导方程为数值算例,研究了蒙特卡罗方法在求解微分方程边值问题中的应用。结果表明:利用蒙特卡罗方法,不仅可以有效解决给定边界条件的微分方程,对于给定初值条件的微分方程,也可以从时域有限差分方程出发,采用蒙特卡罗方法进行求解。数值模拟和对误差的理论分析均表明,增加蒙特卡罗试验中的模拟粒子点数,可以提高计算结果的精度。     

Axisymmetric convection flow of fractional Maxwell fluid past a vertical cylinder with velocity slip and temperature jump

A novel finite volume method is developed to investigate the axisymmetric convection flow and heat transfer of fractional viscoelastic fluid past a vertical cylinder. Fractional cylindrical governing equations are formulated by fractional Maxwell model and generalized Fourier's law. The velocity slip and temperature jump boundary conditions are considered across the fluid-solid interface. Numerical results are validated by exact solutions of special case with source terms. The effects of fractional derivative parameter and boundary condition parameters on flow and heat transfer characteristics are discussed. The viscoelastic fluid performs evident shear thickening property in the fractional Maxwell constitutive relation. Moreover, the boundary condition parameters have remarkable influence on velocity and temperature distributions.     

Heat Transfe for Power Law Non-Newtonian Fluids

We present a theoretical analysis for heat transfer in power law non-Newtonian fluid by assuming that the thermal diffusivity is a function of temperature gradient. The laminar boundary layer energy equation is considered as an example to illustrate the application. It is shown that the boundary layer energy equation subject to the corresponding boundary conditions can be transformed to a boundary value problem of a nonlinear ordinary differential equation when similarity variables are introduced. Numerical solutions of the similarity energy equation are presented.     

Alternating Direction Implicit Orthogonal Spline Collocation on Non-Rectangular Regions

The alternating direction implicit (ADI) method is a highly efficient technique for solving multi-dimensional time dependent initial-boundary value problems on rectangles. When the ADI technique is coupled with orthogonal spline collocation (OSC) for discretization in space, we not only obtain the global solution efficiently, but the discretization error with respect to space variables can be of an arbitrarily high order. In [2], we used a Crank Nicolson ADI OSC method for solving general nonlinear parabolic problems with Robin's boundary conditions on rectangular polygons and demonstrated numerically the accuracy in various norms. A natural question that arises is: Does this method have an extension to non-rectangular regions? In this paper, we present a simple idea of how the ADI OSC technique can be extended to some such regions. Our approach depends on the transfer of Dirichlet boundary conditions in the solution of a two-point boundary value problem (TPBVP). We illustrate our idea for the solution of the heat equation on the unit disc using piecewise Hermite cubics.     

Stochastic radiative transfer in a finite plane-parallel medium with general boundary conditions

The time-independent radiative transfer problem in a scattering and absorbing planar random medium with general boundary conditions and internal energy source is considered. The medium is assumed to consist of two randomly mixed immiscible fluids, with the mixing statistics described as a two-state homogeneous Markov process. The problem is solved in terms of the solution of the corresponding free-source problem with simple boundary conditions which is solved using Pomraning-Eddington approximation in the deterministic case. A formalism, developed to treat radiative transfer in statistical mixtures, is used to obtain the ensemble-averaged solution. The average partial heat fluxes are calculated in terms of the albedoes of the source-free problem. Results are obtained for isotropic and anisotropic scattering for specular and diffused reflecting boundaries.     

Coupling heat conduction and radiative transfer

A method to compute combined heat conduction and radiative transfer in a semitransparent medium is proposed. The heat equation is solved by an implicit finite-difference scheme and the radiation is simulated using a ray-tracing method. Stability and convergence for the coupling of both methods is shown. The crucial point for ray-tracing as well as for discrete-ordinate methods is the selection of the transfer directions and their associated weights. Using a model problem different direction sets proposed in the literature are compared with respect to their accuracy and the computational effort.     

振荡流共轭换热格子——Boltzmann模拟

振荡流共轭换热现象广泛存在于热声热机等工程应用中.基于双分布格子-Boltzmann模型,对平行平板间振荡流共轭换热进行了数值模拟.通过假定共轭界面处流体和固体的未知内能分布函数均为对应的平衡态滑移修正格式,提出了一种处理共轭换热边界的新方法.模拟结果表明,该方法可以保证共轭界面上温度连续和热流连续.分析了不同流体与固体导热系数比情况下振荡流共轭换热的速度场、温度场以及热流分布的特点.     

Radiation transfer in an isotropically scattering homogeneous solid sphere

Radiative heat transfer in an absorbing, emitting, isotropically scattering sphere with a diffusively reflecting and emitting boundary is solved by the Galerkin method. By using the expressions given in this work, the radiation heat flux, the angular distribution of radiation intensity and the divergence of the radiation heat flux anywhere in the sphere can be determined with a high degree of accuracy for all values of the single scattering albedo, from small to moderately large optical radii. Results are presented for representative cases to illustrate the accuracy and the convergence of the method.     

Cubic Auto-Catalysis Reactions in Three-Dimensional Nanofluid Flow Considering Viscous and Joule Dissipations Under Thermal Jump

In this problem, simultaneous effects of Joule and viscous dissipationin three-dimensional flow of nanoliquid have been addressed in slip flow regime under timedependent rotational oscillations. Silver nanoparticles are submerged in the base fluid (water)due to their chemical and biological features. To increment the novelty, effects of cubicautocatalysis chemical reactions and radiative heat transfer have been incorporated in therelated boundary layer equations. Dimensionless partial differential system is solved byemploying the proposed implicit finite difference approach. Convergence conditions andstability criteria are obtained to ensure the convergence and accuracy of solutions.A comparative analysis is proposed for no-slip nanofluid flow (NSNF) and slip nanofluid flow(SNF). Variations in skin-friction coefficients, Sherwood and Nusselt numbers against physicalparameters are tabulated. It is investigated that velocity slip and temperature jump significantlycontrol drag forces and rate of heat transfer.