Stabilized Finite Elements

Ph.D. Defense
Elie Hachem

KTH Mechanics, Stockholm, Sweden

Thierry Coupez
TSFPICI-
6, Ecole
June Centrale
22-24, 2009, Seoul, Korea
de Nantes

( Acknowledgments: Elie Hachem Mines-ParisTech)

1
Plan

Introduction

Équation de convection-diffusion-reaction

Formulation variationnelle et discrétisation par éléments finis

Méthodes de stabilisation : SUPG – RFB – SCPG – GGLS - EMI

Équation de Navier-Stokes incompressible

Formulation variationnelle et discrétisation par éléments finis

Méthodes de stabilisation : VMS

Couplage et applications

Conclusion

Thierry Coupez 2
Introduction: domaines d’applications

Acknowledge:
EDF - Arcelor Mittal - Snecma - Aubert & Duval - Terreal - Sfarsteel - Manoir Industries

Thierry Coupez 3
Introduction aux SFEM: Stabilised Finite Element Method

Formulation forte:
 u  f dans 

 .  ku   f dans 

u
c  .  ku   f dans  0, tn 
t

v.T  
. kT   f dans 

 T 
c 
 t
 v .T   . kT

   f dans  0, tn 

 T 
c 
 t

 v .T   . kT  r T

  f dans  0, tn 
 c est la capacité calorifiqu e
  est la masse volumique
 T la température inconnue
 k la conduc t i vité du matériau
 v la vitesse du fluide
 r coefficien t de reaction

Thierry Coupez 4
2-Stabilized finite element method for CDR equations
u
0
n
An academic test: transient convection

Unstructured coarse grid

u0 u
Skew convection a  (1,0.7) 0
n

u0

Galerkin SUPG SCPG

Profile of the solutions at t=0.5s, t=10-3

SUPG: Streamline-Upwind Petrov Galerkin

SCPG: Shock-Capturing Petrov Galerkin

Thierry Coupez 5
2-Stabilized finite element method for CDR equations
Difficultés de la résolution du problème thermique instationnaire:
Méthode de résolution numérique la plus connue et la plus classique : méthode de type Galerkin
Standard (GS) (Zienkiewicz et al.-2001, Soyris et al.-1992)

En injection (application industrielle), les grandes difficultés numériques rencontrées pour

résoudre ce problème ,proviennent du caractère spécifique des polymères : forte énergie de
dissipation, faible conductivité thermique, chocs thermiques au contact du moule et dépendance
thermique de la rhéologie

Limitations de la méthode de Galerkin :

Apparition d’instabilités (convection dominante)

Traitement du choc thermique ( diffusion pure)

1
v (1, 0) : vitessedeconvection
6
  5.10 : diffusivité
f  0 : terme source
T  100C sur 1

Thierry Coupez 6
Introduction aux SFEM: Stabilised Finite Element Method

Transient Convection-diffusion-reaction equation

Find u  V such that :

ut  L u  f in   (0, T )
u  0 on   (0, T )
u (, 0)  u0 in 
L u : a u  (u)   u
Variational formulation Vh  vh  H 01 (), vh / K islinear for K T h 
(t uh , vh )  b(uh , vh )  l (vh )  vh Vh
b(uh , vh ) : (a uh , vh )  (uh , vh )  ( uh , vh )
l (v): f , vh

Stabilized finite element method for convection dominated problem (high Peclet number)
( t uh , vh )  b(uh , vh )   ( K R uh , Lstab vh ) K  l (vh )  vh Vh
K

Where Lstab is a differential operator acting on weighting function (SUPG : Lstab  a vh )

Thierry Coupez 7
2-Stabilized finite element method for CDR equations

Formulation matricielle: U
M  CU  DU  RU  S
t
Nel
M ij    N i N j d e  matrices de masse
e
e 1
Nel
Cij     .N j N i d e  matrices de convection
e
e 1
Nel
Dij    N j .N i d e  matrices de diffusion
e
e 1
Nel
Rij    rN j N i d e  matrices de réaction
e
e 1
Nel
Si    fN i d e  vecteurs des sources
e
e 1

Forme générale:
U
M  KU  S
t

Thierry Coupez 8
2-Stabilized finite element method for CDR equations

U U U
M  KU  S  M n1 n    Sn1  KU n1   (1   )  Sn  KU n 
t t
Les schémas à deux pas de temps: θ –schémas [Raviart-Thomas 1998]

Euler explicite (θ=0)

Euler implicite (θ=1)
Crank-Nicholson (θ=1/2)

Les schémas à trois pas de temps, on trouve par exemple ceux de Richardson,
de Dupont ou de Lees. [ Soyris, 1990]

les schémas de Taylor Galerkin [Donea, 1984]

Thierry Coupez 9
2-Stabilized finite element method for CDR equations

Convection dominante : Apparitions des oscillations numériques!

Une excellente référence pour ces méthodes: Codina-1998

Le principe de base de ces méthodes est le même, il consiste à partir d’une

formulation de type GS à ajouter un terme stabilisateur au schéma.

Le principal aspect négatif de ces méthodes réside dans leur construction même. En
effet en ajoutant un terme à nos équations discrètes, nous ne résolvons plus le même
problème physique. La difficulté devient alors de pouvoir « mesurer » la pertinence du
modèle utilisé, ce qui est généralement purement empirique.
uh
( , wh )  a(uh , wh )  l ( wh )  0 wh Wh
t

 r (uh , wh ) Terme de stabilisation

nelt
où r (uh , wh )     Re (uh ) wh d e
e 1 e

Thierry Coupez 10
2-Stabilized finite element method for CDR equations

Différentes méthode:

Méthode SUPG (Streamline Upwind Petrov Galerkin). Brooks et al.-1992

Méthode ST-GLS (Space-Time/ Galerkin Least Squares). Codina-1998

Méthode SGS (Subgrid scale method). Hughes-1995

Méthode CG (Characteristic Galerkin method). Codina-1998

Méthode TG (Taylor-Galerkin method). Donea-1984

Méthode RFB (Residual Free Bubbles). Brezzi et al.-1998

Thierry Coupez 11
2-Stabilized finite element method for CDR equations

wh  wh   v .wh
Th
( , wh )  a(Th , wh )  l ( wh )  0 wh Wh
t

a(T , wh )  v .Twh d   T .wh d    v .T v .wh d 

   

 h 
l ( w )  f wh    f v .wh d 


La méthode SUPG est consistante et conservative!

Thierry Coupez 12
2-Stabilized finite element method for CDR equations

Residual Free Bubbles (RFB)

wh  wL  wB Wh  WL  WB
wB  w
K Th
K
B , wBK  B K

Trouver uh  uL  uB WL  WB tel que,


wL  WL , K  Th , et wB  B
K K


a(uL  uB , wL )  l ( wL ) •Enrichir localement l’espace éléments finis.

a(u  u K , wK )  l ( wK ) •Éliminer les inconnues bulles par condensation.

 L B B K B K
•Approche éléments finis classique

Trouver uL WL tel que,


a(u , w ) 
 L L  ˆK  (v.uL  f )(v.wL )  ( f , wL ) wL W

K
K Th

La méthode RFB est consistante et conservative!

Thierry Coupez 13
2-Stabilized finite element method for CDR equations

v hk
Pek  nombre de Peclet

hk
K  taux de stabilisation
2v

v (1, 0) : vitessedeconvection
6
  5.10 :diffusivité
f  0 : terme source
T  100C sur 1

Galerkin SUPG

Thierry Coupez 14
 2-Stabilized finite element method for CDR equations

( t uh , vh )  b(uh , vh )   ( K R uh , Lstab vh ) K  l (vh )  vh Vh

K [Link]. Comparison of some finite element
Where Lstab is a differential operator acting on weighting function methods for solving the diffusion convection-
reaction equation. Comput. Methods. in
and  K
is a local stabilizing parameter Applied Mech. and Eng., 156:185-210, 1998.

Streamline-Upwind Petrov Galerkin Shock-Capturing Petrov Galerkin

Flow
Galerkin u
a//
a┴
vh  ...   kc a// vh
i i+1 a
i-1 SUPG

 R uh
 uh if uh  0
vh  vh  k a vh a//   uh 2
2

0 if uh  0


...  S SUPG (uh , vh )   (R uh , K a vh ) K ...  S SCPG (uh , vh )   (R uh , Kc a// vh ) K

K K

...  (a u hn 1 , vh )  ...  ...  (a u hn 1 , vh ) Galerkin term

  K (a u hn 1 , a vh ) K Control along the streamline direction

K

  Kc (a// u hn 1 , a// vh ) K  ... Control along the direction of the solution gradient
K

Thierry Coupez 15
2-Stabilized finite element method for CDR equations

Usual stability coefficients

1
 4k 2 a 
K   2    
 hK hK 
 Where hK is an appropriate measure of the size of the mesh
 No time dependency

1. σ’: Temporal discretization (e.g. Euler implicit)

u hn
  1/ t  u n 1
h

, vh  (a u h , vh )  (k u h , vh )   (R uh , Lstab vh ) K  (
n 1 n 1

K t
, vh )  ( f , vh )
modified reaction term

Volker John, Ellen Schmeyer, Finite element methods for time-dependent convection–diffusion–reaction equations with small diffusion
Computer Methods in Applied Mechanics and Engineering, Volume 198, Issues 3-4, 15 December 2008, Pages 475-494

2. hK : characteristic dimension of the triangle in the streamline direction

a
1 ay
 ne ai N 
h K  2   ax
  a xi 
  h

The support length in the streamline direction

S. Micheletti, S. Perotto, M. Picasso, Some remarks on the stability coefficients and bubble stabilization of FEM on anisotropic grids, MOX
Report 06, MOX--Modeliug and Scientific Computing, Department of Mathematics "F. Brioschi", Politecnico di Milano, 2002

Thierry Coupez 16
2-Stabilized finite element method for CDR equations

Varying the time step and comparisons

0 0,2 0,4 0,6 0,8 1

0 0,2 0,4 0,6 0,8 1
1,4
1,4
1,2
1,2
1
1 t  0.01
t  0.1
0,8
0,8
0,6
0,6
0,4
0,4
t=0.5 ; x=0.75 0,2
0,2 t=0.5 ; x=0.75
0
0
-0,2
-0,2

0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1

1,4 1,4

1,2 1,2

1 1
t  0.001 t  0.0005
0,8 0,8

0,6 0,6

0,4 0,4

0,2 0,2
t=0.5 ; x=0.75 t=0.5 ; x=0.75
0 0

-0,2 -0,2

Profile of the solutions: Galerkin (solid), SUPG (dotted) and SCPG (dashed)

Thierry Coupez 17
 2-Stabilized finite element method for CDR equations

Forced convection heat transfer between two parallel plates

Mesh 1 Mesh 2 Mesh 3
T
Pe 60 - Mesh 1 0 Pe 60 - Mesh 2 Pe 60 - Mesh 3
n
1,1 1,1 1,1
(0, 0.5) (1, 0.5)
0,9 0,9 0,9

Pe 20
0,7 0,7 0,7
Temperature

Temperature
Temperature
0,5 0,5
0,5

T 1 T 0
0,3 0,3
V 0,3

0,1 0,1
0,1

-0,1 -0,1
-0,1
0 0,2 0,4 X 0,6 0,8 1 0 0,2 0,4 X 0,6 0,8 1
0 0,2 0,4 X 0,6 0,8 1

1,1 Pe 200 - Mesh 1

1,1 Pe 200 - Mesh 2 1,1 Pe 200 - Mesh 3
(0, 0) (1, 0)
0,9 T
0 0,9 0,9

0,7 n Mesh 1 Mesh 2 Mesh 3

Temperature

0,7 0,7

Pe 200

Temperature
Temperature

0,5
0,5 0,5

0,3
0,3 0,3

0,1
0,1 0,1

-0,1
-0,1 -0,1
0 0,2 0,4 X 0,6 0,8 1
0 0,2 0,4 X 0,6 0,8 1 0 0,2 0,4 X 0,6 0,8 1

1,1 Pe 2000 - Mesh 1 1,1 Pe 2000 - Mesh 2 1,1 Pe 2000 - Mesh 3

0,9 0,9 0,9

0,7 0,7 0,7

Temperature
Temperature

Temperature

Pe 2000
0,5 0,5 0,5

0,3 0,3 0,3

0,1 0,1 0,1

-0,1 -0,1 -0,1

0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 X 0,6 0,8 1 0 0,2 0,4 X 0,6 0,8 1
X

Thierry Coupez 18
2-Stabilized finite element method for CDR equations

Heat transfer at the solid level

u
Transient conduction equation: cp  .  ku   f dans  0, tn 
t
Thermal shocks occur for: t   c p h / k
2 1200

1000

Temperature °C
800
800 °C
Numerical example 600
Galerkin solution
Analytical solution
20 °C
400

200

0
0 1 2 3 4 5
Time (s)

400 °C

20 °C

Thierry Coupez 19
2-Stabilized finite element method for CDR equations

Stabilizing the transient conduction equation:

The use of enrichment method:
- The basic idea is to split the unknown u into coarse and fine scale:

~
Find u  uh  u~  V  Vh  V such that
~
b(uh , vh )  b(u , vh )  l (vh )
 v Vh
 ~)  b(u~, v~)  l (v~) ~

b (u h , v v~ V uh
u h*

u h

1. Applied to stationary diffusion equation: the gain is null!

t1 t2 t*

t
2. Discretized in time using the backward Euler method: t *

t 1 1
 (1  K )  ,  0
3. Time interpolation: t *
t 

C1hK2
K 
cp
hK2  kC2
t
T n  T n 1
 
cp
t
v d  K
K
 K k T .v d  
n

 
f v d  
K
K
( K  1) k T n 1.v d   
K
K
 K f .v d 

Thierry Coupez 20
2-Stabilized finite element method for CDR equations

An academic test: transient diffusion

1200

1000

800
EM-interpolation

Temperature °C
600 Galerkin
Analytical
400

200

0
0 1 2 3 4 5
Time (s)

Evolution of the temperature at a node

situated 2h far from the cooled side

Varying the conductivity : k= 10-4 and k=1

T°c t = 0.001s T°c t t==0.1s
0.1s T°c t t==0.5s
0.5s T°c t t==1s
1s
Exact solution
Galerkin
EM-I Exact solution
1 11 GGLS 1 Galerkin 1 Exact solution
EM-I Galerkin
GGLS EM-I
Exact
Exactsolution
solution Exact solution Exact solution Exact solution
GGLS
0,5 Galerkin
Galerkin 0,5 Galerkin 0,5 Galerkin 0,5 Galerkin
0,5
EM-I
EM-I EM-I EM-I EM-I
GGLS
GGLS GGLS GGLS GGLS
0 0 0 0
0
0 0,5 1 0 0,5 1 0 0,5 1 0 0,5 1
x 0 0,5
x 1 x x

Ilinca F. Hétu J-F, Galerkin gradient least-squares formulation for transient conduction heat transfer,
Comput. Methods Appl. Mech. Engrg. 191 (2002) 3073-3097

Thierry Coupez 21
Flow computations

Thierry Coupez 22
2-Stabilized finite element method for CDR equations

Incompressible Navier-Stokes problem

Find u and p such that :
 ( t u  u u )  2    (v)  p  f in 
v  0 in 
v  0 on 
Variational formulation
Find v and p such that where
a(v , w ) :  2  v  :   w 
a ( v , w )  b ( w , p )  c ( v, w )   f , w  w  V
b(v , q ) :     v , q 
b(v , q )  0  q  Q
c(v , w ) : (  t v , w )  (  v v , w )

Stable Mini-element (Arnold-Brezzi-Fortin) mixed formulation

(Stokes problem)

Find v  v  v and p such that

Find v and p such that
a (v , w )  b ( w , p )   f , w  w  V static condensation
a (v , w )  b ( w , p )   f , w  w  V
a(v , w )  b( w , p )   f , w   w  V 
b(v , q )  (v , q )  0  q  Q
b (v , q )  0  q  Q
Build naturally by a matrix form

Thierry Coupez 23
3-Stabilized finite element method for NS equations

Incompressible Navier-Stokes problem

 ( t u  u u )  2    (u )  p  f in  The instabilities in advection-dominated regime (?)

u  0 in 
The velocity-pressure compatibility condition (?)
u  0 on 
Stable Mini-element mixed formulation
Large scale

Find uh  V = Vh  Vh and ph  Qh such that :

  uh 

  t h     uh uh , vh    2  uh  :   vh     ph ,   vh   f , vh 
Small scale
 , v
  
  
 
 , vh    2  uh  :   vh     ph ,   vh   f , vh
 t  
V = u , u

 P1 ( K )  H 1 ( K ) , K  T , i  1,..., D 
   (uh  uh ), qh   0 h

h hK
i
i 0 i h



Matrix formulation

A 0 A u   B 
T A AT   u   B 
 vv vp 
 h  v  vv vp
 h    v 
 0 AT   u    B  A A   ph   B p 

A
bb bp 
 h  b Static condensation
 vp pp   
 A A 0   ph   B p  [Basset 06]
 vp bp    [Hachem 09] A   A A1 AT and B 
pp bp bb bp p
 AT A1B
bp bb b

Thierry Coupez 24
3-Stabilized finite element method for NS equations

Find uh Vh  Vh and ph  Qh  Q such that  (v , q)  V  Q

From a variational multiscale point of view: A  (u  u , p  p);(v, q)   ( f , v )
A  (u  u , p  p);(v , q)   ( f , v )

1- Coarse-scale equations

  (uh  uh ) 
 , vh     (uh  uh (uh  uh ), vh    ph  ph ,   vh    2  uh  :   vh     f , vh  vh  Vh
 t 
   (uh  uh ), qh   0 qh  Qh

2- Fine-scale equations

  (uh  uh ) 
 , vh     (uh  uh ) (uh  uh ), vh    ph  ph ,   vh    2  uh  :   vh     f , vh  vh  V
 t 
   (uh  uh ), qh   0 qh  Q

Assumptions to be made: p    C   uh

Quasi-Static’ subscales [R.. Codina 2002]

Approximation of the nonlinear term using only the large-scale part

only the effect of the small-scale pressure Poisson equation on the large-scale equation must be retained

Thierry Coupez 25
3-Stabilized finite element method for NS equations

Find u  uh  u and p  ph  p such that  (v, q )  V  Q

u : velocity
 B(U h , Vh )  B(U , Vh )  L(Vh )   vh , qh   Vh , Qh  p : pressure
  : fluid viscosity
 B(U h , V )  B(U , V )  L(V )   v , q   V , Q    : fluid density

where
 u 
B(U ,V ) :   , v     u u , v    p,   w    2  u  :   v       u , q 
 t 
L(V ) :  f , v 

 uh 
  , vh     uh uh , vh    2  uh  :   vh     ph ,   vh      uh , qh    f , vh 
 t 
Galerkin terms

  
K T h
K R uh ,  uh vh    
K T h
K R uh , qh 
Convection stabilization terms Pressure stabilization terms

  
K T h
C   uh ,   vh   0 vh V h , q  Q h h

grad  div stabilization term

Thierry Coupez 26
3-Stabilized finite element method for NS equations

Matrix equivalent form:

A AT AT   u   B 
 vv vb vp 
 h  v A AT   u   B 
A A AT   u    B   vv vp
 h    v 
 bv bp 
 h  b A A   ph   B p 
bb
 A A 0   ph   B p   vp pp   
 vp bp   

Two additional contributions

A. Masud, R.A. Khurram, A multiscale finite element method for the incompressible
Navier–Stokes equations , Computer Methods in Applied Mechanics and Engineering,
1) Upwind bubble functions Volume 195, (2006)

  
K T h
K R uh , uc vh  T.J.R. Hugues et al., The variational multiscale method - a paradigm for computational
mechanics , Computer Methods in Applied Mechanics and Engineering, Nov (1998)
Upwind stabilization terms
L.P. Franca, A. Nesliturk, On a two-level finite element method for the incompressible
2) Modelling the fine-scale pressure Navier–Stokes equations, Int. J. Numer. Methods Engrg. (2001)

  
K T h
C   uh ,   vh  R. Codina, Stabilized finite element approximation of transient incompressible flows using
orthogonal subscales, Comput. Methods Appl. Mech. Engrg. (2002)
grad  div stabilization term

Thierry Coupez 27
2-Stabilized finite element method for NS equations

2D Lid driven cavity:

High Reynolds number

Comparisons with the reference
boundary conditions (left) , coarse mesh 64x64(center) , 180x180 fine mesh (right)

0,6 0,6 0,6 0,6

Re=1000 Re=5000 Re=10000 Re=20000
0,4 0,4 0,4
0,4
0,2 0,2
0,2 0,2

uy

uy
uy

0
uy

0
0 0 0,5 1 0
0 0,5 1
-0,2 -0,2 0 0,5 1
0 0,5 1
-0,2 -0,2
-0,4 -0,4

-0,4 -0,6 -0,4

-0,6
x x x x
-0,6 -0,8 -0,8 -0,6

0,6 0,6 0,6 0,6

Re=1000 Re=5000 Re=10000 Re=20000
0,4 0,4 0,4
0,4
0,2 0,2 0,2
0,2
uy
uy

uy
0 0 0
uy

0 0 0,5 1 0 0,5 1 0 0,5 1

-0,2 -0,2 -0,2
0 0,5 1
-0,2 -0,4 -0,4
-0,4
-0,4 -0,6 -0,6 -0,6
x x x x
-0,6 -0,8 -0,8 -0,8

Velocity profile for uy along y=0.5

____ coarse mesh ___ fine mesh ◊ reference (601x601)
Thierry Coupez 28
3-Stabilized finite element method for NS equations

Test case: driven cavity

- Increasing number of vortices
- Apparition of new resolved vortices

14
13
12
Number of vortices

10
10
8
8
6
6
4 5
2 3
Re
0
0 10000 20000 30000 40000 50000

0,6

0,4 Re=33000
Re=50000
0,2 Streamline on colored velocity distribution from top-left to bottom-right: Re=1000, 5000, 10000, 20000, 33000 and 50000 with 180×180 mesh
uy

0
0 0,5 1
-0,2

-0,4
x
-0,6

Velocity profile for uy along y=0.5

Thierry Coupez 29
3-Stabilized finite element method for NS equations

Test case: BFS (backward facing step)

Flow computations using both isotropic
and anisotropic meshes

- High Reynolds number (44000)

- Anisotropic mesh adaptation (~19000 nodes)

E. Hachem, T. Kloczko, H. Digonnet and T. Coupez,, ” Stable mixed-finite element

method for incompressible flows with high Reynolds number” preprint for INJME

Thierry Coupez 30
3-Stabilized finite element method for NS equations

3D Lid driven cavity:

Re=1000, 3200 and 12000

Re = 1000

Re = 3200

Re = 12000

Thierry Coupez 31
3-Stabilized finite element method for NS equations

3D Lid driven cavity:

Re=1000, 3200 and 12000

36, 282 nodes and 192, 080 238, 580 nodes and 1, 229, 089
linear tetrahedral elements linear tetrahedral elements

Re = 1000

Re = 3200

Re = 12000

Vortex formation and characteristic of the flow

Velocity profiles for ux and uz along the centerline [Y. Zang 1993] [R. Bouffanais 2007]
Thierry Coupez 32
4-Coupled Heat Transfer and Fluid Flow

Thierry Coupez 33
4-Coupled Heat Transfer and Fluid Flow

Thierry Coupez 34
2-Stabilized finite element method for CDR equations

3m de rayon Pièce à chauffer

Coupe du maillage
Brûleurs à 1000°C

Thierry Coupez 35
5-Radiative Transport Equation (RTE)

The P1 model :

Very simple to implement

Diffusion-reaction equation
Volumetric source term

Diffusive model
Does not account for better directional influence
(see Thibauld’s work)

Hat-shaped disk

G :incident radiation
k :absorbtion coefficient

Thierry Coupez 36
5-Radiative Transport Equation (RTE)

Comparisons with experimental data

Thierry Coupez 37
4-Coupled Heat Transfer and Fluid Flow

Thierry Coupez 38
6- Applications

Industrial furnace: Terreal

-Inlet velocity : 75m/s

-Temparature: 1100°C

-Dimension: 1.5x1.3x1.0

-Reynolds ~ 106

-32 2.4 Ghz Opteron cores in parallel

-5 days of computations

-155015 nodes and 896539 linear tetrahedral elements

E. Hachem, H. Digonnet, T. Coupez*, ”Stabilized

finite element method for heat transfer and fluid flow
inside industrial furnaces”
MS17-2 Symposium Celebrating the 65th Birthday of
T.J.R. Hugghes, Friday, 15:00,, Room A

Thierry Coupez 39
6- Applications

Industrial furnace: mesh and setup

Thierry Coupez 40
6- Flow characteristics

Thierry Coupez 41
6- Applications

Snapshots: flow characteristics

Thierry Coupez 42
6- Applications

Results and comparisons

E. Hachem, E. Massoni, H. Lemercier and T. Coupez,, ” Numerical simulation in a full-

scale industrial furnace” preprint for INJME

Thierry Coupez 43
6- Applications

Industrial furnace:

-10 burners (average velocity 40m/s)

-Temparature: 1000°C

-Dimension:6m diameter x 2.7m height

-Reynolds ~ 106

-32 2.4 Ghz Opteron cores in parallel

-~800000 linear tetrahedral elements

Thierry Coupez 44
6- Applications

Industrial furnace:

E. Hachem, T. Kloczko, H. Digonnet and T. Coupez, ” Finite element solution to

handle complex heat and fluid flows in industrial furnaces” preprint for INJMF

Thierry Coupez 45
6- Applications

Industrial furnace:

E. Hachem, H. Digonnet, T. Coupez****, ”Stabilized finite element method for heat transfer and
fluid flow inside industrial furnaces” FEF 2009, Tokyo, Japan
MS17-2 Symposium Celebrating the 65th Birthday of T.J.R. Hugghes

Thierry Coupez 46
6- Perspective and future works

More test…for higher accuracy!

Convection-Conduction
IMV: immersed volume method

Thierry Coupez 47
6- Perspective and future works

Results and comparisons

Thierry Coupez 48
6- Perspective and future works

Air-Cooling: Titanium (Snecma)

Thierry Coupez 49
6- Perspective and future works

Results and Comparisons

900 900
900
800 800
800
700 700
700
600 600
600
500 TC1 500 TC2
500 TC3
400 T1 400 T2
400 T3
300 300
300
200 200
200
100 100
100
0 0
0
0 100 200 300 400 500 600 700 0 100 200 300 400 500 600 700
0 100 200 300 400 500 600 7

900 900 900

800 800 800
700 700 700
600 600 600

500 TC5 500 TC6 500 TC7

400 T5 400 T6 400 T7

300 300 300

200 200 200

100 100 100

0 0 0
0 100 200 300 400 500 600 700 0 100 200 300 400 500 600 700 0 100 200 300 400 500 600 7

Thierry Coupez 50
 Experimental validation :

G. François, E. Hachem, T. Coupez, ” Méthodes

éléments finis pour le remplissage à haut nombre
de Reynolds” CFM, Marseille, 2009

Thierry Coupez 51
 Mixing , twin screw extrusion
E. Foudrinier et R. Valette :

Representation of immersed domain (moving screws)

by a Level Set approach

To compute complex flows

inside extruders (free
surfaces, velocities, strain
rates, pressures, filling, …)

Immersion
Thierry Coupez 5252
 Zienkiewicz O.C. et Taylor R.-L. 1989. Finite Element Method: The Basis. Volume 1. Mc Graw-Hill Book Compagny
Soyris N., Tronel Y. et Chenot J. 1992. Finite element calculation of thrmo-coupled large deformation in hot forging. Conduction, Radiation and Phase
Change, 493-511.
Aliaga C. 2000. Simulation numérique par éléments finis en 3D du comportement thermomécanique au cours du traitement thermique d'aciers :
application à la trempe de pièces forgées ou coulées. Thèse de doctorat, Ecole des Mines de paris, 2000.
Brooks A.N. et Hughes T.J.R. 1982. Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the
incompressible Navier-Stokes equations. Comput. Methods Appl. Mech, 32, 199-259.
Magnin B. 1994. Modélisation du remplissage des moules d’injection pour polymères thermoplastiques par une méthode eulérienne-lagrangienne
arbitraire. Thèse de Doctorat de l’Ecole Nationale Supérieure des Mines de Paris.
Pittman J.F.T. et Mahmoudzadeh H. 1983. Performance of melt thermocouples. Part 2: The upstream-pointing parallel-to-flow type. Polymer Engineering
Review, 3, 75.
Ghoneim H. 1987. Simulation of balancing multi-cavity molds. International Polymer Processing 1, 166.
Vila Real P.M.M. et Magalhaes Oliveira C.A. 1989. Use of polynomial functions for hierarchical improvement of the finite element solution for transient
heat conduction problems with high thermal gradients. Numerical Methods in Thermal Problems, 6, 227-240.
Gaston L. 1997. Simulation numérique par éléments finis bidimensionnels du remplissage de moules de fonderie et étude expérimentale sur maquette
hydraulique. Thèse de Doctorat de l’Ecole Nationale Supérieure des Mines de Paris.
Codina R. 1998. Comparison of some finite element methods for solving the diffusion-convection-reaction equation. Comput. Methods Appl. Mech, 156,
185-210.
Hughes T.J.R. 1995. Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of
stabilized methods. Comput. Methods Appl. Mech, 127, 387-401.
Donea J. 1984. A Taylor-Galerkin method for convection transport problems. Int. J. Numer. Methods Engrg., 20, 101-119.
Brezzi F., Marini D. et Russo A. [Link] of the pseudo residual free bubblesto the stabilizationof the convection-diffusion problems. Comput.
Methods Appl. Mech, 166, 3-24.
Cockburn B., Karniadakis G. et Shu C.-W. 2000. The development of discontinuous Galerkin methods. Lecture note in Computational Science and
engineering, 11, Springer.
Ern A. et Guermond J.-L. 2004. Theory and practice of finite elements. Applied Mathematical Sciences, 159, Springer.
Shu C.-W. et Cockburn B. 2001. Runge-Kutta discontinuous Galerkin methods for convection-dominated problems. Journal of scientific computing, 16, 3,
173-261.
Bernacki M. 2005. Méthodes de type Galerkin discontinu pour la propagation des ondes en aéroacoustique. Thèse de Doctorat de l’Ecole Nationale des
Ponts et Chaussés.
Bernacki M., Lanteri S. et Piperno S. 2006. Time-domain parallel simulation of heterogeneous wave propagation on unstructured grids using explicit, non
diffusive, discontinous Galerkin methods. Journal of Computational Acoustics, Vol. 14, No. 1 (2006) 57-81, 2006.
Pichelin E. et Coupez T. 1999. A Taylor discontinuous Galerkin method for the thermal solution in 3D mold filling. Comput. Methods Appl. Mech, 178, 153-
169.
Batkam S., Bruchon J. et Coupez T. 2002. A space-time discontinuous Galerkin method for convection and diffusion in injection moulding. Inter. J. of
Forming Processes.
Silva, L., Gruau, C., Agassant, J.-F., Coupez. T. et Mauffrey, J. 2005. Advanced Finite Element Model for 3D Injection Molding. International Polymer
Processing. Vol.20, n°3, pp. 265-273

Thierry Coupez 53