Mixed Hybid Finite Element nd Itetive Metods fo Flo in Poous Medi E. Moue, C. Le Potie, P. Mugis,, L.V. Benet. Commissit à l'enegie Atomique, C.E. de Sly, Gif su Yvette Cedex, Fne Emil: : emmnuel.moue@e.f
Summy Cst3M Code Rid s eqution & MHFE Fomultion Itetive esolution Illusttions Ai Wte migtion + Tempetue Conlusion
Cst3M (CASTEM2000 ttp: //-st3m.e.f Finite Element Code (CEA PDE : Stutul menis, Fluid Menis, Temis, Objet oiented ode : Objet2 = Opeto Objet (options «Toll Box» (500 opetos 2 lnguges : use s (Gibine & develope s (Esope Pe & Post poessos FE, MHFE & FV
Cst3M (CASTEM2000 Poous medi : Dy (2d & 3d, Tnspot equtions (dv. disp. diff., (See eent ppe in omputtionl geosienes 2004, Gilles Bend Miel et l, bout te COUPLEX test se Bsi bik fo multipse flo : F A + BF +.J J = C F A, B,C ons tn ts = S(x, t Itetive Solutions Fo multipse flo F, G,
Rid s Eqution, Itetive Resolution Mixed Hybid Finite Element fomultion fo te flo in unstuted poous medi», C.Le Potie et l., CMWR XII, 998 C( =.U U = K( θ( + z MHFE : ed pessue, U : Dy veloity, C : Cpilly pity, K : pemebility, _ te ontent FE Mes (QUA4 K Ω Ω U.dΩ= C( Ω vdω= MHFE Mes EFMH (QUAF.dΩ Ω Ω.UvdΩ.ndΓ Element ente, ed Fe ente Veloity, ed
Rid s Eqution, Itetive Resolution Impliit time disetiztion & Pid lgoitm Time step : n Itetive step : i H = + z nd C( U,i,i+ H,i = K( + Δt,i n H =.U,i+ H,i+
Rid s Eqution, Itetive Resolution Time step sttegy : Δt = U Fe.n Fe sufe( Ω N Fes N volume( Ω k F k Pmete (X= C & K omogeneiztion : Diffeent types of mens : funtions of X vlues on te fes, X=F(X(Hfe, XA (itmeti, XG (geometi, XH (moni o enteed vlue X = X(Hente Element ente; X = X(Hente MHFE Mes MHFE (QUAF Fe ente; X =F(X(Hfe
Infilttion in Heteogeneous Soil Infilttion Infilttion nd Rege of n quife in eteogeneous soil 30mn Wte ontent Snd Cly lens Not t Sle dy Oveflo
Stution Rinfll on Slope (Cf DYNAS Rinfll on slope (Runoff, Infilttion, Rege, Oveflo, Lolized Rinfll Unifom Rinfll T0 Aquife T Runoff
Rid s Eqution, Itetive Resolution. Impotne of fist time step in te onvegene poess. If ell seleted, onvegene is ieved in 5-0 itetions. 2. Pmete omogeneiztion : Depending on te sitution : Aitmeti men o enteed vlue. Aitmeti seems to led to bette peision?!. 3. No poblem it te onvetive (gvity tem. 4. Heteogeneous medi my led to toug situtions : flo fom n impevious medium (lo yduli diffusivity, D=K/C to pevious medium (ig D 5. No poblem it te unstuted - stuted tnsition (see 4.
Ai Wte Migtion ( ( ( + θ ω = θ = θ α = ρ µ = = + ρ + θ + θ ρ = + θ α α α α α P P ( P, ( g P. k U 0 U H U. t H ( ( 0.U t : te, : i, H : Heny s onstnt, _ : density, _ : poosity
Ai Wte Migtion ( ( ( ( ( ( + θ ρ + θ + ρ θ ρ = ρ θ + ω θ + + θ = z ( K H. ( K H ( K. t C( H ( t H ( z ( K. t C( : te, : i, H : Heny s onstnt, _ : density, _ : poosity Reitten in tems of pilly ed & i ed pessues
Ai Wte Migtion It looks like te mtix system A A, 2, =.(U + A 2,2, + U,2 =.(U 2, + U 2,2 And solved te sme y s fo Rid s eqution A A,i,,i 2,2,i+,i Δt + Δt n n +.U +.U,i+,,i+ 2,2 =.U =.U,i,2,i+ 2, A,i 2,,i+ Δt n
Toum nd Vulin expeiment Exp. nd num. Anlysis of to-pse infilttion in ptilly stuted soil, TIPM ( 986 «Vetil infilttion in sndy olumn it no ltel i flo nd it i flo» Ai flo No Ai flo. Te te infilttion is dstilly sloed don by i
Ai mss onsevtion poblem THM nd Geo. Beviou of ly bie in diotive ste epositoies, Volket et l, CCE Repot EUR 6744 en 996 «Vetil infilttion in sndy olumn it no i flo nd soluble i» Ai mss onsevtion Ai Pessue
EVEGAS Euopen Pojet Cnnot find te efeene! «Podution of Hydogen Bubble in stuted poous medi» H2 H2 Does you ode see te Ct s es?
+ Tempetue : Temo Hydulis Pysis Numeis A A A A A A, 2, 3, C C C I+ N+, I C, t I+ N+, I A 2, 2 N+, I 3, 3 t T t I+ + A + A + A,2 2,2 3,2 T +.(U A A + A + A 2,3 3,3, + U T T,2 + U +.(U +.(U,3 2, 3, + U + U + U, +. U =.( U + U + U A,4 2,2 3,2 N+ I+ N N N N+, I,, 2, 3, 4, 2 = 0 + U + U T t 2,3 3,3 + U + U I N, I N N N N N, I C N +. U2, 2 =.( U2, + U2, + U2, 3 + U2, 4 A 2, A, t I I N+, I+ N N N N+, I C +, +. U3, 3 =.( U3, + U3, 2 + U3, 4 A 3, A 3, 2 t t I 2,4 3,4 = 0 = 0 + + + +, I 2 3 N I A T I t
Pollok D.W., WRR (22 986 «Simultion of Fluid Flo & Enegy Tnspot Poesses ssoited it HLW Disposl in Unstuted Alluvium» Pollok s esults Stution 00-000 y. Stution 0-00 y.
Conlusion Itetive metods ok quite ell. Not time onsuming, s omped to «globl» metods (benmking it THM odes. Sometimes gndm s tiks (oie of good vibles must be intodued My beome vey toug it some medi su s unstuted flo in ftued medi, geotemy, If you do it, ve fun! A ppe is in peption!