A new efficient and accurate procedure for solving heat condution problems

In this paper, a new technique of differential quadature method called the upwind difference - differential quadature method (UDDQM) for solving two-dimensional heat transfer (convection-diffusion) problems is proposed. Also, investigated the effects of physical quantities on behavior of flow problems, and combined effects of upwind difference mechanism together with differential quadrature method to modified the numerical solutions of heat transfer problems are presented. To validate our proposed UDDQM, two convection-diffusion problems ((i) Steady-state incompressible flow problem has exact solution and (ii) Natural convection motion of the incompressible fluid flow problem hasn't exact solution) are solving numerically. Graphical results on the effects of parameter variation on velocity, temperature, Peclet number, Grashof number, and Prandtl number are presented and discussed. Numerical experiments are conducted to test its accuracy and convergence and compare it with the standard DQM and other numerical methods that are available in literature. The numerical results show the efficiency of the proposed method to handle the problems, and it is more accurate and convergent than other methods.

Natural convection in a square cavity containing a nanofluid and an adiabatic square block at the center

Numerical evaluation of the latent heat thermal energy storage performance enhancement by installing longitudinal fins

Preface. Solving heat transfer problems by the dual reciprocity BEM. Transient problems using time-dependent fundamental solutions. Solving linear heat conduction problems by the multiple reciprocity method. Solving nonlinear heat transfer problems using the boundary element method. Coupled conduction-convection problems. Solving coupled problems involving conduction, convection and thermal radiation. Advanced thermoelastic analysis. Integral equation analyses of natural convection problems in fluid flow. Improperly posed problems in heat transfer. Acknowledgments. References.

Four problems of fluid flow and heat transfer were designed in which non‐orthogonal, boundary‐fitted grids were to be used. These are proposed to serve as test cases for testing new solution methods. This paper presents solutions of the test problems obtained by using a multigrid finite volume method with grids of up to 320 × 320 control volumes. Starting from zero fields, iterations were performed until the sum of the absolute residuals had fallen seven orders of magnitude, thus ensuring that the variable values did not change to six most significant digits. By comparing the solutions for successive grids at moderate Reynolds and Rayleigh numbers, the discretization errors were estimated to be lower than 0·1%. The results presented in this paper may thus serve for comparison purposes as bench‐mark solutions.

The solution of the conjugate problem of convective-conductive heat transfer in the channels of power plants is presented. The problem of convective heat transfer in the gas phase is solved by the integral theory of heat transfer. A one-dimensional problem of conductive heat transfer in the material of the wall is solved by the finite-difference method. Relative laws of heat and mass transfer and friction are obtained by numerical integration with respect to the boundary-layer thickness. The effect of the material and geometry of the wall on the level of problem “conjugation” is studied.

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

The main goal of this paper is to demonstrate the general characteristics of the educational
user-friendly CFD Studio package for CFD teaching. The package was designed for
teaching 2D fluid mechanics and heat transfer process, including conduction, coupled
conduction/convection, natural and forced convection, external and internal flows, among
other phenomena. The finite volume methodology and its related topics can also be taught
using the software. Therefore, general aspects of the three main modules, pre-processor,
solver and post-processor are discussed aiming to show the generality of the tool. These
modules are integrated in the application by a so-called “numerical problem project”
which guide the student through the steps to obtain the solution. To approximate the
partial differential equations the finite volume approach is employed using a fully-implicit
formulation with the interpolation schemes CDS, UDS and WUDS. Mesh editing and nonorthogonal
boundary-fitted mesh generation, using algebraic interpolation and elliptic
equations, are important features of the package. Coupled heat transfer problems are
handled using the “solid-block” formulation and the pressure-velocity coupling uses the
SIMPLE and SIMPLEC methods with non-staggered grids. To demonstrate the
capabilities two fluid flow and heat transfer “problem projects” are presented.

A primary goal of the Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is to develop the next generation of simulations tools to support the design and deployment of commercial nuclear power. The mission of the NEAMS Center of Excellence for Thermal-fluids Applications in Nuclear Energy is to advance this goal by providing leadership, best practices, research, and support and training for computational thermal-fluids applications in nuclear energy. This is accomplished by leveraging expertise in the NEAMS tools to solve challenging problems in fluid flow and heat transfer. The NEAMS tools provide simulation capabilities covering a range of temporal and spatial scales and the Center connects stakeholders in industry with the expertise to use these tools to their fullest. This addresses a pressing need, facilitating advanced reactor development and commercialization.

The problem of non-Newtonian fluid flow has taken considerable interest and has been the subject of several work in latest years due to its various requests in different fields of engineering, in particular the interest in the problems of heat transfer in non-Newtonian liquids, such as lubrication, hot rolling, cooling problem and drag reduction. Here, mixed convection heat transport and its related entropy production in a porous channel with circular obstacle saturated via non-Newtonian power law liquid has been scrutinized. The influences on entropy production of the power law index, the Reynolds number, the Rayleigh number and the Darcy number are investigated. Being a novelty of this work, an optimization study of the thermodynamic irreversibility as a function of the channel inclination angle and the power law index is undertaken. The governing equations of the problem are solved employing the COMSOL software. Outcomes illustrate that the governing parameters strongly affect the entropy production. The thermal entropy generation is maximal at low value of power law index and high value of Reynolds number. The effect of Reynolds number become insignificant at relatively high power law index. At fixed Reynolds number value, a rise in the power index (n) leads to a reduce in the thermal entropy. This decrease is tiny, at low value of Reynolds number (Re) and turn into increasingly considerable as Re rises. The streamlines show the existence of two recirculation zones just after the circular obstacle, whose existence depends on both Re and power law index. Results show that the greatest variation relating to the inclination angle is for power law index equal to 0.4. Results indicate also that, at low Darcy number and relatively high power law index, the intrinsic effect of the modified Darcy number on Darcy viscous irreversibility become pronounced giving a sharp increase in the total entropy production.

A general procedure for infrared thermography heat transfer measurements in hypersonic wind tunnels

This paper explores the application of the mesh-free Local Radial Basis Function Collocation Method (LRBFCM) [1] in solution of coupled heat transfer and fluid flow problems associated with solidification of a binary alloy. The involved temperature, velocity, species and pressure fields are represented on overlapping sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBF’s. The energy and momentum equations are solved through explicit time stepping. The pressure-velocity coupling is calculated iteratively, with pressure correction, predicted from the local continuity equation violation [2,3]. The solution procedure is demonstrated in simulation of solidification of 8 wt % Pb-Sn alloy and 5 wt % Sn-Pb alloy in rectangular cavity at conditions, defined in [4]. The numerical examples include studies with different uniform discretization and its influence on final macrosegregation map. The solution is assessed by comparison with the reference results of the finite volume and finite element method. The advantages of the method are simplicity, accuracy and straightforward applicability in non-uniform node arrangements. The methodology will be use in complementing our heat transfer model [5] of the continuous casting of steel.

Errors in one-dimensional fin optimization problem for convective heat transfer

In this work an empirical estimator is used to estimate the iteration error based on the convergence rate of the variable of interest. Problems of heat transfer and of fluid mechanics are solved by the finite-difference and finite-volume methods using various iteration methods. In the initial iterations the accuracy of the empirical estimator is usually low; when the number of iterations is high, round-off errors predominate over iteration errors, but even so, the accuracy is relatively good; and in the interval between these two extremes, the accuracy tends to improve as the number of iterations increases.

We couple a node-centered finite volume method to a high order finite difference method to simulate dynamic earthquake ruptures along nonplanar faults in two dimensions. The finite volume method is implemented on an unstructured mesh, providing the ability to handle complex geometries. The geometric complexities are limited to a small portion of the overall domain and elsewhere the high order finite difference method is used, enhancing efficiency. Both the finite volume and finite difference methods are in summation-by-parts form. Interface conditions coupling the numerical solution across physical interfaces like faults, and computational ones between structured and unstructured meshes, are enforced weakly using the simultaneous-approximation-term technique. The fault interface condition, or friction law, provides a nonlinear relation between fields on the two sides of the fault, and allows for the particle velocity field to be discontinuous across it. Stability is proved by deriving energy estimates; stability, accuracy, and efficiency of the hybrid method are confirmed with several computational experiments. The capabilities of the method are demonstrated by simulating an earthquake rupture propagating along the margins of a volcanic plug.