Search results

The purpose of this paper is to compute flow effects of the transition from adherence-to-slip in two-dimensional flows, for a polymer melt obeying a memory-integral viscoelastic equation, in isothermal and non-isothermal cases.

Temperature dependence is expressed by Arrhenius and William-Landel-Ferry models. A coupling approach is defined. For the dynamic equations, the Stream-Tube Method (STM) is used with finite differences in a mapped rectangular domain of the real domain, where streamlines are parallel and straight. STM avoids particle-tracking problems and allows simple formulae to evaluate stresses resulting from the constitutive equation. For the temperature field, a finite-element method is carried out to solve the energy equation in the real domain.

The approach avoids numerical problems arising with classical formulations and proves to be robust and efficient. Large elasticity levels are attained without convergence and refinement difficulties that may arise close to the “stick-slip” transition section. The method highlights the role of temperature conditions and reveals interesting differences for the ducts considered.

The results of the study are of interest for polymer processing where slip at the wall can be encountered, in relation with the physical properties of the materials.

The paper presents a simple approach that limits considerably numerical problems coming from stick-slip boundary conditions and avoids particle-tracking. Results are obtained at flow rates encountered in industrial conditions.

This paper aims to develop a hybrid finite difference‐finite element method and apply it to solve the three‐dimensional energy equation in non‐isothermal fluid flow past over a tube.

To implement the hybrid scheme, the tube length is partitioned into uniform segments by choosing grid points along its length, and a plane perpendicular to the tube axis is drawn at each of the points. Subsequently, the Taylor‐Galerkin finite element technique is employed to discretize the energy equation in the planes; while the derivatives along the tube are discretized using the finite difference method.

To demonstrate the validity of the proposed numerical scheme, three‐dimensional test cases have been solved using the method. The variation of L²‐norm of the error with mesh refinement shows that the numerical solution converges to the exact solution with mesh refinement. Moreover, comparison of the computational time duration shows that the proposed method is approximately three times faster than the 3D finite element method. In the non‐isothermal fluid flow around a tube for Re=250 and Pr=0.7, the results show that the Nusselt number decreases with the increase in the tube length and, for the tube length greater than six times the tube diameter, the average Nusselt number converges to the value for the two‐dimensional case.

A hybrid finite difference‐finite element method has been developed and applied to solve the 3D transient energy equation for different test cases. The proposed method is faster, and computationally more efficient, compared with the 3D finite element method.

The purpose of this work is to study the effect of surface-radiation on the phenomenon of natural convection flow of a Newtonian fluid in a non-Darcian porous media cavity. The study is mainly focused on the interaction between the inertial resistance of the fluid layers and the surface radiation.

For numerical simulation of transient vorticity transport and energy equations, the paper uses the alternate direct implicit method. Forward Time Central Space descretization is used for the transient and diffusion terms in the alternate direct implicit method, whereas for the convective terms, the method is modified using second upwind differencing technique. ADI method is adopted here, since this technique is unconditionally stable as a complete sweep and is second-order accurate in time for low velocity changes. The stream function equation is solved using the successive over relaxation technique with residual tolerance of order 10-5.

It was found that despite the reduction of flow, the heat transfer increases as the Forschheimer resistance is increased. Further, with the increase in the Planck number, the heat transfer from the bottom radiating wall increases. Darcy drag parameter did not have a significant impact on flow properties except a slight reduction in the flow. Nevertheless, the increase in temperature ratio has a significant impact on flow properties.

The analysis is valid for unsteady, two-dimensional natural convection flow in a fluid-saturated non-Darcy porous medium enclosed by non-isothermal walls. As a first case, the study is conducted for square cavity. An extension to three-dimensional flow case and the study of Darcy-Forschheimer medium with effect of viscous dissipation is left as a part of future work.

The approach is applicable to the modeling of geothermal systems where the inertial resistance to flow also comes into act with the non-uniform temperature distribution. The method is very useful to analyze solar receiver systems, fire research, electronic cooling, brake housing of an aircraft and many environmental geothermal processes.

The study may be of some interest to engineers interested in heat transfer in ventilated rooms or enclosures, the industrial waste, water and atmospheric pollution.

This paper aims to focus on linear and non-linear convection in a lid-driven square cavity with isothermal and non-isothermal bottom surface.

It is assumed that the top moving wall is adiabatic and the bottom wall is heated in two modes, and the rest of the walls are maintained at uniform cold temperature. The coupled governing non-linear partial differential equations are solved numerically with MAC algorithm for conducting a parametric study with uniform and non-uniform temperature bottom wall.

The numerical results are depicted in the form of streamlines, temperature contours and variation of local Nusselt number. The local Nusselt number at the bottom wall of the cavity increases in presence of non-linear temperature parameter as compared with linear temperature parameter and heat transfer reduces with increasing of Ha for uniform and non-uniform heating of bottom wall.

The numerical investigation is conducted for unsteady, two-dimensional natural convective flow in a square cavity. An extension of the present study with the effect of inclination of cavity, wavy walls and triangular cavity will be the interest of future work.

This work studies the effect of magnetic field in the presence of linear convection and non-linear convection. This study might be useful to cooling of electronic components, alloy casting, crystal growth and fusion reactors, etc.

Modelling of fiber suspension in injection molding cavities is very complex, with fluid flow, fiber orientation, and heat transfer effects taking place at the same time. Moreover, the flow is modified by the presence of fibers and vice versa. Therefore, the aim of the paper is to develop a Computational Fluid Dynamics (CFD) model to simulate and characterise the fiber suspension flow in two dimensional mold cavities. The model is intended to describe the fluid flow and heat transfer aspects of the suspension, and to predict the fiber orientation. The Navier-Stokes equations and the Jeffery (1922) equation are the governing equations for the velocity field and fiber motion respectively. The flow is considered to be two-dimensional incompressible, non-isothermal, transient and behave as non-Newtonian fluid containing suspension of short-fibers. The Finite Volume Method (FVM) combined with Control Volume Method is used to simulate the flow field by solving the momentum, energy and fiber orientation equations. To validate the numerical model, the numerical results are compared with available experimental findings. A good agreement between the numerical results and the experimental data is achieved. Since the behaviour of fiber suspension has great significance on the quality of the final product, this study has wide background of engineering application.

The purpose of this paper is to present an upscale theory of the thermal-mechanical coupling particle simulation for non-isothermal problems in two-dimensional quasi-static system, under which a small length-scale particle model can exactly reproduce the same mechanical and thermal results with that of a large length-scale one.

The objective is achieved by extending the upscale theory of particle simulation for two-dimensional quasi-static problems from an isothermal system to a non-isothermal one.

Five similarity criteria, namely geometric, material (mechanical and thermal) properties, gravity acceleration, (mechanical and thermal) time steps, thermal initial and boundary conditions (Dirichlet/Neumann boundary conditions), under which a small-length-scale particle model can exactly reproduce both the mechanical and thermal behavior with that of a large length-scale model for non-isothermal problems in a two-dimensional quasi-static system are proposed. Furthermore, to test the proposed upscale theory, two typical examples subjected to different thermal boundary conditions are simulated using two particle models of different length scale.

The paper provides some important theoretical guidances to modeling thermal-mechanical coupled problems at both the engineering length scale (i.e. the meter scale) and the geological length scale (i.e. the kilometer scale) using the particle simulation method directly. The related simulation results from two typical examples of significantly different length scales (i.e. a meter scale and a kilometer scale) have demonstrated the usefulness and correctness of the proposed upscale theory for simulating non-isothermal problems in two-dimensional quasi-static system.

The purpose of this paper is to examine the macroscopic and microscopic fields of fiber suspensions in the non‐isothermal situations, also to examine the effect of fiber on this non‐isothermal system.

Control equations are coupled and simultaneously solved by collocated finite volume method on fully triangular meshes.

Temperature dependence and wall temperature have significant effect on both macroscopic and microscopic fields of fiber suspensions. Moreover, the influence of fiber on the non‐isothermal system is similar to that of the isothermal system.

This is the first time that the microstructures of both molecules and fibers are presented in the non‐isothermal condition and it is hoped that the results will provide more insight into the microscopics of complex flows.

In this article we report on progress in the development of software tools for fluid flow prediction in the polymer processing industry. This involves state‐of‐the‐art numerical techniques and the study of a number of non‐trivial model flow problems, in an effort to investigate realistic transient problems relevant to industrial processes. Here we study particularly the effects of variations in non‐Newtonian and heat transfer properties of the flowing materials in the flows, both throughout the transient development period and at steady‐state.

This paper gives a review of the finite element techniques (FE) applied in the area of material processing. The latest trends in metal forming, non‐metal forming, powder metallurgy and composite material processing are briefly discussed. The range of applications of finite elements on these subjects is extremely wide and cannot be presented in a single paper; therefore the aim of the paper is to give FE researchers/users only an encyclopaedic view of the different possibilities that exist today in the various fields mentioned above. An appendix included at the end of the paper presents a bibliography on finite element applications in material processing for 1994‐1996, where 1,370 references are listed. This bibliography is an updating of the paper written by Brannberg and Mackerle which has been published in Engineering Computations, Vol. 11 No. 5, 1994, pp. 413‐55.

The non‐isothermal filling of a powder/binder mixture in metal injection moulding is simulated by the viscoplastic and the heat conduction finite element methods. Proposes a simplified three‐dimensional scheme for the moulding of products with a non‐uniform thickness distribution. The computing time for the simplified three‐dimensional scheme is of the same order as that for two‐dimensional problems. Deals with complex overlapping between the surfaces of the mixture, resulting from the occurrence of jetting during the moulding, by the use of a remeshing scheme. The material flow in metal injection moulding into a rectangular die with a linear thickness distribution is simulated. The jetting behaviour is remarkably influenced by the thickness distribution of the die.