This work aims to clarify and discuss, in a simultaneously accurate and simple way, some relevant issues on the stability and convergence of numerical solutions that are not usually presented in the literature and available for the students in this form. These include stability and physically realistic solutions of unsteady problems, and why the unconditionally stable Crank-Nicolson scheme can lead

In this article, we developed a computationally efficient multilevel local radial basis function (RBF-FD) mesh-free algorithm. The algorithm provides a new strategy to get good order of accuracy with less computational time, which is most important in the present world. The main idea is the layer-by-layer calculation and then layer-by-layer correction from coarsest level to finest level node points

Natural convection in an inclined cavity filled with saturated porous media has been investigated. A time-dependent solution using a spectral weighted residual method is determined. Preparing the governing equations by the substitution of the stream functions, we obtained a system of nonlinear ordinary differential equations. Integrating in time using the Runge-Kutta method of variable order the final

Turbulent combustion models are substantially important in CFD simulation of combustion. Some simple models frequently cannot well simulate the finite-rate chemistry. Widely recognized models are good only for certain flame types and flame structures. More general is the PDF equation model, but its computation requirement is very large, in particular when using large-eddy simulation (LES). Alternatively

A Galerkin-based finite element recursion relation is used to solve the heat transport equation in two-dimensions. The finite element method (FEM) is a powerful technique that is commonly used for solving complex engineering problems. However, the implementation of the FEM in multi-dimensional problems can be computationally expensive. A finite element recursion algorithm based on bilinear triangular

Numerical simulations were carried out to study the heat transfer and friction characteristics for Stirling engine heater tubes with three-dimensional internal extended micro-rib. During the numerical simulations, phase angle in a cycle ranged from 0° to 360°. The results represent that the friction factor (f) increased with micro-rib length (L) and micro-rib height (H) for enhanced tube within the

A two-level variational multiscale meshless local Petrov-Galerkin (VMS-MLPG) method is presented for incompressible Navier-Stokes equations based on two local Gauss integrations which effectively replace a unit operator (first level) and an orthogonal project operator (second level). The present VMS-MLPG method allows arbitrary combinations of interpolation functions for the velocity and pressure fields

In this article, a new modified skew upwind differencing scheme (MSUDS) for convective flow is presented. This scheme is tested by pure advection test cases such as Step, Sine-square, Semi-ellipse, and Smith-Hutton test profiles. The benchmark lid-driven cavity flow problem is also solved with this MSUDS scheme. All the results are well plotted and validated. All the data regarding overshoots/undershoots

A Compressive volume-of-fluid (VOF) schemes exhibits numerical diffusion which inhibit them in obtaining a numerically sharp and wrinkle free description of fluid interface which is vital for understanding the complex interfacial dynamics. Therefore, the present study introduces a novel compressive VOF scheme capable of capturing sharp abrupt interfaces at stringent Courant conditions while demonstrating

In this article, a finite volume formulation for solving the Navier-Stokes equations using unstructured hybrid grids and a staggered arrangement of variable pressure and velocity is presented. In this manner, a tight spatial pressure-velocity coupling is ensured without compromising geometrical flexibility. A second contribution of this work lays upon proposing a suitable interpolation function for

In this article, we propose a finite volume scheme preserving the discrete maximum principle (DMP) for steady heat conduction equations on distorted meshes. In contrary to these finite volume schemes preserving DMP, our new scheme uses the geometric average (instead of harmonic average) of two one-side numerical heat fluxes, especially it produces a more accurate flux approximation, which is verified

Flow pattern and heat transfer of flow boiling for different flow orientation, mini-channel width and height were presented in this work. The data were obtained by the numerical simulation with the coolant of R141b flow in a vertical mini-channel. Orientation includes upward and downward. A constant heat flux was loaded at the wall of the channel, of which the width ranges from 1 mm to 3 mm, and a

A novel formulation for numerical solution of heat conduction problems using superposition of exact solutions (SES) to represent temperature on sub-elements of a region is described and demonstrated. A simple 1-D linear problem is used to describe the method and highlight potential benefits; however, extensions to higher geometric dimensions and linearization for temperature-dependent properties are

This article presents a novel numerical method for steady-state thermal simulation. This method firstly solves the heat flux efficiently by applying the loop-tree basis functions. Then, the temperature is obtained by finding solutions of the gradient equation. The half boundary Rao-Wilton-Glisson (HBRWG) basis functions are employed for handling arbitrary boundary conditions. In addition, the triangulation-based

The article deals with an implicit formulation of the pressure far field boundary condition, also known as the characteristic boundary condition, in a pressure-based coupled solver. This boundary condition applies to compressible flows over the entire Mach regime, and is derived by invoking the Riemann invariants to implicitly express the flow variables at inlets and outlets in terms of their values

A continuum mixture model is proposed for phase change of binary alloy including the effect of thermal anisotropy. Thermal anisotropy is incorporated by an additional departure source in the conventional isotropic heat transfer-based energy equation. The governing equations for heat, mass, momentum, and species transport are solved using implicit finite volume method. The pressure and velocity in the

The SIMPLE algorithm is devised by interpolating the mass continuity and non-advective momentum equations, provoking apparent simplicity and clarity in the formulation. The SIMPLE variant entitled as the SIMPLE-AC scheme convokes an artificial compressibility (AC) parameter to augment the diagonal dominance of discretized pressure-correction equation. Both methods are characteristically pressure-based

The non-Newtonian fluids presenting viscoelastic flow behavior are found in many engineering applications. The development of a new numerical scheme for solution of this class of problems is the main goal of the present work. The proposed methodology adopts a second-order fully implicit finite difference approximation to discretize the convection and diffusion terms in the governing equations. Besides

In this article, the high-order upwinding combined compact difference scheme developed in a three-point grid stencil is applied to solve the incompressible Navier-Stokes (NS) and energy equations in three dimensions. The time integrator with symplectic property is employed to approximate the temporal derivative term in inviscid Euler equation so as to numerically retain the embedded Hamiltions and

In the article, we solve the inverse problems to recover unknown space-time dependent functions of heat conductivity and heat source for a nonlinear convective-diffusive equation, without needing of initial temperature, final time temperature, and internal temperature data. After adopting a homogenization technique, a set of spatial boundary functions are derived, which satisfy the homogeneous boundary

This article deals with the development of an implicit and conservative method for conjugate heat transfer at solid-fluid interfaces. The technique is applicable for both conformal and non-conformal meshes. The method, which is implemented within a fully coupled in-house code, is symmetric in its treatment of the solid and fluid regions and is shown to be very robust for highly complex configurations

We present a novel hybrid scheme for the large eddy simulation (LES) of turbulent reacting flows. The scheme couples the discontinuous spectral element method (DSEM) solver for the unsteady compressible Navier-Stokes equations with a Monte Carlo particle filtered mass density function (FMDF) solver for the transport of reacting species. The method is capable of high-order simulations on unstructured

Unsteady laminar and turbulent pipe flows subjected to a constant temperature gradient in axial direction are investigated based on asymptotic considerations. The thermal boundary condition is that of a constant wall heat flux for a steady flow. The unsteadiness is induced by a sinusoidal pressure gradient with different amplitudes and frequencies. The asymptotic analysis is performed with the help

In the first part of this two-paper series, a new computational approach is presented for analyzing transient heat conduction problems in anisotropic nonhomogeneous media. The approach consists of a truly meshless Fragile Points Method (FPM) being utilized for spatial discretization, and a Local Variational Iteration (LVI) scheme for time discretization. In the present article, extensive numerical

In this article, a grad-div stabilized projection finite element method is proposed for the double-diffusive natural convection model. Moreover, for the presented method, the stability analysis, and error estimates are deduced. Finally, numerical tests are provided that demonstrate the efficiency of the method. It is found that the presented method can improve the accuracy of the solution compared

A new and effective computational approach is presented for analyzing transient heat conduction problems. The approach consists of a meshless Fragile Points Method (FPM) being utilized for spatial discretization, and a Local Variational Iteration (LVI) scheme for time discretization. Anisotropy and nonhomogeneity do not give rise to any difficulties in the present implementation. The meshless FPM is

Incompressible two-phase flow involving interface evolution was studied using a level set method. To maintain the level set function as a signed distance function from the zero level set and meanwhile preserve the mass conservation, a level set redistancing algorithm of level set methods was developed. Important to the above level set redistancing algorithm was a new idea that keeps the zero level

In this work, one three-dimensional periodic composite structure consisting of alumina frame and inscribed nickel spheres is proposed for designing the micron composite porous structure. The finite-difference time domain method is applied for predicting the radiation spectra of proposed composite structure and distribution of absorbed radiative power. The spectral characteristics of composite structure

In this article, a completely new numerical method called the Local Least-Squares Element Differential Method (LSEDM), is proposed for solving general engineering problems governed by second order partial differential equations. The method is a type of strong-form finite element method. In this method, a set of differential formulations of the isoparametric elements with respect to global coordinates

This work proposes a mechanical analog of the unsteady energy conservation equation for analysis of the stability of its numerical solution. The space discretized energy conservation equation is the analog of the linear momentum equation, from which no relevant information can be extracted concerning the stability of its numerical solution. The corresponding angular momentum equation can be obtained

One of the main factors affecting the heat transfer efficiency of solar collector is that the ordinary fluid in it is in the state of natural convection. Supercritical fluid is expected to improve the heat transfer efficiency of solar collectors due to its dramatic changes in thermal properties, so it is necessary to carry out the research of natural convection of supercritical fluid. Although researchers

In this article, a new interpolation method is developed for the Lagrangian statistics in the non-isothermal gas-solid flow. The Lagrangian statistics are inseparable from interpolation, especially the optimal interpolation. The new interpolation method in this article is developed based on the smooth Dirac delta function, which is also called delta function interpolation method (DFIM). The performance

The micro-channel heat dissipation system has minor specifications and good thermal conductivity per unit, which is the best choice for heat dissipation of micro-chips. By optimizing the cross section of microchannel, the heat exchange efficiency and temperature uniformity can be effectively improved. In this article, a double-layer triangular microchannel heat sink is proposed, which uniquely combines

(2020). Correction. Numerical Heat Transfer, Part B: Fundamentals: Vol. 77, No. 5, pp. 429-429.

Extensive computational investigations have been performed to obtain more detailed information about the peculiar phenomena of turbulent supercritical carbon dioxide (sCO2) flow as ideal heat transfer fluids in various thermal engineering applications. This paper reviews the simulation techniques used and discusses their advantages, shortcomings and applicability. Not only is a comprehensive inspection