Search results

The purpose of this paper is to address unsteady heat conduction in two subsets of ordinary bodies. One subset consists of a large plane wall, a long cylinder and a sphere in one dimension. The other subset consists of a short cylinder and a large rectangular bar in two dimensions. The prevalent assumptions in the two subsets are: constant initial temperature, uniform surface heat flux and thermo-physical properties invariant with temperature. The engineering applications of the unsteady heat conduction deal with the determination of temperature–time histories in the two subsets using electric resistance heating, radiative heating and fire pool heating.

To this end, a novel numerical procedure named the enhanced method of discretization in time (EMDT) transforms the linear one-dimensional unsteady, heat conduction equations with non-homogeneous boundary conditions into equivalent nonlinear “quasi–steady” heat conduction equations having the time variable embedded as a time parameter. The equivalent nonlinear “quasi–steady” heat conduction equations are solved with a finite difference method.

Based on the numerical computations, it is demonstrated that the approximate temperature–time histories in the simple subset of ordinary bodies (large plane wall, long cylinder and sphere) exhibit a perfect matching over the entire time domain 0 < t < ∞ when compared against the rigorous exact temperature–time histories expressed by classical infinite series. Furthermore, using the method of superposition of solutions in the convoluted subset (short cylinder and large rectangular crossbar), the same level of agreement in the approximate temperature–time histories in the simple subset of ordinary bodies is evident.

The performance of the proposed EMDT coupled with a finite difference method is exhaustively assessed in the solution of the unsteady, one-dimensional heat conduction equations with prescribed surface heat flux for: a subset of one-dimensional bodies (plane wall, long cylinder and spheres) and a subset of two-dimensional bodies (short cylinder and large rectangular bar).

The purpose of this study is to address one-dimensional, unsteady heat conduction in a large plane wall exchanging heat convection with a nearby fluid under “small time” conditions.

The Transversal Method of Lines (TMOL) was used to reformulate the unsteady, one-dimensional heat conduction equation in the space coordinate and time into a transformed “quasi-steady”, one-dimensional heat conduction equation in the space coordinate housing the time as an embedded parameter. The resulting ordinary differential equation of second order with heat convection boundary conditions is solved analytically with the method of undetermined coefficients.

Semi-analytical TMOL dimensionless temperature profiles of compact form with/without regressed terms are obtained for the whole spectrum of Biot number (0 < Bi < ∞) in the “small time” sub-domain. In addition, a new “large time” sub-domain is redefined, that is, setting a smaller critical dimensionless time or critical Fourier number τ_cr = 0.18.

The computed dimensionless center, surface and mean temperature profiles in the large plane wall accounting for all Biot number (0 < Bi < ∞) in the “small time” sub-domain τ < τ_cr = 0.18 exhibit excellent quality while carrying reasonable relative errors for engineering applications. The exemplary level of accuracy indicates that the traditional evaluation of the center, surface and mean temperatures with the standard infinite series retaining a large number of terms is no longer necessary.

The purpose of this paper is to present outcomes of efforts made over the last 20 years to extend the applicability of the Richardson extrapolation procedure to numerical predictions of multidimensional, steady and unsteady, fluid flow and heat transfer phenomena in regular and irregular calculation domains.

Pattern-preserving grid-refinement strategies are proposed for mathematically rigorous generalizations of the Richardson extrapolation procedure for numerical predictions of steady fluid flow and heat transfer, using finite volume methods and structured multidimensional Cartesian grids; and control-volume finite element methods and unstructured two-dimensional planar grids, consisting of three-node triangular elements. Mathematically sound extrapolation procedures are also proposed for numerical solutions of unsteady and boundary-layer-type problems. The applicability of such procedures to numerical solutions of problems with curved boundaries and internal interfaces, and also those based on unstructured grids of general quadrilateral, tetrahedral, or hexahedral elements, is discussed.

Applications to three demonstration problems, with discretizations in the asymptotic regime, showed the following: the apparent orders of accuracy were the same as those of the numerical methods used; and the extrapolated results, measures of error, and a grid convergence index, could be obtained in a smooth and non-oscillatory manner.

Strict or approximate pattern-preserving grid-refinement strategies are used to propose generalized Richardson extrapolation procedures for estimating grid-independent numerical solutions. Such extrapolation procedures play an indispensable role in the verification and validation techniques that are employed to assess the accuracy of numerical predictions which are used for designing, optimizing, virtual prototyping, and certification of thermofluid systems.

The purpose of this paper is to describe an extension of the boundary element method (BEM) and the dual reciprocity boundary element method (DRBEM) formulations developed for one- and two-dimensional steady-state problems, to analyse transient convection–diffusion problems associated with first-order chemical reaction.

The mathematical modelling has used a dual reciprocity approximation to transform the domain integrals arising in the transient equation into equivalent boundary integrals. The integral representation formula for the corresponding problem is obtained from the Green’s second identity, using the fundamental solution of the corresponding steady-state equation with constant coefficients. The finite difference method is used to simulate the time evolution procedure for solving the resulting system of equations. Three different radial basis functions have been successfully implemented to increase the accuracy of the solution and improving the rate of convergence.

The numerical results obtained demonstrate the excellent agreement with the analytical solutions to establish the validity of the proposed approach and to confirm its efficiency.

Finally, the proposed BEM and DRBEM numerical solutions have not displayed any artificial diffusion, oscillatory behaviour or damping of the wave front, as appears in other different numerical methods.

The paper aims to report a fire test conducted on a three-dimensional frame in order to investigate the behaviour of bare steel flush end-plate connections with relatively low thickness at elevated temperatures.

A half-scale model was fabricated and exposed to modified (scaled) ISO 834 heating curve using a semi-open furnace. The maximum temperature inside the furnace reached 1,026 °C.

The rotations of connections are reported and compared with those of a previous study on an exactly the same model with thick end-plates. Various modes of failure such as local buckling of the beams flanges and lateral-torsional buckling of beams were observed during the test. Finally, the structure collapsed after 29 min of heating due to the fracture of weld between one of the beams and one of its attached end-plates whilst the other beam had undergone a maximum deflection of 35 cm (≈ 1/6 span length). Other observed failure modes included bolt fracture, bolt thread stripping and large inelastic deformation of the end-plates.

Although the adoption of thin end-plates increased the rotational capacity of the connections, it did not improve the robustness of the structure under fire conditions.

The purpose of this paper is to present a new approach of meshless local B-spline based finite difference (FD) method for solving two dimensional transient heat conduction problems.

In the present method, any governing equations are discretized by B-spline approximation which is implemented in the spirit of FD technique using a local B-spline collocation scheme. The key aspect of the method is that any derivative is stated as neighbouring nodal values based on B-spline interpolants. The set of neighbouring nodes are allowed to be randomly distributed thus enhanced flexibility in the numerical simulation can be obtained. The method requires no mesh connectivity at all for either field variable approximation or integration. Time integration is performed by using the Crank-Nicolson implicit time stepping technique.

Several heat conduction problems in complex domains which represent for extended surfaces in industrial applications are examined to demonstrate the effectiveness of the present approach. Comparison of the obtained results with solutions from other numerical method available in literature is given. Excellent agreement with reference numerical method has been found.

The method is presented for 2D problems. Nevertheless, it would be also applicable for 3D problems.

A transient two dimensional heat conduction in complex domains which represent for extended surfaces in industrial applications is presented.

The presented new meshless local method is simple and accurate, while it is also suitable for analysis in domains of arbitrary geometries.

The purpose of this paper is to examine the effects of Darcian drag force and radiation-conduction on unsteady two-dimensional magnetohydrodynamic flow of viscous, electrically conducting and Newtonian fluid over a vertical plate adjacent to a Darcian regime in presence of thermal radiation and transversal magnetic field. A well-tested, numerically stable Crank-Nicolson finite-difference procedure is employed for the conservation equations. Excellent agreement is obtained for numerical solutions with previously published work.

In this investigation, an efficient, accurate, extensively validated and unconditionally stable finite-difference scheme based on the Crank-Nicolson model is developed to solve the governing coupled, non-linear partial differential equations. The accuracy and effectiveness of the method are demonstrated.

Different numerical results are obtained and presented graphically to explain the effect of various physical parameters on the velocity and temperature profiles, local, as well as average, skin friction and Nusselt number. It is found that, with a rise in Darcian drag force, flow velocity and temperature are reduced, but increased for all times. Both average and local skin frictions are reduced considerably with an increase in Darcian drag force, but reversed behavior is observed for the local Nusselt number. Increasing the thermal radiation effects accelerated the flow velocity as well as the fluid temperature and wall local skin friction in a saturated porous medium, but effectively reduced the local Nusselt number and average Nusselt number at the wall. Comparison with previously published works in the limits shows excellent agreement.

The analysis is valid for unsteady, two-dimensional laminar flow of an optically thick no-gray gas, electrically conducting, and Newtonian fluid past an isothermal vertical surface adjacent to the Darcian regime with variable surface temperature. An extension to three-dimensional flow case is left for future work.

Practical interest of such study includes applications in electromagnetic lubrication, boundary cooling, bio-physical systems and in many branches of engineering and science. It is well known that the effect of thermal radiation is important in space technology and high temperature processes. Thermal radiation also plays an important role in controlling heat transfer process in polymer processing industry.

The purpose of this paper is to present an immersed volume method that accounts for solid conductive bodies (hat‐shaped disk) in calculation of time‐dependent, three‐dimensional, conjugate heat transfer and fluid flow.

The incompressible Navier‐Stokes equations and the heat transfer equations are discretized using a stabilized finite element method. The interface of the immersed disk is defined and rendered by the zero isovalues of a level set function. This signed distance function allows turning different thermal properties of each component into homogeneous parameters and it is coupled to a direct anisotropic mesh adaptation process enhancing the interface representation. A monolithic approach is used to solve a single set of equations for both fluid and solid with different thermal properties.

In the proposed immersion technique, only a single grid for both air and solid is considered, thus, only one equation with different thermal properties is solved. The sharp discontinuity of the material properties was captured by an anisotropic refined solid‐fluid interface. The robustness of the method to compute the flow and heat transfer with large materials properties differences is demonstrated using stabilized finite element formulations. Results are assessed by comparing the predictions with the experimental data.

The proposed method demonstrates the capability of the model to simulate an unsteady three‐dimensional heat transfer flow of natural convection, conduction and radiation in a cubic enclosure with the presence of a conduction body. A previous knowledge of the heat transfer coefficients between the disk and the fluid is no longer required. The heat exchange at the interface is solved and dealt with naturally.

The effect of radiation on unsteady natural convection in a two‐dimensional participating medium between two horizontal concentric and vertically eccentric cylinders is investigated numerically. The equations of transfer are written by using a bicylindrical coordinates system, the stream function, and the vorticity. The finite volume radiation solution method and the control volume approach are used to discretize the coupled equations of radiative transfer, momentum, and energy. Original results are obtained for three eccentricities, Rayleigh number equal to 10⁴, 10⁵, and a wide range of radiation‐conduction parameter. The effects of optical thickness, wall emissivity, and scattering on flow intensity and heat transfer are discussed.

The purpose of this study is to validate whether the local thermal equilibrium for unsteady state is an appropriate assumption for the porous media with closed pores. It also compares the transient temperatures between the pore scale and volume averaged approaches to prove that the volume averaged method is an appropriate technique for the heat transfer in closed-cell porous media. The interfacial heat transfer coefficient for the closed-cell porous media is also discussed in details.

The governing equations for the pore scale and continuum domains are given. They are solved numerically for the pore scale and volume-averaged domains. The results are compared and discussion was done. The performed discussions and explanations are supported with figure and graphics.

A local thermal non-equilibrium exits for the closed-cell porous media in which voids are filled with water during the unsteady heat transfer process. Local thermal non-equilibrium condition exists in the cells under high temperature gradient and it disappears when the heat transfer process becomes steady-state. Although a local thermal equilibrium exists in the porous media in which the voids are filled with air, a finite value for heat transfer coefficient is found. The thermal diffusivity of air and solid phase are close to each other and hence a local thermal equilibrium exists.

The study is done only for the closed-cell porous media and for Rayleigh number till 10⁵. Two common working fluids as water and air are considered.

There are many applications of porous media with closed pores particularly in the industry, such as the closed-cell metal foam or the closed cells in porous materials such as foods and plastic-based insulation material. The obtained results are important for transient heat transfer in closed-cell porous materials.

The obtained results are important from the transient application of heat transfer in the closed-cell material existing in nature and industry.

The authors’ literature survey shows that it is the first time the closed-cell porous media is discussed from local thermal non-equilibrium point of view and it is proved that the local thermal non-equilibrium can exist in the closed-cell porous media. Hence, two equations as solid and fluid equations should be used for unsteady heat transfer in a closed-cell porous medium.