1 / 21

Computational Rheology Isaac Newton Institute Dynamics of Complex Fluids -10 Years on

Computational Rheology Isaac Newton Institute Dynamics of Complex Fluids -10 Years on. Mike Webster. Schlumberger, UNAM-(Mexico), INNFM. Juan P. Aguayo Hamid Tamaddon Mike Webster. Institute of non-Newtonian Fluid Mechanics EPSRC Portfolio Partnership.

adonica
Download Presentation

Computational Rheology Isaac Newton Institute Dynamics of Complex Fluids -10 Years on

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Computational Rheology Isaac Newton InstituteDynamics of Complex Fluids -10 Years on Mike Webster Schlumberger, UNAM-(Mexico), INNFM Juan P. Aguayo Hamid Tamaddon Mike Webster Institute of non-Newtonian Fluid Mechanics EPSRC Portfolio Partnership

  2. To achieve highly elastic, high strain-rate/deformation rate solutions (polymer melts & polymer solutions) To quantitatively predict pressure-drop, as well as flow field structures (vortices, stress distributions) To accurately represent transient flow evolution in complex flows To quantitatively predict multiple-scale response (multi-mode) To achieve compressible viscoelastic representations Computational Rheology – Some Outstanding Challenges

  3. TRANSIENT & STEADYContraction Flows EPTT Oldroyd Planar Axisymmetric Planar Axisymmetric

  4. Pressure-drop vs flow-rate in contractions axisymmetric planar Fluid viscosity = 1.75Pa.s – 8:1 contraction, exit length 7.4mm Fluid viscosity = 1.75Pa.s – 20:1 contraction, exit length 40mm  Newtonian syrup  Boger fluid

  5. Pressure drop (epd) vs. We, 4:1:4 axisymmetric Szabo et al. with FENE-CR J. Non-Newt. Fluid Mech. 72:73-86, 1997 epd We

  6. Schematic diagram for a) 4:1:4 contraction/expansion, b) 4:1 contraction Szabo et al. J. Non-Newt. Fluid Mech. 72:73-86, 1997 Rothstein and McKinley J. Non-Newt. Fluid Mech. 86:61-88, 1999 J. Non-Newt. Fluid Mech. 98:33-63, 2001 Wapperom and Keunings J. Non-Newt. Fluid Mech. 97:267-281, 2001 Excess pressure drop (epd - P ) : Total pressure drop

  7. Pressure-drop (epd) vs. We, Oldroyd-B, a, c) axisymmetric, b, d) planar 4:1:4 4:1 a) c) a) Axisymmetric b) b) d) Planar

  8. Pressure profile around constriction zone, 4:1:4 axisymmetric and planar case Oldroyd-B, =0.9

  9. N1p 3D view – 4:1:4 contraction/expansion Axisymmetric Planar Oldroyd-B, =0.9

  10. (P - PNewt) and stress profiles along wall, 4:1 and 4:1:4 axisymmetric case 4:1 4:1:4 Oldroyd-B, =0.9

  11. (P - PNewt)and stress profiles along wall – 4:1:4 planar and axisymmetric case Planar Axisymmetric Oldroyd-B, =0.9

  12. Pressure-drop (epd) vs. We, 4:1:4 axisymmetric, alternative models epd We epd We

  13. upturn & enhanced epd upturn epd monotonic decrease epd (P - PNewt) profiles along wall – 4:1:4 axisymmetric, increasing 

  14. Alternative differential pressure-drop measure Since &by calibration 

  15. 0 Rate of dissipation & pressure-drop, 4:1:4 definition rate of dissipation , Seeking {P – 1} > 0

  16. Pressure-drop (epd) vs. a)We,b) upstream sampling distance, 4:1:4 axisymmetric epd We epd

  17. upturn & enhanced epd mono-dec epd 4:1:4 axisymmetric vortex cell size, Oldroyd-B, change =0.99 =1/9

  18. Rheological properties: Oldroyd-B, LPTT, EPTT, SXPP a) Oldroyd-B extensional viscosity, b b) Shear and extensional viscosity, b=0.9 c) Shear and extensional viscosity, b=0.99

  19. NEW BOGER fluid modelling & Pressure Drop Axisymmetric contraction Planar contraction

  20. Centreline pressure gradient 4:1:4 axisymmetric, Oldroyd-B =1/9 =0.9 =0.99

More Related