Simulation of Non-isothermal Injection Molding for a Non-Newtonian Fluid by Level Set Method

A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic simulation.The validity of the numerical method is verified by a benchmark problem.The melt interface evolution versus time is captured and the physical quantities such as temperature,velocity and pressure at each time step are obtained with corresponding analysis.A"frozen skin"layer with the thickness increasing versus time during the injection process is found.The fact that the"frozen skin"layer can be reduced by increasing the injection velocity is numerically verified.The fountain flow phenomenon near the melt interface is also captured.Moreover,comparisons with the non-isothermal Newtonian case show that the curvatures of the interface arcs and the pressure contours near the horizontal mid-line of the cavity for the non-Newtonian pseudoplastic case is larger than that for the Newtonian case.The velocity profiles are different at different positions for the non-Newtonian pseudoplastic case,while in the case of Newtonian flow the velocity profiles are parabolic and almost the same at different positions.

THE NUMERICAL SIMULATION OF TURBULENT FLOWS OF NEWTONIAN AND A SORT OF NON-NEWTONIAN FLUID IN 180-DEGREE SQUARE SECTIONED BEND

Using k- model of turbulence and measured wall functions, turbulent flows of Newtonian (pure water) andasort of non-Newtonian fluid (dilute, drag-reduction solution of polymer) in a 180-degree curved bend were simulated numerically. The calculated results agreed well with the measured velocity profiles. On the basis of calculation and measurement, appropriateness of turbulence model to complicated flow in which the large-scale vortex exists was analyzed and discussed.

Implementation of Non-Newtonian Fluid Properties for Compressible Multiphase Flows in OpenFOAM

The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.

Jeffery-Hamel flow of non-Newtonian fluid with nonlinear viscosity and wall friction

A Jeffery-Hamel （J-H） flow model of the non-Newtonian fluid type inside a convergent wedge （inclined walls） with a wall friction is derived by a nonlinear ordinary differential equation with appropriate boundary conditions based on similarity relationships. Unlike the usual power law model, this paper develops nonlinear viscosity based only on a tangential coordinate function due to the radial geometry shape. Two kinds of solutions are developed, i.e., analytical and semi-analytical （numerical） solutions with suitable assumptions. As a result of the parametric examination, it has been found that the Newtonian normalized velocity gradually decreases with the tangential direction progress. Also, an increase in the friction coefficient leads to a decrease in the normalized Newtonian velocity profile values. However, an increase in the Reynolds number causes an increase in the normalized velocity function values. Additionally, for the small values of wedge semi-angle, the present solutions are in good agreement with the previous results in the literature.

ANALYTICAL SOLUTIONS FOR EQUATIONS OF UNSTEADY FLOW OF NON-NEWTONIAN FLUIDS IN TUBE

This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of nun-Newtonian visco-elastie fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.

Unsteady Rotative Flow of Non-Newtonian Fluid in an Annular Pipe

This paper .Studies power law no-Newtonian fluid rotative flow. in an annularpipe. The governing equation is nonlinear one, we linearized the governing equationby assuming that partial factor is at state. With Laplace transform we obtain ananalytical solution of the problem In the paper several groups of curves are given.these curves reflect the temporal change law and. spatial distribution of fluid velocity.In addition.we study the effection of power law index on the flow field the resultindicates that when the power law index n < l. the flow velocity is highly sensitive tothe index. and this fact is importanl in related engineering decisions.

A level set based immersed boundary method for simulation of non-isothermal viscoelastic melt filling process

In this work,the polymer melt filling process is simulated by using a coupled finite volume and levelset based immersed boundary(LSIB)method.Firstly,based on a shape level set(LS)function to represent the mold boundary,a LSIB method is developed to model the complex mold walls.Then the nonisothermal melt filling process is simulated based on nonNewtonian viscoelastic equations with different Reynolds numbers in a circular cavity with a solid core,and the effects of Reynolds number on the flow patterns of polymer melt are presented and compared with each other.And then for a true polymer melt with a small Reynolds number that varies with melt viscosity,the moving interface,the temperature distributions and the molecular deformation are shown and analyzed in detail.At last,as a commonly used application case,a socket cavity with seven inserts is investigated.The corresponding physical quantities,such as the melt velocity,molecular deformation,normal stresses,first normal stress difference,temperature distributions and frozen layer are analyzed and discussed.The results could provide some predictions and guidance for the polymer processing industry.

Modeling study on fluid flow and inclusion motion in 6-strand bloom caster tundishes

The behavior of fluid flow and particle motion in a 6-strand bloom castertundish was investigated by a water model and numerical simulation. Compared with a device withoutflow control, the tundish with flow control has an important effect on the fluid flow pattern andinclusion removal. It is revealed that by non-isothermal process, which is real productioncondition, the fluid flow in tundish shows a strong buoyancy pattern, which drives particles to moveupwards. The particle removal was quantitatively studied by mathematical and physical simulations.

COMPUTATIONALINTELLECTUAL ANALYTICAL THEORYOF COMPUTATIONAL ANALYTICAL APPROACH TOROTATING FLOW OF NON_NEWTONIAN FLUID

A combination of the computational symbolic calculation, mathematical approach and physico-mechanical model lends to a computational intellectual analytical approach developed by the author. There is a principal difference between the computer proof and the computer derivation completed by the computer, also difference between the numerical and symbolic calculations. In this investigation the computational analytical approach is extended, and an unsteady flow of non-Newtonian fluid in the gap between two rotating coaxial cylinders is studied. The Oldroyd fluid B model is used by which the Weissenberg effects are explained in a good comparison with the experiments. The governing equations are reduced to a partial differential equation of 3 rd order for the dimensionless velocity. Using the computer software Macsyma and an improved variational approach the problem with the initial and boundary conditions is then reduced to a problem of an ordinary differential equation for different approximations. The analytical solutions are given for the 1 st, 2 nd and 3 rd approximations. The present investigation shows the ability of the computational symbolic manipulation in solving the problems of non-Newtonian fluid flows. There is a possibility of that to solve the problems in mathematics and mechanics. An important conclusion can be drawn from the results that the transition from a steady state to another steady state is non-unique.

Time Dependent Pressure Gradient Effect on Unsteady MHD Couette Flow and Heat Transfer of a Casson Fluid

The unsteady magnetohydrodynamic flow of an electrically conducting viscous incompressible non-Newto- nian Casson fluid bounded by two parallel non-conducting porous plates has been studied with heat transfer considering the Hall effect. The fluid is acted upon by a uniform and exponential decaying pressure gradient. An external uniform magnetic field is applied perpendicular to the plates and the fluid motion is subjected to a uniform suction and injection. The lower plate is stationary and the upper plate is suddenly set into mo- tion and simultaneously suddenly isothermally heated to a temperature other than the lower plate temperature. Numerical solutions are obtained for the governing momentum and energy equations taking the Joule and viscous dissipations into consideration. The effect of unsteady pressure gradient, the Hall term, the parameter describing the non-Newtonian behavior on both the velocities and temperature distributions have been stud- ied.

Modeling and Numerical Simulation of Yield Viscoplastic Fluid Flow in Concentric and Eccentric Annuli

Numerical solution of yield viscoplastic fluid flow is hindered by the singularity inherent to the Herschel-Bulkley model. A finite difference method over the boundary-fitted orthogonal coordinate system is util- ized to investigate numerically the fully developed steady flow of non-Newtonian yield viscoplastic fluid through concentric and eccentric annuli. The fluid rheology is described with the Herschel-Bulkley model. The numerical simulation based on a continuous viscoplastic approach to the Herschel-Bulkley model is found in poor accordance with the experimental data on volumetric flow rate of a bentonite suspension. A strict mathematical model for Herschel-Bulkley fluid flow is established and the corresponding numerical procedures are proposed. However, only the case of flow of a Herschel-Bulkley fluid in a concentric annulus is resolved based on the presumed flow stnicture by using the common optimization technique. Possible flow structures in an eccentric afinulus are presumed, and further challenges in numerical simulation of the Herschel-Bulkley fluid flow are suggested.

Dynamics of self-organizing single-line particle trains in the channel flow of a power-law fluid

The formation of self-organizing single-line particle train in a channel flow of a power-law fluid is studied using the lattice Boltzmann method with power-law index 0.6≤n≤1.2,particle volume concentration 0.8%≤Φ≤6.4%,Reynolds number 10≤Re≤100,and blockage ratio 0.2≤k≤0.4.The numerical method is validated by comparing the present results with the previous ones.The effect n,Φ,Re and k on the interparticle spacing and parallelism of particle train is discussed.The results showed that the randomly distributed particles would migrate towards the vicinity of the equilibrium position and form the ordered particle train in the power-law fluid.The equilibrium position of particles is closer to the channel centerline in the shear-thickening fluid than that in the Newtonian fluid and shear-thinning fluid.The particles are not perfectly parallel in the equilibrium position,hence IH is used to describe the inclination of the line linking the equilibrium position of each particle.When self-organizing single-line particle train is formed,the particle train has a better parallelism and hence benefit for particle focusing in the shearthickening fluid at highΦ,low Re and small k.Meanwhile,the interparticle spacing is the largest and hence benefit for particle separation in the shear-thinning fluid at lowΦ,low Re and small k.

COMPUTER SIMULATION OF NON-NEWTONIAN FLOW AND MASS TRANSPORT THROUGH CORONARY ARTERIAL STENOSIS

A numerical analysis of Newtonian and non-Newtonian flow in an axi-symmetric tube with a local constriction simulating a stenosed artery under steady and pulsatile flow conditions war carried out. Bared on these results, the concentration fields of LDL ( (low-density lipoprotein) and Albumin were discussed. According to the results, in great details the macromolecule transport influences of wall shear stress, non-Newtonian fluid character and the scale of the molecule etc are given. The results of Newtonian fluid flow and non-Newtonian fluid flow, steady flow and pulsatile flow are compared. These investigations can provide much valuable information about the correlation between the flow properties, the macromolecule transport and the development of atherosclerosis.

Turbulent Characteristic of Liquid Around a Chain of Bubbles in Non-Newtonian Fluid

The turbulence behavior of gas-liquid two-phase flow plays an important role in heat transfer and mass transfer in many chemical processes. In this work, a 2D particle image velocimetry (PIV) was used to investigate the turbulent characteristic of fluid induced by a chain of bubbles rising in Newtonian and non-Newtonian fluids. The instantaneous flow field, turbulent kinetic energy (TKE) and TKE dissipation rate were measured. The results demonstrated that the TKE profiles were almost symmetrical along the column center and showed higher values in the central region of the column. The TKE was enhanced with the increase of gas flow and decrease of liquid viscosity. The maximum TKE dissipation rate appeared on both sides of the bubble chain, and increased with the increase of gas flow rate or liquid viscosity. These results provide an understanding for gas-liquid mass transfer in non-Newtonian fluids.

THE NONLINEAR STABILITY OF PLANE POISEUILLE FLOW FOR NON－NEWTONIAN POWER SERIES FLUID

Non-Newtonan fluid is a kind of fluid whose components of stresstensor aren’t theliear funtions of compoents of the strain rate tensor.Non-Newtonianfluid is beingprocessed in many kinds of modern industry,Stability of flows for Non- Newtonianfluid is of important applicatuib,In this article we calculate subcritical thrdshold of flow which oecurs in polymer-processing when the melting substance is driven throughtwo parallel fixed boundaries.

Analytical Study of the Electroosmotic Flow of Two Immiscible Power-Law Fluids in a Microchannel

The multilayer microchannel flow is a promising tool in microchannel-based systems such as hybrid microfluidics. To assist in the efficient design of two-liquid pumping system, a two-fluid electroosmotic flow of immiscible power-law fluids through a microtube is studied with consideration of zeta potential difference near the two-liquid interface. The modified Cauchy momentum equation in cylindrical coordinate governing the two-liquid velocity distributions is solved where both peripheral and inner liquids are represented by power-law model. The two-fluid velocity distribution under the combined interaction of power-law rheological effect and circular wall effect is evaluated at different viscosities and different electroosmotic characters of inner and peripheral power-law fluids. The velocity of inner flow is a function of the viscosities, electric properties and electroosmotic characters of two power-law fluids, while the peripheral flow is majorly influenced by the viscosity, electric property and electroosmotic characters of peripheral fluid. Irrespective of the configuration manner of power-law fluids, the shear thinning fluid is more sensitive to the change of other parameters.

A VARIATIONAL ANALYTICAL APPROACH TO TIME DEPENDENT FLOW OF OLDROYD B FLUID IN CIRCULAR TUBE

In the present investigation the time dependent flow of an Oldroyd fluid B in ahorizontal eylindrical pipe is stuided by the variational analytical approach developedby author. The tome dependent problem is mathematically reduced to a partialdifferential equation of third order. Using the improved variational approach due toKantorovich the partial differential equation can be reduced to a system of ordinarydifferential equations for different approximations. The ordinary differential equationsare solved by the method of the Laplace transform which is led to an analytical formof the solutions.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Study on the Flow Field around Two Parallel Moving Bubbles and Interaction Between Bubbles Rising in CMC Solutions by PIV

The flow fields surrounding two parallel moving bubbles rising from two identical orifices submerged in non-Newtonian fluid of carboxymethylcellulose （CMC） solution of three different mass concentration were measured experimentally by the use of particle image velocimetry （PIV）. The influences of gas flowrate, solution mass concentration, orifice interval and the angle between two bubble centers line and vertical direction on the flow field surrounding bubbles were discussed respectively by analyzing the velocity vector, velocity contours as well as individual velocity components. The results show that the liquid velocity both in front of two bubbles and behind increases with gas flowrate duo to shear-thinning effect of previous bubbles, whereas decreases with the increase of CMC concentration due to the increase of drag force acting on bubbles. The effect of the orifice interval on the flow field around two moving bubbles becomes gradually obvious as the interval becomes closer. Moreover, two adjacent side-by-side bubbles repulse each other during rising, leading to the practical interval between them increased somewhat above the orifice interval. When the distance between bubbles is less than the orifice interval 10 mm, the interaction between two neighboring bubbles changed from mutual repellence to attraction with the decrease of the angle of the line of linking two bubble centers to the vertical direction.

Fractal Analysis of Power-Law Fluid in a Single Capillary

The fractai expressions for flow rate and hydraulic conductivity for power-law fluids in a single capillary are derived based on the fractai nature of tortuous capillaries. Every parameter in the proposed expressions has clear physical meaning. The flow rate and hydraulic conductivity for power-law fluids are found to be related to the tortuosity fractal dimension and the power-law index. The flow rate for power-law fluids increases with the increasing power-law index but decreases with the increasing tortuosity fractal dimension. Good agreement between the model predictions for flow in a fractai capillary and in a converging-diverging duct is obtained. The results suggest that the fractal capillary model can be used to model the power-law fluids with different rheologicai properties.