On Numerical Simulations of Complex Flows of Viscoplastic Materials

2000 ◽  
Author(s):  
Paulo R. Souza Mendes ◽  
Mônica F. Naccache ◽  
Harry T. M. Vinagre

Abstract The performance of a typical numerical simulation for complex flows of viscoplastic materials was examined. The inertialess flow of viscoplastic materials through an axisymmetric channel formed by an abrupt expansion followed by a contraction was employed with this purpose. Flow visualization experiments were performed with a well characterized Carbopol aqueous solution. Numerical solutions of the mass and momentum balance equations were obtained, using the Generalized Newtonian Liquid model with a biviscosity function. The flow visualization results showed that the flow pattern is essentially Newtonian for large expansion lengths. For smaller expansion lengths, however, flow is observed only in an inner axisymmetric region whose diameter is approximately the same as the one of the inlet and outlet tubes. Outside this region the flow is stagnant, and a slip interface between these two regions seems to occur. The corresponding numerical solution was not capable of predicting the observed flow pattern.

2001 ◽  
Author(s):  
Erick F. Quintella ◽  
Paulo R. Souza Mendes ◽  
Luis F. A. Azevedo ◽  
Mônica F. Naccache

Abstract Stenosis is a constriction or narrowing of a duct or passage. This word is often used to refer to constrictions that occur in blood vessels, especially arteries. A stenotic vessel generally causes blood flow stagnation and hence formation of thrombus, which may be rather harmful to the human physiology. In this work we perform flow visualization experiments and finite-volume computation to study the flow of an yield-stress (or highly shear-thinning) liquid through a tube in the neighborhood of an axisymmetric constriction. The rheology of the liquid was chosen to be close to the one of blood with 55% hematocrit Reynolds and Yield number values characteristic to different arteries have been studied, as well as different levels of constriction. Results for a Newtonian liquid were also obtained, for comparison purposes. Numerical solutions of the mass and momentum balance equations were obtained. In these solutions it was assumed that the materials behave like a Generalized Newtonian Liquid with a biviscosity function, which mimics closely the Herschel-Bulkley equation but allows deformation below the yield stress limit. Among other findings, it has been observed that the non-Newtonian rheology causes a significant change in the flow pattern in the neighborhood of the stenosis. Because the conditions for thrombus formation are directly related to the flow pattern, studies of the flow through stenotic vessels that assume a Newtonian rheology may lead to erroneous conclusions.


Introduction .—In nearly all the previous determinations of the ratio of the specific heats of gases, from measurements of the pressures and temperature before and after an adiabatic expansion, large expansion chambers of fror 50 to 130 litres capacity have been used. Professor Callendar first suggests the use of smaller vessels, and in 1914, Mercer (‘Proc. Phys. Soc.,’ vol. 26 p. 155) made some measurements with several gases, but at room temperature only, using volumes of about 300 and 2000 c. c. respectively. He obtained values which indicated that small vessels could be used, and that, with proper corrections, a considerable degree of accuracy might be obtained. The one other experimenter who has used a small expansion chamber, capacity about 1 litre, is M. C. Shields (‘Phys. Rev.,’ 1917), who measured this ratio for air and for hydrogen at room temperature, about 18° C., and its value for hydroger at — 190° C. The chief advantage gained by the use of large expansion chambers is that no correction, or at the most, a very small one, has to be made for any systematic error due to the size of the containing vessels, but it is clear that, in the determinations of the ratio of the specific heats of gases at low temperatures, the use of small vessels becomes a practical necessity in order that uniform and steady temperature conditions may be obtained. Owing, however, to the presence of a systematic error depending upon the dimensions of the expansion chamber, the magnitude of which had not been definitely settled by experiment, the following work was undertaken with the object of investigating the method more fully, especially with regard to it? applicability to the determination of this ratio at low temperatures.


2012 ◽  
Vol 26 (29) ◽  
pp. 1250146 ◽  
Author(s):  
BHASKAR SEN GUPTA ◽  
SHANKAR P. DAS

The renormalized dynamics described by the equations of nonlinear fluctuating hydrodynamics (NFH) treated at one loop order gives rise to the basic model of the mode coupling theory (MCT). We investigate here by analyzing the density correlation function, a crucial prediction of ideal MCT, namely the validity of the multi step relaxation scenario. The equilibrium density correlation function is calculated here from the direct solutions of NFH equations for a hard sphere system. We make first detailed investigation for the robustness of the correlation functions obtained from the numerical solutions by varying the size of the grid. For an optimum choice of grid size we analyze the decay of the density correlation function to identify the multi-step relaxation process. Weak signatures of two step power law relaxation is seen with exponents which do not match predictions from the one loop MCT. For the final relaxation stretched exponential (KWW) behavior is seen and the relaxation time grows with increase of density. But apparent power law divergences indicate a critical packing fraction much higher than the corresponding MCT predictions for a hard sphere fluid.


ASAIO Journal ◽  
1998 ◽  
Vol 44 (5) ◽  
pp. M714-M718 ◽  
Author(s):  
Masako Nakata ◽  
Toru Masuzawa ◽  
Eisuke Tatsumi ◽  
Yoshiyuki Taenaka ◽  
Takashi Nishimura ◽  
...  

1999 ◽  
Author(s):  
Alexander V. Kasharin ◽  
Jens O. M. Karlsson

Abstract The process of diffusion-limited cell dehydration is modeled for a planar system by writing the one-dimensional diffusion-equation for a cell with moving, semipermeable boundaries. For the simplifying case of isothermal dehydration with constant diffusivity, an approximate analytical solution is obtained by linearizing the governing partial differential equations. The general problem must be solved numerically. The Forward Time Center Space (FTCS) and Crank-Nicholson differencing schemes are implemented, and evaluated by comparison with the analytical solution. Putative stability criteria for the two algorithms are proposed based on numerical experiments, and the Crank-Nicholson method is shown to be accurate for a mesh with as few as six nodes.


2018 ◽  
Vol 844 ◽  
pp. 61-91 ◽  
Author(s):  
Weihua Li ◽  
Satish Kumar

The coating of discrete objects is an important but poorly understood step in the manufacturing of a broad variety of products. An important model problem is the flow of a thin liquid film on a rotating cylinder, where instabilities can arise and compromise coating uniformity. In this work, we use lubrication theory and flow visualization experiments to study the influence of surfactant on these flows. Two coupled evolution equations describing the variation of film thickness and concentration of insoluble surfactant as a function of time, the angular coordinate and the axial coordinate are solved numerically. The results show that surface-tension forces arising from both axial and angular variations in the angular curvature drive flows in the axial direction that tend to smooth out free-surface perturbations and lead to a stable speed window in which axial perturbations do not grow. The presence of surfactant leads to Marangoni stresses that can cause the stable speed window to disappear by driving flow that opposes the stabilizing flow. In addition, Marangoni stresses tend to reduce the spacing between droplets that form at low rotation rates, and reduce the growth rate of rings that form at high rotation rates. Flow visualization experiments yield observations that are qualitatively consistent with predictions from linear stability analysis and the simulation results. The visualizations also indicate that surfactants tend to suppress dripping, slow the development of free-surface perturbations, and reduce the shifting and merging of rings and droplets, allowing more time for solidifying coatings in practical applications.


Materials ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1099 ◽  
Author(s):  
Hanxue Cao ◽  
Chao Shen ◽  
Chengcheng Wang ◽  
Hui Xu ◽  
Juanjuan Zhu

Although numerical simulation accuracy makes progress rapidly, it is in an insufficient phase because of complicated phenomena of the filling process and difficulty of experimental verification in high pressure die casting (HPDC), especially in thin-wall complex die-castings. Therefore, in this paper, a flow visualization experiment is conducted, and the porosity at different locations is predicted under three different fast shot velocities. The differences in flow pattern between the actual filling process and the numerical simulation are compared. It shows that the flow visualization experiment can directly observe the actual and real-time filling process and could be an effective experimental verification method for the accuracy of the flow simulation model in HPDC. Moreover, significant differences start to appear in the flow pattern between the actual experiment and the Anycasting solution after the fragment or atomization formation. Finally, the fast shot velocity would determine the position at which the back flow meets the incoming flow. The junction of two streams of fluid would create more porosity than the other location. There is a transition in flow patterns due to drag crisis under high fast shot velocity around two staggered cylinders, which resulted in the porosity relationship also changing from R1 < R3 < R2 (0.88 m/s) to R1 < R2 < R3 (1.59 and 2.34 m/s).


Sign in / Sign up

Export Citation Format

Share Document