Thermal entry flow problem for Giesekus fluid inside an axis-symmetric tube through isothermal wall condition: a comparative numerical study between exact and approximate solution

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Muhammad Waris Saeed Khan ◽  
Nasir Ali
Keyword(s):  
Nusselt Number ◽  
Energy Equation ◽  
Flow Problem ◽  
Numerical Study ◽  
Fluid Model ◽  
Weissenberg Number ◽  
Mean Flow ◽  
Mobility Parameter ◽  
Giesekus Fluid ◽  
The Impact

Abstract The thermal entry flow problem also known as the Graetz problem is investigated for a Giesekus fluid model. Both analytical (exact) and approximate solutions for velocity are obtained. The nondimensional pressure gradient is numerically obtained via the mean flow rate relation. The energy equation along with the Giesekus fluid velocity is analytically solved for the constant wall temperature case by using the classical separation of variable method. This method transforms the energy equation into a Sturm–Liouville (SL) boundary value problem. The MATLAB solver bvp5c is employed to compute the eigenvalues and the related eigenfunctions numerically. The impact of mobility parameter and Weissenberg number on local Nusselt number, mean temperature, and average Nusselt number is discussed and displayed graphically. It is also found that the presence of the Weissenberg number elevates the Nusselt numbers. Further, the presence of the mobility parameter of the Giesekus fluid model delays the prevalence fully developed conditions in both entrance and fully developed regions. The comparison between approximate and exact solution is also presented. It reveals that both solutions have an exact match with each other for smaller values of mobility parameter and Weissenberg number. However, there is a deviation for larger values of both parameters.

scholarly journals A numerical study on the impact of nonlinear interactions on the amplitude of the migrating semidiurnal tide

2006 ◽  
Vol 24 (12) ◽  
pp. 3241-3256 ◽  
Author(s):  
C. M. Huang ◽  
S. D. Zhang ◽  
F. Yi
Keyword(s):  
Numerical Study ◽  
Steady Solution ◽  
Nonlinear Interactions ◽  
Semidiurnal Tide ◽  
Mean Flow ◽  
Structure Variation ◽  
Amplitude Difference ◽  
The Mean ◽  
The Impact ◽  
Mlt Region

Abstract. To quantitatively study the effects of nonlinear interactions on tide structure, a nonlinear numerical tidal model is developed, and the reliability and convergence of the adopted algorithm and coding are checked by numerical experiments. Under the same conditions as those employed by the GSWM-00 (Global Scale Wave Model 2000), our model provides the nonlinear quasi-steady solution of the migrating semidiurnal tide, which differs from the GSWM-00 result (the linear steady solution) in the MLT region, especially above 100 km. Additionally, their amplitude difference displays a remarkable month-to-month variation, and its significant magnitudes occur during the month with strong semidiurnal tide. A quantitative analysis suggests that the main cause for the amplitude difference is that the initial migrating 12-h tide will interact with the mean flow as well as the nonlinearity-excited 6-h tide, and subsequently yield a new 12-h tidal part. Furthermore, our simulations also show that the mean flow/tidal interaction will significantly alter the background wind and temperature fields. The large magnitudes of the tidal amplitude difference and the background alteration indicate that the nonlinear processes involved in tidal propagations should be comprehensively considered in the description of global atmospheric dynamics in the MLT region. The comparisons among our simulations, the GSWMs and some observations of tides suggest that the nonlinearity-induced tidal structure variation could be a possible mechanism to account for some discrepancies between the GSWMs and the observations.

Numerical approach for the calendering process using Carreau-Yasuda fluid model

2021 ◽  
pp. 875608792098874
Author(s):  
Muhammad Asif Javed ◽  
Nasir Ali ◽  
Sabeen Arshad ◽  
Shahbaz Shamshad
Keyword(s):  
Stream Function ◽  
Power Law ◽  
Numerical Study ◽  
Fluid Model ◽  
Pressure Profile ◽  
Numerical Approach ◽  
Weissenberg Number ◽  
Model Parameters ◽  
Flow Equations ◽  
Power Law Index

This paper presents a numerical study of the calendering mechanism. The calendered material is represented using the Carreau-Yasuda fluid model. The governing flow equations in the calendering process are made first dimensionless then the lubrication approximation theory (LAT) is used to simplify them. The simplified flow equations are transformed into stream function and then are numerically solved. A numerical method is constructed with Matlab’s built-in-bvp4c routine to find the stream function and pressure gradient. We use the Runge-Kutta algorithm to calculate the pressure and mechanical quantities related to the calendering process. In this analysis the pressure distribution increases with increasing Weissenberg number, however the pressure domain length decreases as the Weissenberg number increases. The pressure inside the nip region decreases from its Newtonian value when the power law index is less than one (shear thinning), and the pressure profile increases from its Newtonian pressure when the power law index is greater than one(shear thickening). How the Carreau-Yasuda fluid model parameters influence the velocity and related calendering process quantities are also discussed via graphs.

scholarly journals Numerical Study on Natural Convection of Alumina-water Nanofluid in a Square Cavity with Two Localized Heat Sources on Adjacent Surfaces

2021 ◽  
Author(s):  
Farooq Shaik ◽  
Vinay Kumar Domakonda
Keyword(s):  
Heat Transfer ◽  
Nusselt Number ◽  
Natural Convection ◽  
Vertical Velocity ◽  
Volume Fraction ◽  
Numerical Study ◽  
Al2o3 Nanoparticles ◽  
Heat Sources ◽  
Square Cavity ◽  
The Impact

In the present work, results of a numerical study carried out using finite volume method, to investigate the fluid flow and heat transfer characteristics of Alumina ( Al2O3 ) nanoparticles in the base fluid (water) in a square cavity under natural convection mode are presented. The Semi Implicit Method for Pressure Linked Equations (SIMPLE) algorithm was used to solve the discretized momentum and energy equations. Constant temperature heat sources of same strength are placed on bottom and left vertical surfaces whereas the right surface was kept cold, while the top surface was maintained as adiabatic. The impact of Rayleigh number (RaN) ( 1000 to 106 ) and nanoparticles volume fraction (Φ = 0 %, 5 %, 10 %, 15 % and 20 %) on fluid and heat flow characteristics were numerically investigated and presented in the form of streamlines, isothermal lines, mid line horizontal and vertical velocity components, local Nusselt number ( Nuloc ) and average Nusselt number ( Nuavg ). The obtained results indicate, for lower RaN ( i.e; 103 ), conduction dominates over convection near heated surfaces and results in lower fluid velocities and poor heat transfer. For higher values of RaN ( RaN = 105 and 106 ) and volume fraction of nanoparticles, there was a significant increase in mid horizontal and vertical velocity components, Nuloc and Nuavg due to increase in convective heat transfer and thermal conductivity of nanofluid.

Numerical analysis of the calendering process by using Giesekus fluid model

2019 ◽  
Vol 36 (2) ◽  
pp. 167-190 ◽  
Author(s):  
Muhammad Asif Javed ◽  
Nasir Ali ◽  
Sabeen Arshad
Keyword(s):  
Pressure Gradient ◽  
Power Function ◽  
Numerical Study ◽  
Fluid Model ◽  
Sheet Thickness ◽  
Deborah Number ◽  
Flow Equations ◽  
Giesekus Fluid ◽  
Engineering Parameters ◽  
Roll Separating Force

A numerical study of the calendering process is presented. The material to be calendered is modeled by using Giesekus constitutive equation. The flow equations are first presented in dimensionless forms and then simplified by incorporating the lubrication approximation theory. The resulting equations are analytically solved for the stream function. The pressure gradient, pressure, and other engineering parameters related to the calendering process, such as roll-separating force, power function, and entering sheet thickness, are numerically calculated by using Runge–Kutta algorithm. The influence of the Giesekus parameter and the Deborah number on the velocity profile, pressure gradient, pressure, power function, roll-separating force, and exiting sheet thickness are discussed in detail with the help of various graphs. The present analysis indicates that the pressure in the nip region decreases with increasing Giesekus parameter and Deborah number. The power function and the roll-separating force exhibit decreasing trends with increasing Deborah number. The exiting sheet thickness decreases up to a certain entering sheet thickness, as compared to the Newtonian case. Beyond this entering sheet thickness, the exiting sheet thickness increases with increasing entering sheet thickness.

Effect of Viscous Dissipation on MHD Williamson Nanofluid Flow in a Porous Medium

2019 ◽  
Vol 8 (3) ◽  
pp. 5795-5802 ◽  
Keyword(s):  
Porous Medium ◽  
Steady State ◽  
Numerical Solution ◽  
Viscous Dissipation ◽  
Flow Problem ◽  
Numerical Study ◽  
Steady State Flow ◽  
Nanofluid Flow ◽  
Shooting Technique ◽  
Order Four

The main objective of this paper is to focus on a numerical study of viscous dissipation effect on the steady state flow of MHD Williamson nanofluid. A mathematical modeled which resembles the physical flow problem has been developed. By using an appropriate transformation, we converted the system of dimensional PDEs (nonlinear) into coupled dimensionless ODEs. The numerical solution of these modeled ordinary differential equations (ODEs) is achieved by utilizing shooting technique together with Adams-Bashforth Moulton method of order four. Finally, the results of discussed for different parameters through graphs and tables.

scholarly journals On the flow of MHD generalized maxwell fluid via porous rectangular duct

Open Physics ◽  
2020 ◽  
Vol 18 (1) ◽  
pp. 989-1002
Author(s):  
Aamir Farooq ◽  
Muhammad Kamran ◽  
Yasir Bashir ◽  
Hijaz Ahmad ◽  
Azeem Shahzad ◽  
...  
Keyword(s):  
Fluid Model ◽  
Maxwell Fluid ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Rectangular Duct ◽  
Modified Darcy’S Law ◽  
Initial Boundary Value ◽  
Generalized Maxwell Fluid ◽  
Limiting Case ◽  
The Impact

Abstract The purpose of this proposed investigation is to study unsteady magneto hydrodynamic (MHD) mixed initial-boundary value problem for incompressible fractional Maxwell fluid model via oscillatory porous rectangular duct. Considering the modified Darcy’s law, the problem is simplified by using the method of the double finite Fourier sine and Laplace transforms. As a limiting case of the general solutions, the same results can be obtained for the classical Maxwell fluid. Also, the impact of magnetic parameter, porosity of medium, and the impact of various material parameters on the velocity profile and the corresponding tangential tensions are illuminated graphically. At the end, we will give the conclusion of the whole paper.

scholarly journals Model of thermal buoyancy in cavity-ventilated roof constructions

2021 ◽  
pp. 174425912098418
Author(s):  
Toivo Säwén ◽  
Martina Stockhaus ◽  
Carl-Eric Hagentoft ◽  
Nora Schjøth Bunkholt ◽  
Paula Wahlgren
Keyword(s):  
Analytical Model ◽  
Flow Rate ◽  
Numerical Study ◽  
Air Flow ◽  
Wind Pressure ◽  
Air Flow Rate ◽  
Test Set ◽  
Air Cavity ◽  
Thermal Buoyancy ◽  
The Impact

Timber roof constructions are commonly ventilated through an air cavity beneath the roof sheathing in order to remove heat and moisture from the construction. The driving forces for this ventilation are wind pressure and thermal buoyancy. The wind driven ventilation has been studied extensively, while models for predicting buoyant flow are less developed. In the present study, a novel analytical model is presented to predict the air flow caused by thermal buoyancy in a ventilated roof construction. The model provides means to calculate the cavity Rayleigh number for the roof construction, which is then correlated with the air flow rate. The model predictions are compared to the results of an experimental and a numerical study examining the effect of different cavity designs and inclinations on the air flow rate in a ventilated roof subjected to varying heat loads. Over 80 different test set-ups, the analytical model was found to replicate both experimental and numerical results within an acceptable margin. The effect of an increased total roof height, air cavity height and solar heat load for a given construction is an increased air flow rate through the air cavity. On average, the analytical model predicts a 3% higher air flow rate than found in the numerical study, and a 20% lower air flow rate than found in the experimental study, for comparable test set-ups. The model provided can be used to predict the air flow rate in cavities of varying design, and to quantify the impact of suggested roof design changes. The result can be used as a basis for estimating the moisture safety of a roof construction.

scholarly journals A numerical study of gravity-driven instability in strongly coupled dusty plasma. Part 1. Rayleigh–Taylor instability and buoyancy-driven instability

2021 ◽  
Vol 87 (2) ◽  
Author(s):  
Vikram S. Dharodi ◽  
Amita Das
Keyword(s):  
Dusty Plasma ◽  
Numerical Study ◽  
Fluid Model ◽  
Density Region ◽  
Strongly Coupled ◽  
Rising Bubble ◽  
Inhomogeneous Density ◽  
Gravity Driven ◽  
The Individual ◽  
Strongly Coupled Dusty Plasma

Rayleigh–Taylor (RT) and buoyancy-driven (BD) instabilities are driven by gravity in a fluid system with inhomogeneous density. The paper investigates these instabilities for a strongly coupled dusty plasma medium. This medium has been represented here in the framework of the generalized hydrodynamics (GHD) fluid model which treats it as a viscoelastic medium. The incompressible limit of the GHD model is considered here. The RT instability is explored both for gradual and sharp density gradients stratified against gravity. The BD instability is discussed by studying the evolution of a rising bubble (a localized low-density region) and a falling droplet (a localized high-density region) in the presence of gravity. Since both the rising bubble and falling droplet have symmetry in spatial distribution, we observe that a falling droplet process is equivalent to a rising bubble. We also find that both the gravity-driven instabilities get suppressed with increasing coupling strength of the medium. These observations have been illustrated analytically as well as by carrying out two-dimensional nonlinear simulations. Part 2 of this paper is planned to extend the present study of the individual evolution of a bubble and a droplet to their combined evolution in order to understand the interaction between them.

scholarly journals Numerical Simulation of the Impact of the Heat Source Position on Melting of a Nano-Enhanced Phase Change Material

Nanomaterials ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 1425
Author(s):  
Tarek Bouzennada ◽  
Farid Mechighel ◽  
Kaouther Ghachem ◽  
Lioua Kolsi
Keyword(s):  
Heat Transfer ◽  
Phase Change ◽  
Phase Change Material ◽  
Numerical Study ◽  
Melting Rate ◽  
Heat Transfer Fluid ◽  
Commercial Code ◽  
Tube Position ◽  
The Impact ◽  
Change Material

A 2D-symmetric numerical study of a new design of Nano-Enhanced Phase change material (NEPCM)-filled enclosure is presented in this paper. The enclosure is equipped with an inner tube allowing the circulation of the heat transfer fluid (HTF); n-Octadecane is chosen as phase change material (PCM). Comsol-Multiphysics commercial code was used to solve the governing equations. This study has been performed to examine the heat distribution and melting rate under the influence of the inner-tube position and the concentration of the nanoparticles dispersed in the PCM. The inner tube was located at three different vertical positions and the nanoparticle concentration was varied from 0 to 0.06. The results revealed that both heat transfer/melting rates are improved when the inner tube is located at the bottom region of the enclosure and by increasing the concentration of the nanoparticles. The addition of the nanoparticles enhances the heat transfer due to the considerable increase in conductivity. On the other hand, by placing the tube in the bottom area of the enclosure, the liquid PCM gets a wider space, allowing the intensification of the natural convection.

scholarly journals Generation of Reverse Meniscus Flow by Applying An Electromagnetic Brake

2021 ◽  
Author(s):  
Alexander Vakhrushev ◽  
Abdellah Kharicha ◽  
Ebrahim Karimi-Sibaki ◽  
Menghuai Wu ◽  
Andreas Ludwig ◽  
...  
Keyword(s):  
Magnetic Field ◽  
Lorentz Force ◽  
Applied Magnetic Field ◽  
Numerical Study ◽  
Flow Direction ◽  
Experimental Studies ◽  
Submerged Entry Nozzle ◽  
Mean Flow ◽  
Slab Caster ◽  
The Mean

AbstractA numerical study is presented that deals with the flow in the mold of a continuous slab caster under the influence of a DC magnetic field (electromagnetic brakes (EMBrs)). The arrangement and geometry investigated here is based on a series of previous experimental studies carried out at the mini-LIMMCAST facility at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR). The magnetic field models a ruler-type EMBr and is installed in the region of the ports of the submerged entry nozzle (SEN). The current article considers magnet field strengths up to 441 mT, corresponding to a Hartmann number of about 600, and takes the electrical conductivity of the solidified shell into account. The numerical model of the turbulent flow under the applied magnetic field is implemented using the open-source CFD package OpenFOAM®. Our numerical results reveal that a growing magnitude of the applied magnetic field may cause a reversal of the flow direction at the meniscus surface, which is related the formation of a “multiroll” flow pattern in the mold. This phenomenon can be explained as a classical magnetohydrodynamics (MHD) effect: (1) the closure of the induced electric current results not primarily in a braking Lorentz force inside the jet but in an acceleration in regions of previously weak velocities, which initiates the formation of an opposite vortex (OV) close to the mean jet; (2) this vortex develops in size at the expense of the main vortex until it reaches the meniscus surface, where it becomes clearly visible. We also show that an acceleration of the meniscus flow must be expected when the applied magnetic field is smaller than a critical value. This acceleration is due to the transfer of kinetic energy from smaller turbulent structures into the mean flow. A further increase in the EMBr intensity leads to the expected damping of the mean flow and, consequently, to a reduction in the size of the upper roll. These investigations show that the Lorentz force cannot be reduced to a simple damping effect; depending on the field strength, its action is found to be topologically complex.

Sign in / Sign up

Export Citation Format

Share Document