Maximum density effects on convective instability of horizontal plane poiseuille flows in the thermal entrance region

A linear stability analysis is used to study the conditions marking the onset of secondary flow in the form of longitudinal vortices for plane Poiseuille flow of water in the thermal entrance region of a horizontal parallel-plate channel by a numerical method. The water temperature range under consideration is 0∼30°C and the maximum density effect at 4°C is of primary interest. The basic flow solution for temperature includes axial heat conduction effect and the entrance temperature is taken to be uniform at far upstream location jackie=−∞ to allow for the upstream heat penetration through thermal entrance jackie=0. Numerical results for critical Rayleigh number are obtained for Peclet numbers 1, 10, 50 and thermal condition parameters (λ ₁, λ ₂) in the range of −2.0≤λ ₁≤−0.5 and −1.0≤λ ₂≤1.4. The analysis is motivated by a desire to determine the free convection effect on freezing or thawing in channel flow of water.     

The thermoconvective instability of plane poiseuille flow heated from below: A proposed benchmark solution for open boundary flows

The transient two-dimensional Navier-Stokes and energy equations have been solved numerically for flow in a horizontal channel heated from below in the Boussinesq limit. For the set of dimensionless parameters chosen, the flow consists of periodic transverse travelling waves resulting from a convective instability. The solution is proposed as a benchmark for the application of outflow boundary conditions (OBC) in time-dependent flows with strong buoyancy effects. Richardson extrapolation in both time and space is used in obtaining the solution. Field plots and profiles of velocity, temperature, vorticity and streamfunction at selected axial positions and times are also presented from the finest grid and smallest time step calculation. The calculations have been made on an extended domain so that the effects of OBC used in the present study would be negligible in the test region.     

Analytical and numerical solution of the poiseuille flow of a Johnson-Segalman fluid

The axisymmetric Poiseuille flow of a Johnson-Segalman fluid is selected as a specific one-dimensional example for testing numerical algorithms, while analytical results are also made available. A particular emphasis is put on the appearance of limit points with some numerical algorithms.     

On pressure and thermal flux in gas-mixture flows and in two-phase flows

To begin with, two different definitions of pressure, thermal flux, etc. in the diffusion model and two-fluid model are given. Then the physical interpretations of the pressure and the thermal flux are provided by introducing the momentum and energy fluxes,M and ε, through a surface dS in the flow field. The quantities defined in the diffusion model are suggested when the motion of the mixture is studied as a whole, while the quantities defined in the two-fluid model are suggested when the motion of individual species is studied. The collision pressure and thermal flux in dense gas-mixtures are also discussed in detail, i.e. their origin, their expressions in the momentum and energy equations, and their distinctions from the normal partial pressure and thermal flux. A gas-particle flow can be treated as a flow of dense gas-mixtures. The long-standing controversy whether the “inertial coupling term” should exist in the momentum equation can be clarified by the two different definitions of pressure.     

Shock slip-relations for thermal and chemical nonequilibrium flows

This paper appears to be the first where the multi-temperature shock slip-relations for the thermal and chemical nonequilibrium flows are derived. The derivation is based on analysis of the influences of thermal nonequilibrium and viscous effects on the mass, momentum and emergy flux balance relations at the shock wave. When the relaxation times for all internal energy modes tend to zero, the multi-tmperature shock slip-relations are converted into single-temperature ones for thermal equilibrium flows. The present results can be applied to flow over vehicles of different geometries with or without angles of attack. In addition, the present single-temperature shock slip-relations are compared with those in the literature, and some defects and limitations in the latter are clarified. The project supported by the National Natural Science Foundation of China and the National Defence Science and Industry Commission of China.     

Numerical solution of subsonic and transonic cascade flows

In this work a study of the application of the finite element method to transonic flows in axial turbomachines is undertaken. Solution techniques capable of accurately predicting flows from the incompressible regime up to the establishment of shocks in the transonic regime are presented. In the subsonic and shockless transonic regimes a local linearization method capable of very rapid convergence is used. In the full transonic regime the artificial compressibility method is employed to exclude downstream influences in the supersonic regions. The two approaches can be combined in a unified package and appropriate switches introduced to select the relevant method in any flow regime.     

Influence of thermal aging on tensile and creep behavior of thermoplastic polyurethane

Changes in mechanical and physical properties of polyurethane thermoplastic during aging at 70 °C and 90 °C were investigated. The loss weight response was analyzed by gravimetric measurements under these temperatures. Changes in appearance and morphology of TPU after thermal aging were revealed by optical microscopy. The prolongation of the thermal exposure time, up to 270 days, leads to a progressive increase in tensile strength. In fact, elastic modulus and stress at 200% of strain were increased with thermal exposure time. These results can be explained by the increase of thermal stability due to the increase of material rigidity and the decrease in chain mobility. The evolution of the mechanical properties from tensile tests seems to be well correlated to the creep behavior. Finally, Scanning Electron Microscopy (SEM) revealed the modification of TPU morphology fracture surface after thermal aging.     

Multi-dimensional hydrodynamic code SOVA for interfacial flows: Application to the thermal layer effect

Eulerian, three-dimensional, numerical code, which conserves mass, momentum and energy simultaneously both in the Lagrangian and remap steps, is developed. The use of special form of linear viscosity provides a weaker time step restriction as compared with the Courant condition. The code is designed to investigate the multi-material problems, including dusty flows. The performance of the code is illustrated by the modeling of shock wave interaction with a dusty thermal layer. Received 15 June 1998/ Accepted 6 January 1999     

Massively parallel mesh adaptation and linear system solution for multiphase flows

ABSTRACT

In this paper, a work performed to allow massively parallel finite element flow computations is presented. It includes the development and optimisation of two particular features of a finite element multiphase computational fluid dynamics software, which are mesh generation and linear system solution, using anisotropic adaptation and multigrid preconditioning. Parallel performances on supercomputers are shown, where the largest generated mesh (on 65 536 Intel Xeon or 261 144 Power PC cores) had 33.4 billions of nodes, leading to a 100 billion of unknowns linear system solution. Final applications concern, between others, image-based flow simulations.     

Comparison of different shear rheometers with regard to creep and creep recovery measurements

The performance of two shear rheometers with regard to creep and creep recovery measurements is investigated. The first one is a commercially available stress-controlled rheometer that uses an air bearing, and the second one is a magnetic bearing torsional creep apparatus that was built at the Institute of Polymer Materials in Erlangen. The creep and creep recovery measurements were performed in the linear-viscoelastic regime of two polyethylene melts at a temperature of 150°C. The creep compliances of the polyethylenes measured by both rheometers are in excellent agreement. The recoverable compliances of both polyethylene melts, however, have lower values in the case of the commercial rheometer than in the case of the magnetic bearing rheometer. The experimental parameters of creep recovery experiments and the features of the two bearings that are responsible for the different results of the two rheometers are discussed. The influence of the level of the applied shear stress on the short time behavior is investigated as well as the residual torque of the bearings which influences the long time region of the recoverable compliances. Also addressed are the influence of the momentum of inertia of the rotor and the bearing friction which is different for the two rheometers.     

Numerical solution scheme for inert,disperse, and dilute gas-particle flows

A multi-velocity formulation is proposed for the solution of an Eulerian representation of an inert, disperse, and dilute particle-phase of a gas-particle flow. Single-velocity formulations are capable of predicting regions of zero particle concentration but are problematic with crossing particle trajectories or compression waves. The multi-velocity formulation described here can account for crossing particle trajectories by splitting the particle-phase into distinct velocity families which are transported separately in the flow. Switching of the particle families at solid boundaries and due to momentum transfer with the gas-phase is conducted in a manner that enforces conservation of mass, momentum, and energy. This numerical method is combined with a parallel block-based adaptive mesh refinement algorithm that is very effective in treating problems with disparate length scales. The block-based data structure lends itself naturally to domain decomposition and thereby enables efficient and scalable implementations of the algorithm on distributed-memory multi-processor architectures. Numerical results are described to demonstrate the capabilities of the approach for predicting gas-particle flows.     

Modeling the turbulent heat and momentum transfer in flows under different thermal conditions

Two-equation turbulence models for velocity and temperature (scalar) fields are developed to calculate wall shear flows under various flow conditions and related turbulent heat transfer under various wall thermal conditions. In the present models, we make the modified dissipation rates of both turbulent energy and temperature variance zero at a wall, though the wall limiting behavior of velocity and temperature fluctuations is reproduced exactly. Thus, the models assure computational expediency and convergence. Also, the present k- model is construted using a new type of expression for the Reynolds stress proposed by Abe et al. [Trans. JSME B 61 (1995) 1714–1721], whose essential feature lies in introducing the explicit algebraic stress model concept into the nonlinear k- formulation, and the present two-equation heat transfer model is constructed to properly take into account the effects of wall thermal conditions on the eddy diffusivity for heat. The models are tested with five typical velocity fields and four typical thermal fields. Agreement with experiment and direct simulation data is quite satisfactory.     

Stabilized finite element solution to handle complex heat and fluid flows in industrial furnaces using the immersed volume method

We consider the numerical simulation of conjugate heat transfer, incompressible turbulent flows for multicomponents systems using a stabilized finite element method. We present an immersed volume approach for thermal coupling between fluids and solids for heating high‐alloy steel inside industrial furnaces. It consists in considering a single 3D grid of the furnace and solving one set of equations with different thermal properties. A distance function enables to define precisely the position and the interface of any objects inside the volume and to provide homogeneous physical and thermodynamic properties for each subdomain. An anisotropic mesh adaptation algorithm based on the variations of the distance function is then applied to ensure an accurate capture of the discontinuities that characterize the highly heterogeneous domain. The proposed method demonstrates the capability of the model to simulate an unsteady three‐dimensional heat transfers and turbulent flows in an industrial furnace with the presence of three conducting solid bodies. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element solution of three‐dimensional turbulent flows applied to mold‐filling problems

This paper presents a finite element solution algorithm for three‐dimensional isothermal turbulent flows for mold‐filling applications. The problems of interest present unusual challenges for both the physical modelling and the solution algorithm. High‐Reynolds number transient turbulent flows with free surfaces have to be computed on complex three‐dimensional geometries. In this work, a segregated algorithm is used to solve the Navier–Stokes, turbulence and front‐tracking equations. The streamline–upwind/Petrov–Galerkin method is used to obtain stable solutions to convection‐dominated problems. Turbulence is modelled using either a one‐equation turbulence model or the κ–ε two‐equation model with wall functions. Turbulence equations are solved for the natural logarithm of the turbulence variables. The change of dependent variables allows for a robust solution algorithm and good predictions even on coarse meshes. This is very important in the case of large three‐dimensional applications for which highly refined meshes result in untreatable large numbers of elements. The position of the flow front in the mold cavity is computed using a level set approach. Finally, equations are integrated in time using an implicit Euler scheme. The methodology presents the robustness and cost effectiveness needed to tackle complex industrial applications. Copyright © 2000 John Wiley & Sons, Ltd.     

Comparing analytical and numerical solution of nonlinear two and three‐dimensional hydrostatic flows

New test cases for frictionless, three‐dimensional hydrostatic flows have been derived from some known analytical solutions of the two‐dimensional shallow water equations. The flow domain is a paraboloid of revolution and the flow is determined by the initial conditions, the nonlinear advective terms, the Coriolis acceleration and by the hydrostatic pressure. Wetting and drying is also included. Some specific properties of the exact solutions are discussed under different hypothesis and relative importance of the forcing terms. These solutions are proposed for testing the stability, the accuracy and the efficiency of numerical models to be used for simulating environmental hydrostatic flows. The computed solutions obtained with a semi‐implicit finite difference—finite volume algorithm on unstructured grid are compared with the corresponding analytical solutions in both two and three space dimension. Excellent agreement are obtained for the velocity and for the resulting water surface elevation. Comparison of the computed inundation area also shows a good agreement with the analytical solution with degrading accuracy observed when the inundation area becomes relatively large and for long simulation time. Copyright © 2006 John Wiley & Sons, Ltd.     

On the numerical solution of the turbulence energy equations for wave and tidal flows

This paper deals with the numerical solution, using finite difference methods, of the hydrodynamic and turbulence energy equations which describe wind wave and tidally induced flow. Calculations are performed using staggered and non-staggered finite difference grids in the vertical, with various time discretizations of the production and dissipation terms in the turbulence energy equations. It is shown that the time discretization of these terms can significantly influence the stability of the solution. The effect of time filtering on the numerical stability of the solution is also considered. The form of the mixing length is shown to significantly influence the bed stress in wind wave problems. A no-slip condition is applied at the sea bed, and the associated high-shear bottom boundary layer is resolved by transforming the equations onto a logarithmic or log-linear co-ordinate system before applying the finite difference scheme. A computationally economic method is developed which remains stable even when a very fine vertical grid (over 200 points) is used with a time step of up to 30 min.     

A segregated implicit solution algorithm for 2D and 3D laminar incompressible flows

A segregated algorithm for the solution of laminar incompressible, two- and three-dimensional flow problems is presented. This algorithm employs the successive solution of the momentum and continuity equations by means of a decoupled implicit solution method. The inversion of the coefficient matrix which is common for all momentum equations is carried out through an approximate factorization in upper and lower triangular matrices. The divergence-free velocity constraint is satisfied by formulating and solving a pressure correction equation. For the latter a combined application of a preconditioning technique and a Krylov subspace method is employed and proved more effecient than the approximate factorization method. The method exhibits a monotonic convergence, it is not costly in CPU time per iteration and provides accurate solutions which are independent of the underrelaxation parameter used in the momentum equations. Results are presented in two- and three-dimensional flow problems.     

Poiseuille flow and thermal creep in a capillary tube on the basis of the kinetic R-model

Flow of a diatomic rarefied gas in a capillary tube of infinite length and an arbitrary cross-section under a given small pressure gradient (Poiseuille flow) or a small temperature gradient (thermal creep) is studied on the basis of a kinetic model that takes account for the rotational degrees of freedom of molecules (R-model). Numerical investigation is carried out for flows between parallel flat plates and in a circular capillary tube at the gas parameters corresponding to nitrogen. The main calculated quantity is the gas flow rate through a tube cross-section. The results are compared with the corresponding data obtained on the basis of the S-model.     

Polymer contribution to the thermal conductivity and viscosity in a dilute solution (Fraenkel Dumbbell model)

The phase-space kinetic theory for polymeric liquid mixtures is used to obtain an expression for the polymer contribution to the thermal conductivity of a nonflowing, dilute solution of polymers, where the polymer molecules are modeled as Fraenkel dumbbells. This theory takes into account three mechanisms for the energy transport: diffusion of kinetic energy (including the Öttinger-Petrillo term), diffusion of intramolecular energy, and the work done against the intramolecular forces. This paper is an extension of previous developments for the Hookean dumbbell model and the finitely-extensible dumbbell model. A comparison among the dumbbell results suggests that the thermal conductivity increases with chain stiffness. In addition, the zero-shear-rate viscosity and first normal-stress coefficient are also given for the Fraenkel dumbbell model.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.