head loss estimation in diaphragm type configuration at ... · pressure drop profiles in unsteady...

1
A.GRAND-CLÉMENT Contextualization and Mesh configurations to optimize computational time Physical phenomena Comparison results for different mesh generators with Code_Saturne Difficulties encountered Similar results compared with previous computational calculations 4% deviation with experimental results for the tetrahedral mesh with fine discretization at the diaphragm surface and boundary layers Pressure drop profiles in unsteady calculation Comparison between Fluent and Code_Saturne (present study) results Head loss estimation in diaphragm type configuration at EDF hydraulic engineering Block division for the hexa mesh generation Flow recirculation after the bent Small unsteadyness after diaphragm Concerns : The surge chamber allows to decrease the high pressure level in the inlet tunnel and penstock. The high pressure drop in the diaphragm reduces the water level variation in the surge shaft The sizing of this device is an important part to control water level and to avoid overflowing in the valley. K hexa mesh K tetra mesh Salomé K=0.0375 K=0.028 Element number 56889 1226000 Ansys K=0.0345 K=0.025 Element number 40000 231823 Experimental K=0.0240 Lowest deviation (%) 4% 10 15 20 25 30 52 54 56 58 60 62 64 temps (s) pertes de charge en (mce) Fluent 6.3.26 Code_Saturne 84m 84m y+ too high in the interest domain even for k-eps model Difficulties to handle the meshes defaults Divergence solver for boundary layers with hexaedral mesh Method for extracting value on a surface Steady algorithm Unsteady algorithm K [m-5.s-2] Head losses coefficient 0.027 0.028 Surge chamber geometry Horseshoe-shaped diaphragm mesh Scalar diffusion to control water level Distance to the boundary 84m 84m t=13s t=26s t=27s Comparison between steady and unsteady algorithm Tetra mesh with refinement and boundary layers Hexa mesh without refinement and boundary layers Diaphragm V=cste Outlet p=0 Surge shaft Inlet tunnel t=28s t=13s t=13.5s t=14s t=15s

Upload: others

Post on 16-Mar-2020

12 views

Category:

Documents


0 download

TRANSCRIPT

A.GRAND-CLÉMENT

Contextualization and Mesh configurations to optimize computational time

Physical phenomena

Comparison results for different mesh generators with Code_Saturne

Difficulties encountered

Similar results compared withprevious computational calculations

4% deviation with experimental results for the tetrahedral mesh with fine discretization at the

diaphragm surface and boundary layers

Pressure drop profiles in unsteady calculationComparison between Fluent and Code_Saturne

(present study) results

Head loss estimation in diaphragm type configuration at EDF hydraulic engineering

• Block division for the hexa mesh generation

Flow recirculation after the bentSmall unsteadyness after diaphragm

Concerns : The surge chamberallows to decrease the high pressure level in the inlet tunnel and penstock.

The high pressure drop in the diaphragm reduces the water level

variation in the surge shaft

The sizing of this device is an important part to control water level

and to avoid overflowing in the valley.

Khexa mesh

K tetra mesh

Salomé K=0.0375 K=0.028

Element number 56889 1226000

Ansys K=0.0345 K=0.025

Element number 40000 231823

Experimental K=0.0240

Lowest deviation(%)

4%

10 15 20 25 30

52

54

56

58

60

62

64

temps (s)

pe

rte

s d

e c

har

gee

n (

mce

) Fluent 6.3.26

Code_Saturne

84m

84m

• y+ too high in the interest domaineven for k-eps model

Difficulties to handle the meshesdefaults

Divergence solver for boundary layerswith hexaedral mesh

Method for extracting value on a surface

Steadyalgorithm

Unsteadyalgorithm

K [m-5.s-2]Head losses coefficient

0.027 0.028

Surge chamber geometry

Horseshoe-shaped diaphragm mesh

Scalar diffusion to control water level

Distance to the boundary 84m

84m

t=13s

t=26s

t=27s Comparison between steady and unsteadyalgorithm

Tetra mesh withrefinement and boundary layers

Hexa meshwithout

refinement and boundary layers

Diaphragm

V=cste

Outletp=0

Surge shaft

Inlet tunnel

t=28s

t=13s

t=13.5s

t=14s

t=15s