Application of Lattice Boltzmann Methods in Heat and Moisture Transfer in Frozen Soil

Although water only takes a little percentage in the total mass of soil, it indeed plays an important role to the strength of structure. Moisture transfer can be carried out by many different mechanisms which may involve heat and mass transfer, thermodynamic phase change, and the interplay of various forces such as viscous, buoyancy, and capillary forces. The continuum models are not well suited for describing those phenomena in which the connectivity of the pore space or the fracture network, or that of a fluid phase, plays a major role. However, Lattice Boltzmann methods (LBMs) are especially well suited to simulate flows around complex geometries. Lattice Boltzmann methods were initially invented for solving fluid flows. Recently, fluid with multicomponent and phase change is also included in the equations. By comparing the numerical result with experimental result, the Lattice Boltzmann methods with phase change will be optimized.

Heat Transfer Analysis of Rectangular Channel Plate Heat Sink

In order to improve the simulation effects of space cold black environment, this paper described a rectangular channel plate heat sink. By using fluid mechanics theory and finite element method, the internal fluid flow and heat transfer in heat sink was numerically simulated to analyze the impact of channel structural on fluid flow and heat transfer. The result showed that heat sink temperature uniformity is well, and the impact of channel structural on the heat sink temperature uniformity is not significant. The channel depth and spacing are important factors which affect the fluid flow and heat transfer in the heat sink. The two factors of heat transfer and resistance need to be considered comprehensively to determine the optimal flow structure parameters.

CFD Flow and Heat Transfer Simulation for Empty and Packed Fixed Bed Reactor in Catalytic Cracking of Naphtha

This work aims to test the application of computational fluid dynamics (CFD) modeling to fixed bed catalytic cracking reactors. Studies of CFD with a fixed bed design commonly use a regular packing with N=2 to define bed geometry. CFD allows us to obtain a more accurate view of the fluid flow and heat transfer mechanisms present in fixed bed equipment. Naphtha was used as feedstock and the reactor length was 80cm. It is divided in three sections that catalyst bed packed in the middle section of the reactor. The reaction scheme was involved one primary reaction and 24 secondary reactions. Because of high CPU times in these simulations, parallel processing have been used. In this study the coke formation process in fixed bed and empty tube reactor was simulated and coke in these reactors are compared. In addition, the effect of steam ratio and feed flow rate on coke formation was investigated.

Effect of Helium-Argon Mixtures on the Heat Transfer and Fluid Flow in Gas Tungsten Arc Welding

A transient finite element model has been developed to study the heat transfer and fluid flow during spot Gas Tungsten Arc Welding (GTAW) on stainless steel. Temperature field, fluid velocity and electromagnetic fields are computed inside the cathode, arc-plasma and anode using a unified MHD formulation. The developed model is then used to study the influence of different helium-argon gas mixtures on both the energy transferred to the workpiece and the time evolution of the weld pool dimensions. It is found that the addition of helium to argon increases the heat flux density on the weld axis by a factor that can reach 6.5. This induces an increase in the weld pool depth by a factor of 3. It is also found that the addition of only 10% of argon to helium decreases considerably the weld pool depth, which is due to the electrical conductivity of the mixture that increases significantly when argon is added to helium.

Laminar Free Convection of Nanofluid Flow in Horizontal Porous Annulus

A numerical study has been carried out to investigate the heat transfer by natural convection of nanofluid taking Cu as nanoparticles and the water as based fluid in a three dimensional annulus enclosure filled with porous media (silica sand) between two horizontal concentric cylinders with 12 annular fins of 2.4mm thickness attached to the inner cylinder under steady state conditions. The governing equations which used are continuity, momentum and energy equations under an assumptions used Darcy law and Boussinesq-s approximation which are transformed to dimensionless equations. The finite difference approach is used to obtain all the computational results using the MATLAB-7. The parameters affected on the system are modified Rayleigh number (10 ≤Ra*≤ 1000), fin length Hf (3, 7 and 11mm), radius ratio Rr (0.293, 0.365 and 0.435) and the volume fraction(0 ≤ ¤ò ≤ 0 .35). It was found that the average Nusselt number depends on (Ra*, Hf, Rr and φ). The results show that, increasing of fin length decreases the heat transfer rate and for low values of Ra*, decreasing Rr cause to decrease Nu while for Ra* greater than 100, decreasing Rr cause to increase Nu and adding Cu nanoparticles with 0.35 volume fraction cause 27.9% enhancement in heat transfer. A correlation for Nu in terms of Ra*, Hf and φ, has been developed for inner hot cylinder.

Thermal Performance Analysis of Nanofluids in Microchannel Heat Sinks

In the present study, the pressure drop and laminar convection heat transfer characteristics of nanofluids in microchannel heat sink with square duct are numerically investigated. The water based nanofluids created with Al2O3 and CuO particles in four different volume fractions of 0%, 0.5%, 1%, 1.5% and 2% are used to analyze their effects on heat transfer and the pressure drop. Under the laminar, steady-state flow conditions, the finite volume method is used to solve the governing equations of heat transfer. Mixture Model is considered to simulate the nanofluid flow. For verification of used numerical method, the results obtained from numerical calculations were compared with the results in literature for both pure water and the nanofluids in different volume fractions. The distributions of the particles in base fluid are assumed to be uniform. The results are evaluated in terms of Nusselt number, the pressure drop and heat transfer enhancement. Analysis shows that the nanofluids enhance heat transfer while the Reynolds number and the volume fractions are increasing. The best overall enhancement was obtained at φ=%2 and Re=100 for CuO-water nanofluid.

CFD Simulation of Non-Newtonian Fluid Flow in Arterial Stenoses with Surface Irregularities

CFD simulations are carried out in arterial stenoses with 48 % areal occlusion. Non-newtonian fluid model is selected for the blood flow as the same problem has been solved before with Newtonian fluid model. Studies on flow resistance with the presence of surface irregularities are carried out. Investigations are also performed on the pressure drop at various Reynolds numbers. The present study revealed that the pressure drop across a stenosed artery is practically unaffected by surface irregularities at low Reynolds numbers, while flow features are observed and discussed at higher Reynolds numbers.

Investigation on Fluid Flow Characteristics of the Orifice in Nuclear Power Plant

The present paper represents a methodology for investigating flow characteristics near orifice plate by using a commercial computational fluid dynamics code. The flow characteristics near orifice plate which is located in the auxiliary feedwater system were modeled via three different levels of grid and four different types of Reynolds Averaged Navier-Stokes (RANS) equations with proper near-wall treatment. The results from CFD code were compared with experimental data in terms of differential pressure through the orifice plate. In this preliminary study, the Realizable k-ε and the Reynolds stress models with enhanced wall treatment were suitable to analyze flow characteristics near orifice plate, and the results had a good agreement with experimental data.

A New Solution for Natural Convection of Darcian Fluid about a Vertical Full Cone Embedded in Porous Media Prescribed Wall Temperature by using a Hybrid Neural Network-Particle Swarm Optimization Method

Fluid flow and heat transfer of vertical full cone embedded in porous media is studied in this paper. Nonlinear differential equation arising from similarity solution of inverted cone (subjected to wall temperature boundary conditions) embedded in porous medium is solved using a hybrid neural network- particle swarm optimization method. To aim this purpose, a trial solution of the differential equation is defined as sum of two parts. The first part satisfies the initial/ boundary conditions and does contain an adjustable parameter and the second part which is constructed so as not to affect the initial/boundary conditions and involves adjustable parameters (the weights and biases) for a multi-layer perceptron neural network. Particle swarm optimization (PSO) is applied to find adjustable parameters of trial solution (in first and second part). The obtained solution in comparison with the numerical ones represents a remarkable accuracy.

Numerical Investigation of Nozzle Shape Effect on Shock Wave in Natural Gas Processing

Natural gas flow contains undesirable solid particles, liquid condensation, and/or oil droplets and requires reliable removing equipment to perform filtration. Recent natural gas processing applications are demanded compactness and reliability of process equipment. Since conventional means are sophisticated in design, poor in efficiency, and continue lacking robust, a supersonic nozzle has been introduced as an alternative means to meet such demands. A 3-D Convergent-Divergent Nozzle is simulated using commercial Code for pressure ratio (NPR) varies from 1.2 to 2. Six different shapes of nozzle are numerically examined to illustrate the position of shock-wave as such spot could be considered as a benchmark of particle separation. Rectangle, triangle, circular, elliptical, pentagon, and hexagon nozzles are simulated using Fluent Code with all have same cross-sectional area. The simple one-dimensional inviscid theory does not describe the actual features of fluid flow precisely as it ignores the impact of nozzle configuration on the flow properties. CFD Simulation results, however, show that nozzle geometry influences the flow structures including location of shock wave. The CFD analysis predicts shock appearance when p01/pa>1.2 for almost all geometry and locates at the lower area ratio (Ae/At). Simulation results showed that shock wave in Elliptical nozzle has the farthest distance from the throat among the others at relatively small NPR. As NPR increases, hexagon would be the farthest. The numerical result is compared with available experimental data and has shown good agreement in terms of shock location and flow structure.

Transient Hydrodynamic and Thermal Behaviors of Fluid Flow in a Vertical Porous Microchannel under the Effect of Hyperbolic Heat Conduction Model

The transient hydrodynamics and thermal behaviors of fluid flow in open-ended vertical parallel-plate porous microchannel are investigated semi-analytically under the effect of the hyperbolic heat conduction model. The model that combines both the continuum approach and the possibility of slip at the boundary is adopted in the study. The Effects of Knudsen number , Darcy number , and thermal relaxation time  on the microchannel hydrodynamics and thermal behaviors are investigated using the hyperbolic heat conduction models. It is found that as  increases the slip in the hydrodynamic and thermal boundary condition increases. This slip in the hydrodynamic boundary condition increases as  increases. Also, the slip in the thermal boundary condition increases as  decreases especially the early stage of time.

Simulation of Fluid Flow and Heat Transfer in Inclined Cavity using Lattice Boltzmann Method

In this paper, Lattice Boltzmann Method (LBM) is used to study laminar flow with mixed convection heat transfer inside a two-dimensional inclined lid-driven rectangular cavity with aspect ratio AR = 3. Bottom wall of the cavity is maintained at lower temperature than the top lid, and its vertical walls are assumed insulated. Top lid motion results in fluid motion inside the cavity. Inclination of the cavity causes horizontal and vertical components of velocity to be affected by buoyancy force. To include this effect, calculation procedure of macroscopic properties by LBM is changed and collision term of Boltzmann equation is modified. A computer program is developed to simulate this problem using BGK model of lattice Boltzmann method. The effects of the variations of Richardson number and inclination angle on the thermal and flow behavior of the fluid inside the cavity are investigated. The results are presented as velocity and temperature profiles, stream function contours and isotherms. It is concluded that LBM has good potential to simulate mixed convection heat transfer problems.

Evaluation of Stiffness and Damping Coefficients of Multiple Axial Groove Water Lubricated Bearing Using Computational Fluid Dynamics

This research details a Computational Fluid Dynamics (CFD) approach to model fluid flow in a journal bearing with 8 equispaced semi-circular axial grooves. Water is used as the lubricant and is fed from one end of the bearing to the other, under pressure. The geometry of the bearing is modeled using a commercially available modeling software GAMBIT and the flow analysis is performed using a dedicated CFD analysis software FLUENT. The pressure distribution in the bearing clearance is obtained from FLUENT for various whirl ratios and is used to calculate the hydrodynamic force components in the radial and tangential direction of the bearing. These values along with the various whirl speeds can be used to do a regression analysis to determine the stiffness and damping coefficients. The values obtained are then compared with the stiffness and damping coefficients of a 3 Axial groove water lubricated journal bearing and those obtained from a FORTRAN code for a similar bearing.

CFD Analysis of Incompressible Turbulent Swirling Flow through Circle Grids Space Filling Plate

Circle grid space filling plate is a flow conditioner with a fractal pattern and used to eliminate turbulence originating from pipe fittings in experimental fluid flow applications. In this paper, steady state, incompressible, swirling turbulent flow through circle grid space filling plate has been studied. The solution and the analysis were carried out using finite volume CFD solver FLUENT 6.2. Three turbulence models were used in the numerical investigation and their results were compared with the pressure drop correlation of BS EN ISO 5167-2:2003. The turbulence models investigated here are the standard k-ε, realizable k-ε, and the Reynolds Stress Model (RSM). The results showed that the RSM model gave the best agreement with the ISO pressure drop correlation. The effects of circle grids space filling plate thickness and Reynolds number on the flow characteristics have been investigated as well.

Pulsating Flow of an Incompressible Couple Stress Fluid Between Permeable Beds

The paper deals with the pulsating flow of an incompressible couple stress fluid between permeable beds. The couple stress fluid is injected into the channel from the lower permeable bed with a certain velocity and is sucked into the upper permeable bed with the same velocity. The flow between the permeable beds is assumed to be governed by couple stress fluid flow equations of V. K. Stokes and that in the permeable regions by Darcy-s law. The equations are solved analytically and the expressions for velocity and volume flux are obtained. The effects of the material parameters are studied numerically and the results are presented through graphs.

Numerical Simulation of Heat Transfer in Primary Surface with Corrugations Recuperators

Study fluid flow and heat transfer characteristics of microchannel in a primary Cross-corrugated(CC) surface recuperators with corrugations and without corrugations, using CFD method. The pitch-over-height ratios P/H of Cross-corrugated (CC) surface is from 1.5 to 4.0, included angles β=75º. The study was performed using CFD software FLUENT to create unit model and simulate fluid temperature, velocity, heat transfer coefficient and other parameters. The results from these simulations were compared to experimental data. It is concluded that, when the Reynolds number is constant, if increase P/H, j/f will decrease, also the decreasing trend will become weak. Under the condition of P/H=2.2, if increase the inlet velocity j/f will decrease; in addition, the heat transfer performance in surface with corrugation will increase 10% compared to that without corrugation. The study results can provide the basis to optimize the design, select the type of heat transfer surface, the scale structure, and heat-transfer surface arrangement for recuperators.

A Meshfree Solution of Tow-Dimensional Potential Flow Problems

In this paper, mesh-free element free Galerkin (EFG) method is extended to solve two-dimensional potential flow problems. Two ideal fluid flow problems (i.e. flow over a rigid cylinder and flow over a sphere) have been formulated using variational approach. Penalty and Lagrange multiplier techniques have been utilized for the enforcement of essential boundary conditions. Four point Gauss quadrature have been used for the integration on two-dimensional domain (Ω) and nodal integration scheme has been used to enforce the essential boundary conditions on the edges (┌). The results obtained by EFG method are compared with those obtained by finite element method. The effects of scaling and penalty parameters on EFG results have also been discussed in detail.

Effects of Operating Conditions on Calcium Carbonate Fouling in a Plate Heat Exchanger

The aim of this work is to investigate on the internalflow patterns in a plate heat exchanger channel, which affect the rate of sedimentation fouling on the heat transfer surface of the plate heat exchanger. The research methodologies were the computer simulation using Computational Fluid Dynamics (CFD) and the experimental works. COMSOL MULTIPHYSICS™ Version 3.3 was used to simulate the velocity flow fields to verify the low and high flow regions. The results from the CFD technique were then compared with the images obtained from the experiments in which the fouling test rig was set up with a singlechannel plate heat exchanger to monitor the fouling of calcium carbonate. Two parameters were varied i.e., the crossing angle of the two plate: 55/55, 10/10, and 55/10 degree, and the fluid flow rate at the inlet: 0.0566, 0.1132 and 0.1698 m/s. The type of plate “GX-12" (the surface area 0.12 m2, the depth 2.9 mm, the width of fluid flow 215 mm and the thickness of stainless plate of 0.5 mm) was used in this study. The results indicated that the velocity distribution for the case of 55/55 degree seems to be very well organized when compared with the others. Also, an increase in the inlet velocity resulted in the reduction of fouling rate on the surface of plate heat exchangers.

Mathematical Modelling of Partially Filled Fluid Coupling Behaviour

Modelling techniques for a fluid coupling taken from published literature have been extended to include the effects of the filling and emptying of the coupling with oil and the variation in losses when the coupling is partially full. In the model, the fluid flow inside the coupling is considered to have two principal velocity components; one circumferentially about the coupling axis (centrifugal head) and the other representing the secondary vortex within the coupling itself (vortex head). The calculation of liquid mass flow rate circulating between the two halves of the coupling is based on: the assumption of a linear velocity variation in the circulating vortex flow; the head differential in the fluid due to the speed difference between the two shafts; and the losses in the circulating vortex flow as a result of the impingement of the flow with the blades in the coupling and friction within the passages between the blades.

Numerical Simulations of Shear Driven Square and Triangular Cavity by Using Lattice Boltzmann Scheme

In this paper, fluid flow patterns of steady incompressible flow inside shear driven cavity are studied. The numerical simulations are conducted by using lattice Boltzmann method (LBM) for different Reynolds numbers. In order to simulate the flow, derivation of macroscopic hydrodynamics equations from the continuous Boltzmann equation need to be performed. Then, the numerical results of shear-driven flow inside square and triangular cavity are compared with results found in literature review. Present study found that flow patterns are affected by the geometry of the cavity and the Reynolds numbers used.