Turbulence Modulation by Slender Fibers

In this paper, we numerically investigate the turbulence modulation produced by long flexible fibres in channel flow. The simulations are based on an Euler–Lagrangian approach, where fibres are modelled as chains of constrained, sub-Kolmogorov rods. A novel algorithm is deployed to make the resoluti...

Full description

Bibliographic Details
Main Authors: Davide Di Giusto, Cristian Marchioli
Format: Article
Language:English
Published: MDPI AG 2022-07-01
Series:Fluids
Subjects:
Online Access:https://www.mdpi.com/2311-5521/7/8/255
_version_ 1827628233308241920
author Davide Di Giusto
Cristian Marchioli
author_facet Davide Di Giusto
Cristian Marchioli
author_sort Davide Di Giusto
collection DOAJ
description In this paper, we numerically investigate the turbulence modulation produced by long flexible fibres in channel flow. The simulations are based on an Euler–Lagrangian approach, where fibres are modelled as chains of constrained, sub-Kolmogorov rods. A novel algorithm is deployed to make the resolution of dispersed systems of constraint equations, which represent the fibres, compatible with a state-of-the-art, Graphics Processing Units-accelerated flow-solver for direct numerical simulations in the two-way coupling regime on High Performance Computing architectures. Two-way coupling is accounted for using the Exact Regularized Point Particle method, which allows to calculate the disturbance generated by the fibers on the flow considering progressively refined grids, down to a quasi-viscous length-scale. The bending stiffness of the fibers is also modelled, while collisions are neglected. Results of fluid velocity statistics for friction Reynolds number of the flow <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>R</mi><msub><mi>e</mi><mi>τ</mi></msub><mo>=</mo><mn>150</mn></mrow></semantics></math></inline-formula> and fibers with Stokes number <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>S</mi><mi>t</mi></mrow></semantics></math></inline-formula> = 0.01 (nearly tracers) and 10 (inertial) are presented, with special regard to turbulence modulation and its dependence on fiber inertia and volume fraction (equal to <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>ϕ</mi><mo>=</mo><mn>2.12</mn></mrow></semantics></math></inline-formula>·<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msup><mn>10</mn><mrow><mo>−</mo><mn>5</mn></mrow></msup></semantics></math></inline-formula> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mn>2.12</mn></mrow></semantics></math></inline-formula>·<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msup><mn>10</mn><mrow><mo>−</mo><mn>4</mn></mrow></msup></semantics></math></inline-formula>). The non-Newtonian stresses determined by the carried phase are also displayed, determined by long and slender fibers with fixed aspect ratio <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mi>λ</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>200</mn></mrow></semantics></math></inline-formula>, which extend up to the inertial range of the turbulent flow.
first_indexed 2024-03-09T13:27:49Z
format Article
id doaj.art-fed22a819f2b4da8a73649b19cbfcf9c
institution Directory Open Access Journal
issn 2311-5521
language English
last_indexed 2024-03-09T13:27:49Z
publishDate 2022-07-01
publisher MDPI AG
record_format Article
series Fluids
spelling doaj.art-fed22a819f2b4da8a73649b19cbfcf9c2023-11-30T21:22:40ZengMDPI AGFluids2311-55212022-07-017825510.3390/fluids7080255Turbulence Modulation by Slender FibersDavide Di Giusto0Cristian Marchioli1Department of Engineering and Architecture, University of Udine, 33100 Udine, ItalyDepartment of Engineering and Architecture, University of Udine, 33100 Udine, ItalyIn this paper, we numerically investigate the turbulence modulation produced by long flexible fibres in channel flow. The simulations are based on an Euler–Lagrangian approach, where fibres are modelled as chains of constrained, sub-Kolmogorov rods. A novel algorithm is deployed to make the resolution of dispersed systems of constraint equations, which represent the fibres, compatible with a state-of-the-art, Graphics Processing Units-accelerated flow-solver for direct numerical simulations in the two-way coupling regime on High Performance Computing architectures. Two-way coupling is accounted for using the Exact Regularized Point Particle method, which allows to calculate the disturbance generated by the fibers on the flow considering progressively refined grids, down to a quasi-viscous length-scale. The bending stiffness of the fibers is also modelled, while collisions are neglected. Results of fluid velocity statistics for friction Reynolds number of the flow <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>R</mi><msub><mi>e</mi><mi>τ</mi></msub><mo>=</mo><mn>150</mn></mrow></semantics></math></inline-formula> and fibers with Stokes number <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>S</mi><mi>t</mi></mrow></semantics></math></inline-formula> = 0.01 (nearly tracers) and 10 (inertial) are presented, with special regard to turbulence modulation and its dependence on fiber inertia and volume fraction (equal to <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>ϕ</mi><mo>=</mo><mn>2.12</mn></mrow></semantics></math></inline-formula>·<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msup><mn>10</mn><mrow><mo>−</mo><mn>5</mn></mrow></msup></semantics></math></inline-formula> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mn>2.12</mn></mrow></semantics></math></inline-formula>·<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msup><mn>10</mn><mrow><mo>−</mo><mn>4</mn></mrow></msup></semantics></math></inline-formula>). The non-Newtonian stresses determined by the carried phase are also displayed, determined by long and slender fibers with fixed aspect ratio <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mi>λ</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>200</mn></mrow></semantics></math></inline-formula>, which extend up to the inertial range of the turbulent flow.https://www.mdpi.com/2311-5521/7/8/255multi-phase flowsfibersdrag reduction
spellingShingle Davide Di Giusto
Cristian Marchioli
Turbulence Modulation by Slender Fibers
Fluids
multi-phase flows
fibers
drag reduction
title Turbulence Modulation by Slender Fibers
title_full Turbulence Modulation by Slender Fibers
title_fullStr Turbulence Modulation by Slender Fibers
title_full_unstemmed Turbulence Modulation by Slender Fibers
title_short Turbulence Modulation by Slender Fibers
title_sort turbulence modulation by slender fibers
topic multi-phase flows
fibers
drag reduction
url https://www.mdpi.com/2311-5521/7/8/255
work_keys_str_mv AT davidedigiusto turbulencemodulationbyslenderfibers
AT cristianmarchioli turbulencemodulationbyslenderfibers