Search results

This paper aims to present a novel numerical technique for solving the incompressible multiphase mixture model.

The multiphase mixture model contains a set of momentum and continuity equations for the mixture phase, a second phase continuity equation and the algebraic equation for the relative velocity. For solving continuity equation for the second phase and advection term of momentum, an improved approach fine grid advection-multiphase mixture flow (FGA-MMF) is developed. In the FGA-MMF method, the continuity equation for the second phase is solved with higher-order schemes in a two times finer grid. To solve the advection term of the momentum equation, the advection fluxes of the volume fraction in the continuity equation for the second phase are used.

This approach has been used in various tests to simulate unsteady flow problems. Comparison between numerical results and experimental data demonstrates a satisfactory performance. Numerical examples show that this approach increases the accuracy and stability of the solution and decreases non-monotonic results.

The solver for the multi-phase mixture model can only be adopted to solve the incompressible fluid flow.

The paper developed an innovative solution (FGA-MMF) to find multi-phase flow field value in the multi-phase mixture model. Advantages of the FGA-MMF technique are the ability to accurately determine the phases interpenetrating, decreasing the numerical diffusion of the interface and preventing instability and non-monotonicity in solution of large density variation problems.

In this work, a new two-phase version of the finite element-based Artificial Compressibility (AC) Characteristic-Based Split (CBS) algorithm is developed and applied for the first time to heat and mass transfer phenomena in porous media with associated phase change. The purpose of this study is to provide an alternative for the theoretical analysis and numerical simulation of multiphase transport phenomena in porous media. Traditionally, the more complex Separate Flow Model was used in which the vapour and liquid phases were considered as distinct fluids and mathematically described by the conservation laws for each phase separately, resulting in a large number of governing equations.

Even though the adopted mathematical model presents analogies with the conventional multicomponent mixture flow model, it is characterized by a considerable reduction in the number of the differential equations for the primary variables. The fixed-grid numerical formulation can be applied to the resolution of general problems that may simultaneously include a superheated vapour region, a two-phase zone and a sub-cooled liquid region in a single physical domain with irregular and moving phase interfaces in between. The local thermal non-equilibrium model is introduced to consider the heat exchange between fluid and solid within the porous matrix.

The numerical model is verified considering the transport phenomena in a homogenous and isotropic porous medium in which water is injected from one side and heated from the other side, where it leaves the computational domain in a superheated vapour state. Dominant forces are represented by capillary interactions and two-phase heat conduction. The obtained results have been compared with the numerical data available in the scientific literature.

The present algorithm provides a powerful routine tool for the numerical modelling of complex two-phase transport processes in porous media.

For the first time, the stabilized AC-CBS scheme is applied to the resolution of compressible viscous flow transport in porous materials with associated phase change. A properly stabilized matrix inversion-free procedure employs an adaptive local time step that allows acceleration of the solution process even in the presence of large source terms and low diffusion coefficients values (near the phase change point).

The purpose of this study is to experimentally and numerically scrutinize the heat transfer enhancement in pool boiling and condensation by changing the hydrophilicity or hydrophobicity properties of the working fluid, i.e. by use of nanofluid solution.

For specifying the effects of nanoparticle concentration on heat transfer properties, two different nanofluid solutions (h-BN/DCM and SiO₂/DCM) at three different volumetric concentrations were prepared and tested under different heat flux conditions. Boiling curves, alterations in pressure with heat flux and heat transfer coefficients for both boiling and condensation processes were obtained and viscosity measurements were performed for dichloromethane (DCM) and each working fluid was prepared. In addition, a series of numerical simulations, via computational fluid dynamics approach, was performed for specifying the evaporation–condensation phenomena and temperature and velocity distributions.

Nanoparticle addition inside the base fluid increased the thermal characteristics of the base fluid significantly. For the experimental results of h-BN/DCM nanofluid, the increment rate in heat transfer coefficient for saturation boiling, after-saturation boiling and condensation processes was found as 27.59%, 14.44% and 15%, respectively.

The novelty of this comparison study is that there is no such experimental and numerical comparison study in literature for DCM fluid, which concentrates on thermal performance enhancement and compares the effect of different kinds of nanoparticles on heat transfer characteristics for boiling–condensation processes.

A novel numerical formulation of the two‐phase macroscopic balance equations governing the flow field in incompressible porous media is presented. The numerical model makes use of the weighted average flux method and total variation diminishing flux limiting techniques, and results in a second‐order accurate scheme. A shock tube study was carried out to examine the interaction of a normal shock wave with a thin layer of porous, incompressible cellular ceramic foam. Particular attention was paid to the transmitted and reflected flow fields. The numerical model was used to simulate the experimental test cases, and their results compared with a view to validate the numerical model. A phenomenological model is proposed to explain the behaviour of the transmitted flow field.

This study aims to investigate numerically the problem of conjugate conduction and mixed convection heat transfer of a nanofluid in a rotational/stationary circular enclosure using a two-phase mixture model.

Hot and cold surfaces on the wall or inside the enclosure (heater and cooler) are maintained at constant temperature of T_h and T_c, respectively, whereas other parts are thermally insulated. To examine the effects of various parameters such as Richardson number (0.01 = Ri =100), thermal conductivity ratio of solid to base fluid (1 = K_r = 100), volume fraction of nanoparticle (0 = φ = 0.05), insertion of conductive covers (C.Cs) around the heater in a different shape (triangular, circular or square), segmentation and arrangement of the conductive blocks (C.Bs) and rotation direction of the enclosure on the flow structure and heat transfer rate, two-dimensional equations of mass, momentum and energy conservation, as well as volume fraction, are solved using finite volume method and Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm.

The results show that inserting C.C around heater can increase or decrease heat transfer rate, and it depends on thermal conductivity ratio of solid to pure fluid. Also, it is found that by the division of C.B and location of its portions in a horizontal configuration, heat transfer rate reduces. Moreover, it is observed that external heating and cooling of the enclosure causes enhancement of heat transfer relative to that of internal heating and cooling. Finally, results illustrate that under the condition that cylinders rotate in the same direction, the heat transfer rate increases as compared to those that rotate in the opposite direction. Hence rotation direction of cylinders can be used as a desired parameter for controlling heat transfer rate.

A comprehensive report of results for the problem of conjugate conduction and mixed convection heat transfer in a circular cylinder containing different shapes of C.C, conducting obstacle and heater and cooler has been presented. An efficient numerical technique has been developed to solve this problem. The achievements of this paper are purely original, and the numerical results were never published by any researcher.

The stringent requirements for environmental protection have induced the extensive applications of water-lubricated journal bearings in marine propulsion. The nonlinear dynamic analysis of multiple grooved water-lubricated bearings (MGWJBs) has not been fully covered so far in the literature. This study aims to conduct the nonlinear dynamic analysis of the instability for MGWJBs.

An attenuation rate interpolation method is proposed for the determination of the critical instability speed. Based on a structured mesh movement algorithm, the transient hydrodynamic force model of MGWJBs is set up. Furthermore, the parameters’ analysis of nonlinear instability for MGWJBs is conducted. The minimum water film thickness, side leakage, friction torque and power loss of friction are fully analyzed.

With the increase of speed, the journal orbits come across the steady state equilibrium motion, sub-harmonic motion and limit circle motion successively. At the limit circle motion stage, the orbits are much larger than that of steady state equilibrium and sub-harmonic motion. The critical instability speed increases when the spiral angle decreases or the groove angle increases. The minimum water film thickness peak is at the rotor speed of 4,000 r/min for the MGWJB with Sa = 0°. As rotor speed increases, the side leakage decreases slightly while the friction torque and the power loss of friction increase gradually.

Present research provides a beneficial reference for the dynamic mechanism analysis and design of MGWJBs.

Corrosion in deadlegs occurs as a result of water separation due to the very low flow velocity. The present work aims to investigate the effect of geometry on flow field oil/water separation in deadlegs in an attempt for obtaining the conditions for avoiding formation of deadleg.

The investigation is based on the solution of the mass and momentum conservation equations of an oil/water mixture together with the volume fraction equation for the secondary phase. A fluid flow model based on the time‐averaged governing equation of 3D turbulent flow has been developed. An algebraic slip mixture model for the calculation of the two immiscible fluids (water and crude oil) is utilized.

Results are obtained for different lengths of the deadleg. The inlet flow velocity is kept unchanged (1.0 m/s) and the deadleg length to diamter ratio (L/D_B) ranges from 1 to 7. The considered fluid mixture contains 90 percent oil and 10 percent water (by volume). The results show that the size of the stagnant fluid region increases with the increase of L/D_B 1≈3D_B.

Deadlegs should be avoided whenever possible in design of piping for fluids containing or likely to contain corrosive substance. When deadlegs are unavoidable, the length of the inactive pipe must be as short as possible to avoid stagnant or low‐velocity flows.

The model solves the continuity and momentum equations for the mixture, and the volume fraction equation for the secondary phase utilizing an algebraic expression for the relative velocity.

This is a 3D numerical study of convective heat transfer through a micro concentric annulus governing non-uniform heat flux boundary conditions employing water-Al₂O₃ nanofluid. The nanofluid is modeled using two-phase mixture model, as it has a good agreement to experimental results.

Half of the inner pipe surface area of the annulus section of a double pipe heat exchanger is exposed to a constant heat flux which two models are considered to divide the exposing surface area to smaller ones considering the fact that in all cases half of the inner pipe surface area has to be exposed to the heat flux: in model (A), the exposing surface area is divided radially to two parts (A1), four parts (A2) and eight parts (A3) by covering the whole length of the annulus and in model (B) the exposing surface area is divided axially to two parts (B1), four parts (B2) and eight parts (B3) by covering half of the annulus radially.

The results reveal that model (B) leads to higher Nusselt numbers compared to model (A); however, at Reynolds number 10, model (A3) exceeds model (B3). The average Nusselt number is increased up to 142 and 83 per cent at models (A3) with Reynolds number 10 and model (B3) with Reynolds number 1000, respectively.

This paper is a two-phase investigation of water-Al₂O₃ nanofluid in a micro concentric annulus under non-uniform heat flux boundary conditions.

The purpose of this study is three-dimensional flow and heat transfer investigation of water/Al₂O₃ nanofluid inside a microchannel with different cross-sections in two-phase mode.

The effect of microchannel walls geometry (trapezoidal, sinusoidal and stepped microchannels) on flow characteristics and also changing circular cross section to trapezoidal cross section in laminar flow at Reynolds numbers of 50, 100, 300 and 600 were investigated. In this study, two-phase water/Al₂O₃ nanofluid is simulated by the mixture model, and the effect of volume fraction of nanoparticles on performance evaluation criterion (PEC) is studied. The accuracy of obtained results was compared with the experimental and numerical results of other similar papers.

Results show that in flow at lower Reynolds numbers, sinusoidal walls create a pressure drop in pure water flow which improves heat transfer to obtain PEC < 1. However, in sinusoidal and stepped microchannel with higher Reynolds numbers, PEC > 1. Results showed that the stepped microchannel had higher pressure drop, better thermal performance and higher PEC than other microchannels.

Review of previous studies showed that existing papers have not compared and investigated nanofluid in a two-phase mode in inhomogeneous circular, stepped and sinusoidal cross and trapezoidal cross-sections by considering the effect of changing channel shape, which is the aim of the present paper.

Flow force is one of the crucial factors affecting the performance of conical throttle valves. The purpose of this paper is to determine the relationship between the flow force and operating parameters of the conical throttle valve.

The flow force of the throttle valve can be obtained by the difference between the axial force and static pressure on the valve spool. In this paper, the internal fluid is divided into two regions and the axial force and static pressure are obtained, respectively. In addition, a two-dimensional axisymmetric simulation model and experimental test are carried out to validate the results of the flow force.

It can be seen that the theoretical, simulation and experimental results exhibit high agreement with each other and the error between them decreases with the increase in the size of the opening. The curves of pressure distribution reveal that the pressure on the spool first decreases then increases when it reaches the minimum pressure at the orifice. Additionally, the minimum pressure decreases with the increase of opening and pressure difference. The results also indicate that the increase in the size of the opening and inlet pressure has a positive effect on the flow force. However, the increase in outlet pressure has a negative effect on the flow force.

In this paper, the flow force calculation model of conical throttle valve is established and the influence of operating parameters on the flow force of conical throttle valve is studied.