Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des...
Transcript of Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des...
![Page 1: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/1.jpg)
Habilitation a diriger des recherches
Modelisation et simulation numeriquede quelques problemes de transport ;multiresolution, problemes inverses
Marie Postel
Laboratoire Jacques-Louis LionsUniversite Pierre et Marie Curie, Paris 6
1er octobre 2008
![Page 2: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/2.jpg)
Short CV1985 Docteur Ingenieur Ecole Centrale des Arts et Manufactures (D.Aubry)
1985-1987 Geophysicist, E. P. Schlumberger
1987-1990 Post-doc Courant Institute
1990-2008 Associate professor LJLL� 1994-1996 visiting scientist IWR Heidelberg
Themes of research1987-1994 Waves in [1,2,3,4,5,6,7]
random media [8,9,10]
1994-2008 Chromatography [31,32,33][34,35,36,37]
1997-2004 Multiresolution [11,12,13][14,15,16,17]
2003-2008 ERT Transportation [20,23],[21,22]of hydrocarbons [24,25,26,27,28,29,30]
![Page 3: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/3.jpg)
Outline
1 Waves in random medium
2 Chromatography
3 MultiresolutionMultiscale analysisAdaptive scheme for hyperbolic PDEs
4 Transportation of hydrocarbonsMathematical model and numerical schemeAdaptive grid and time step refinement
![Page 4: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/4.jpg)
Waves in random medium
Work in collaboration with G. Papanicolaou (Stanford U.), B.White (Exxon), P. Sheng (Hong-Kong), W. Kohler (Virginia Tech),Marc Asch (UPJV), J.P. Fouque (UCSB), J.F. Clouet (CEA)
Propagation of acousticwaves in a layered mediumSeparation of scales :
slow smooth variationsO(1)
pulse width O(ε)
fast random variationsO(ε2)
![Page 5: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/5.jpg)
Numerical simulation of waves in random mediumStatistics of the reflected signal explained by a limittheorem when ε→ 0 (Burridge, Papanicolaou, Sheng, White, 1989),Simulation of the direct problemSpectral analysis techniqueInverse problem: O(1) fluctuations of the mean properties
Plane wavesSpherical waves
(Asch, Kohler, Papanicolaou, P., White. SIAM Rev., 1991)
![Page 6: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/6.jpg)
Outline
1 Waves in random medium
2 Chromatography
3 MultiresolutionMultiscale analysisAdaptive scheme for hyperbolic PDEs
4 Transportation of hydrocarbonsMathematical model and numerical schemeAdaptive grid and time step refinement
![Page 7: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/7.jpg)
Chromatography
Work in collaboration with M. Sepulveda (Concepcion, Chile), W.Jager (Heidelberg, Allemagne), F. James (MAPMO, Orleans), M.Schoenauer (INRIA Futur), A. Auger (INRIA Futur), M. Jebalia(LJLL)
��������������������
��������������������
����������������
����������������
������������
������������
����������������
������������
������������
������������
������������
������������
������������
������������
������������
������������������������������������������������������������������������������������������
������������������������������������������������������������������������������������������
��������������������������������������������������������������������
��������������������������������������������������������������������
������������
������������
������������������������������������������������
������������������������������������������������
����������������
����������������
��������������������������������������������
��������������������������������������������
��������������������������������
��������������������������������
����������������������������������������������������������������������������
����������������������������������������������������������������������������
solvantinjection
A+B BA B
AB
A B
t
Industrial process(pharmacy, oilindustry, ...)
Used to separate oranalysemulticomponentchemical mixtures
![Page 8: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/8.jpg)
Modeling of the Chromatography processconcentrations c ∈ R
p, isotherm h : Rp → R
p
∂t(c + h(c)) + ∂x (uc) = 0, t ∈ [0, T ], x ∈ [0, L],c(0, t) = c inj(t),c(x , 0) = 0,
0
10
20
30
2 4 6
c(g/l
)
time(min)
BA:PE:MBA=1:3:1 vi=1ml
BAPE
MBA
0
10
20
30
2 4 6
c(g/l
)
time(min)
BA:PE=3:1 vi=0.5ml
BAPE
MBA
(James,Sepulveda,1991)
Other models for finite speed adsorption (P., Sep´ulveda, 1997, James,
P., Sepulveda, 2000)
![Page 9: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/9.jpg)
Chromatography parameters identification
Objective function : J(h) =12
∫ T
0||ch(L, t)− cobs(t)||
2 dt
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1040
60
80
100
120
140
160
180
K
N*
� Identificationof h characteristicof the mixture� Ill posed
identification problem
� Gradientestimationfor identificationof h by deterministicmethod, (James, P.,
2007)� Identification of h by stochastic optimisation methods (Jebalia,
Auger, Schoenauer, James, P., 2007)
![Page 10: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/10.jpg)
Outline
1 Waves in random medium
2 Chromatography
3 MultiresolutionMultiscale analysisAdaptive scheme for hyperbolic PDEs
4 Transportation of hydrocarbonsMathematical model and numerical schemeAdaptive grid and time step refinement
![Page 11: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/11.jpg)
Multiresolution method
Work in collaboration with A. Cohen (LJLL), S. M. Kaber (LJLL),M. Campos-Pinto (IRMA), S. Muller (Aachen), N. Dyn (Tel-Aviv)
Context: systems of hyperbolic PDEs
∂tu + DivxF (u) = 0
Difficulties:
singularities
theory : entropy weak solutions
numerical analysis : costly schemes and limited order ofconvergence
![Page 12: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/12.jpg)
Multiscale analysis
Answer: adaptive mesh refinementDifficulties:
implementation : displacement of the singularities → meshconvergence analysis
Existing approaches:AMR (Adaptative Mesh Refinement) (Berger, Oliger, ...)Adaptive multiresolution flux evaluation (Harten, Abgrall,Chiavassa, Donat, ...)Galerkin methods in wavelet spaces (Dahmen, Cohen,Masson, Maday, Perrier, Bertoluzza, ...)
Fully adaptive multiresolution schemeHarten’s discrete multiresolution framework and link withwavelet theoryFinite volume scheme on a time adaptive grid
![Page 13: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/13.jpg)
Discrete multiresolution frameworkIn 1D Discretization of a function u(x) by its� point values uj ,k = u(k2−j)
� mean values uj ,k = 2j∫ (k+1)2−j
k2−ju(x)dx
x=0 x=L∆
∆
x
x/2
Dyadic hierarchy of gridsSj , j = 0, . . . , J
In 2D
![Page 14: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/14.jpg)
Encoding / Decoding (mean value discretization)Projection operator P j
j−1 of Sj on Sj−1: Uj−1 = P jj−1Uj .
j
j−1k
2k2k+1
uj−1,k = 12(uj ,2k + uj ,2k+1)
Prediction operator P j−1j from Sj−1 to Sj : Uj = P j−1
j Uj−1
j−1 j
kk−1 k+1 2k 2k+1
is an approximation of Uj
� local� consistentP j
j−1P j−1j = I
� exact for polynomialsof degree r = 2n
uj ,2k = uj−1,k +
n∑
l=1
γl(uj−1,k+l − uj−1,k−l
)
![Page 15: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/15.jpg)
Multiscale decomposition
Prediction error:uj ,2k − uj ,2k
(consistence ⇒ )= uj ,2k+1 − uj ,2k+1
Details :dj ,k = uj ,2k − uj ,2k
errors on all but onesubdivisions
Details Dj
+Mean values Uj−1
m
Mean values Uj
UJ ⇔ (UJ−1, DJ)
⇔ (UJ−2, DJ−1, DJ)
...
⇔ (U0, D1, . . . , DJ)
= MUJ = (dλ)λ∈∇J
UJ MUJ
Encoding
Decoding
complexity
in O(Card(SJ))
![Page 16: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/16.jpg)
Biorthogonal wavelet framework
u =∑
λ dλΨλ,
dλ =< u, Ψλ >,Ψλ = 2j/2Ψ(2j .− k),λ = (j , k).
Ψ∼
U2
U1
U0
D2
D1
∣∣dj ,k∣∣ ≤ C2−js
Measure of the local smoothness if u ∈ Cs(Supp Ψj ,k ),Polynomial exactness of degree r with s ≤ r .
Thresholding u → TΛu :=∑
λ∈Λ dλΨλ
Set of significant details Λ = Λε := {λ, |dλ| ≥ ε}.
![Page 17: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/17.jpg)
0 0.5 1
−0.8
−0.4
0
0.4
x
U
Function U on finest grid
U
0 0.5 10
2
4
6
x
j
significant details
0 0.5 1x
Adaptive grid
0 0.5 1
−0.8
−0.4
0
0.4
x
U
Uε reconstructed on finest grid
UUε
0 0.5 1
−0.8
−0.4
0
0.4
x
U
Uε on adaptive grid
![Page 18: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/18.jpg)
� Thresholding
Γε = {λ ∈ ∇J , |dλ| ≥ ε|λ|}, εj = 2d(j−J)ε.
Tε(dλ) =
{dλ if λ ∈ Γε
0 otherwise
AεUJ =M−1TεMUJ .
� Adaptive grid Sε = S(Γε)
(dλ)λ∈Γε←→ (Uλ)λ∈Sε
.
Tree structure for the adaptive tree Γε ⇒ Complexity of thecoding / decoding algorithm in O(Card(Γε)).
![Page 19: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/19.jpg)
Control of error
||AεUJ − UJ ||L1 ≤ Cε
ε−17
10−1
10
−1610
−1210
−810
−410
010
f1f2f3f4
0 100 200 300 400 5000
2
4
6
8x 10−3
smooth function
10 levels
←r=0
←r=2←r=4 ←r=6
Error versus compression
0 50 100 1500
0.002
0.004
0.006
0.008
0.01discontinuous function
10 levels←r=0
←r=2←r=4
←r=6
Error versus compression
![Page 20: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/20.jpg)
Extension to multi dimension
� Cartesian grids: tensorization of 1D multiresolution
level j−1
level j
0
1 2
3
0
3
2 1
� Prediction for 2Dtriangular discretizations
� Polynomial reconstructionof order 1 diverges
� Stable reconstruction bylinear combination
8>><>>:
bui,0 = ui−1,0,bui,1 = ui−1,0 + (ui−1,2 + ui−1,3 − 2ui−1,1)/6,bui,2 = ui−1,0 + (ui−1,1 + ui−1,3 − 2ui−1,2)/6,bui,3 = ui−1,0 + (ui−1,1 + ui−1,2 − 2ui−1,3)/6.
Cohen, Dyn, Kaber, P., JCP, 2000
![Page 21: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/21.jpg)
Adaptive scheme for hyperbolic PDEs
∂tu + DivxF (u) = 0.
Reference scheme on finest grid SJ :
Un+1J = EJUn
J
Un+1λ = Un
λ + B(Unν ; ν ∈ Vλ ⊂ SJ), λ ∈ SJ
Harten’s algorithm still on finest grid SJ
Un+1ε,λ = Un
ε,λ + Bε(Unε,ν; ν ∈ Vλ ⊂ SJ), λ ∈ SJ
Fully adaptive algorithm on S(Γnε)
Un+1ε = EεUn
ε .
![Page 22: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/22.jpg)
Evolution of (Unε , Γn
ε) into (Un+1ε , Γn+1
ε )
� Encoding of thesolution at time tn
� Thresholding to obtainthe tree Γn
ε
� Prediction
� Decoding� Evolution
Solution at time tn
0
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 1
Ugrid
Details at time tn
-0.4
-0.2
0
0.2
0.4
0 0.2 0.4 0.6 0.8 1
D3D2D1
![Page 23: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/23.jpg)
Evolution of (Unε , Γn
ε) into (Un+1ε , Γn+1
ε )
� Encoding
� Thresholding toobtain the tree Γn
ε
� Prediction of the treeΓn+1
ε and grid S(Γn+1ε )
(Harten heuristics ormore refined strategy)
� Decoding� Evolution
Details at time tn
-0.4
-0.2
0
0.2
0.4
0 0.2 0.4 0.6 0.8 1
D3D2D1
Tree Γn+1ε and grid S(Γn+1
ε )
0
0.5
1
1.5
2
2.5
3
3.5
4
0 0.2 0.4 0.6 0.8 1
![Page 24: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/24.jpg)
Evolution of (Unε , Γn
ε) into (Un+1ε , Γn+1
ε )
� Encoding
� Thresholding� Prediction
� Decoding on the gridS(Γn+1
ε )� Evolution with fluxescomputed usingreconstructed solution
Tree Γn+1ε and grid S(Γn+1
ε )S
0
0.5
1
1.5
2
2.5
3
3.5
4
0 0.2 0.4 0.6 0.8 1
Solution at time tn+1
0
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 1
Ugrid
![Page 25: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/25.jpg)
Error analysis
TheoremCohen, Kaber, Muller, P., Math. Comp., 2003
||UnJ − Un
ε ||L1 ≤ Cnε
Hypotheses
More demanding rules for the refinement of the tree
Smoothness of the underlying wavelet
Scalar equation, on 1D or cartesian multi-D grid
![Page 26: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/26.jpg)
Space shuttle reentry problem
![Page 27: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/27.jpg)
Outline
1 Waves in random medium
2 Chromatography
3 MultiresolutionMultiscale analysisAdaptive scheme for hyperbolic PDEs
4 Transportation of hydrocarbonsMathematical model and numerical schemeAdaptive grid and time step refinement
![Page 28: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/28.jpg)
Transportation of hydrocarbons
Work in collaboration with F. Coquel (LJLL), Q.-H. Tran (IFP), N.Andrianov (Schlumberger), N. Poussineau, L. Linise, Q.-L.Nguyen. (Equipe de Recherche Technologique)
������������������
������������������
���������������������������������������
���������������������������������������
��������������������������������
������������������
������������������
100 to 300 km
Transient Simulation
������������������������
1000 to 5000 m
Gas
20 to 1000 m
WaterLiquid
![Page 29: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/29.jpg)
Technological bottleneck: slugging.
![Page 30: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/30.jpg)
Drift Flux Model (DFM)Conservation laws
∂t(ρ) + ∂x (ρv) = 0 (1a)
∂t(ρv) + ∂x (ρv2 + P) = ρS (1b)
∂t(ρY ) + ∂x (ρYv + σσσ) = 0 (1c)
thermodynamical and hydrodynamical closure lawsexpressing
P(u) = ρY (1− Y )φ2(u) + p(u) (2a)
σ(u) = ρY (1− Y )φ(u) (2b)
as a function of u = (ρ, ρv , ρY ), with φ(u) the drift (vg − vℓ )
hyperbolic system with eigenvalues :λ1(u)≪ λ2(u)≪ λ3(u)
Abstract form
∂tu + ∂x f(u) = ρs(u) (3)
![Page 31: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/31.jpg)
Modeling difficulties
Eigenvalues on different scalesQuantity of interest: gas mass fraction (slow wave CFL ≃ 1)Fast acoustic waves of less interest CFL ≃ 20
→ semi-implicit scheme with large time step (Faille, Heintze,’99,
Baudin, Coquel, Tran,’05)
Time varying coupling boundary conditions
Expensive hydro and thermo closure laws→ Multiresolution
![Page 32: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/32.jpg)
New answer: coupling multiresolution with largetime step methods
Treatment of thenonlinearities byrelaxation
Implicit/explicit eulerianscheme Andrianov, Coquel, P.,
Tran, 2006
Lagrange Eulerdecomposition forselective implicitation
Multiresolution Coquel, P.,
Poussineau, Tran, 2006
Local time stepping
![Page 33: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/33.jpg)
Taking care of the nonlinearities
Initial system
∂t(ρ) + ∂x(ρv) = 0 (4a)
∂t(ρv) + ∂x(ρv2 + P) = ρS (4b)
∂t(ρY )+ ∂x(ρYv + σσσ) = 0 (4c)
Relaxed system, for λ ≥ 0,
∂t(ρ)λ + ∂x(ρv)λ = 0 (5a)
∂t(ρv)λ + ∂x(ρv2 + Π )λ = ρλS(uλ) (5b)
∂t(ρΠ)λ + ∂x(ρΠv + a2v)λ = λρλ[P(uλ)− Πλ] (5c)
∂t(ρY )λ + ∂x(ρYv + Σ )λ = 0 (5d)
∂t(ρΣ)λ + ∂x(ρΣv + b2Y )λ = λρλ[σσσ(uλ) − Σλ] (5e)
Abstract form :∂tU
λ + ∂xF(Uλ) = ρS(Uλ) + λ[Euλ − Uλ]
![Page 34: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/34.jpg)
Properties of the relaxed system
Theorem ( Baudin-Berthon-Coquel-Masson-Tran)The relaxed system (5a)–(5e) is always hyperbolic. Itseigenvalue fields are all linearly degenerate
v − aτ, v − bτ, v , v + bτ, v + aτ, (6)
with τ = ρ−1.
TheoremThe relaxed system is a good approximation of the initialsystem in the neigborhood of λ→ +∞ under thesub-characteristic conditions
a >
√−Pτ (u) + P2
v (u) et b > |σY (u)| (7)
Linear asymptotic analysis by Chapman-Enskog
![Page 35: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/35.jpg)
ALE decomposition of the relaxed system
Moving frame with speed w .
∂tx|χ = v −w .
Dilatation rate J = ∂χx|t
∂t(J) + ∂χ(w) − ∂χ(v) = 0 (8a)
∂t(ρJ) + ∂χ(ρw) = 0 (8b)
∂t(ρvJ) + ∂χ(ρv w) + ∂χ(Π) = ρJS(u) (8c)
∂t(ρΠJ) + ∂χ(ρΠ w) + ∂χ(a2u) = λρJ[P(u)− Π] (8d)
∂t(ρYJ) + ∂χ(ρY w) + ∂χ(Σ) = 0 (8e)
∂t(ρΣJ) + ∂χ(ρΣ w)︸ ︷︷ ︸ + ∂χ(b2Y )︸ ︷︷ ︸ = λρJ[σ(u) − Σ] (8f)
Euler Lagrange
![Page 36: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/36.jpg)
Graphical interpretation of the ALE decomposition
x
t
v − aτ
v + aτ
v − bτ v + bτv
=
x
t
v − aτ
v + aτ
+
x
t
v − bτ v + bτ
Lagrange∆t∆x |v
n ± bτn| < 12
+
x
t
v
projection∆t∆x |v
n♯| < 12
![Page 37: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/37.jpg)
Analysis of the numerical scheme
Theorem ( Coquel, Nguyen, P., Tran, 2008)The numerical scheme derived from the Lagrange-projectiondecomposition is conservative. It is consistent and positive inthe following sense
ρn+1i > 0 and Y n+1
i ∈ [0, 1] (9)
if∆t∆x
<1
M (un, BCn, ∂tBCn)(10)
where the bound M (un, BCn, ∂t BCn) can be computedexplicitly at time tn and depends on the slow speed.
![Page 38: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/38.jpg)
Selective semi-implicit scheme for theLagrange-projection formulation
Uniform grid implementation
Numerical performances comparable to the pure Euleriansemi-implicit scheme
Improved robustness thanks to ∆t control
Adaptive grid implementation
Local grid refinement near the transport wavediscontinuities
Acoustic waves : not interesting, smoothed out by theimplicit 3-point scheme
Prediction of the grid from one time step to the next basedon the slow transport speed
![Page 39: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/39.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 40: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/40.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 41: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/41.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 42: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/42.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 43: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/43.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 44: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/44.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 45: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/45.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 46: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/46.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 47: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/47.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 48: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/48.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 49: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/49.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 50: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/50.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 51: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/51.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 52: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/52.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 53: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/53.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 54: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/54.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 55: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/55.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 56: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/56.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 57: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/57.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 58: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/58.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 59: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/59.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 60: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/60.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 61: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/61.jpg)
Example on v-shaped pipegrid level , uniform solution and adaptive solution
200
400
600
800
1000
0 2000 4000
x
Density
0
5e+06
1e+07
0 2000 4000
x
Pressure
0.003
0.005
0.007
0 2000 4000
x
Gas mass fraction
0.02
0.04
0 2000 4000
x
Debit gaz
Restart film Skip film
![Page 62: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/62.jpg)
From multiresolution to locally adaptive timestepping
Standard MR scheme
At a given time, same time step used for all cells
Time step ruled by CFL condition
∆tn <minj ∆xj
M (un, BCn, dtBCn).
Goal: use a time-step adapted to the local size of the cellconstant λ = ∆t/∆x(Berger, Colella ’89,Mueller, Stiriba ’06)
∆2∆
2∆∆t
t
x x
![Page 63: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/63.jpg)
Local time stepping implementation
Intermediate time stepping in fine/coarse transition zone :rely on flux conservation to synchronize the differentlevels solution.
Fnj,k−3/2
12 Fn
j,k−3/2
12 Fn
j,k−3/2
Fnj,k−1/2
Fnj+1,2k−1/2
Fn+1/2j+1,2k−1/2
x
t
Partial stencil of influence for the implicit Lagrangian step(ratio between the acoustic and kinetic speeds)
x
t
![Page 64: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/64.jpg)
Test case with smoothly varying initial condition
400
450
500
550
0 10000 20000 30000 0
2
4
6
ρ
x
UNIMRLTSt=0
LTS grid
![Page 65: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/65.jpg)
Error versus computing time gain - simplistic statelaw
![Page 66: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/66.jpg)
Error versus number of calls to state laws gain
1e-07
1e-06
1e-05
1e-04
0.001
0.01
0.1 1 10 100 1000
E(Y)
# State laws gain
MRLTS
![Page 67: Habilitation a diriger des recherches` - UPMC · 2017-05-09 · Habilitation a diriger des recherches` Modelisation et simulation num´ erique´ de quelques problemes de transport](https://reader030.fdocument.pub/reader030/viewer/2022040823/5e6d5136241fb16f2e44f9c8/html5/thumbnails/67.jpg)
Conclusion
Extensive study of the multiresolution schemeperformances : compression, precisionApplication to the industrial problem of hydrocarbontransportation :
Control of positivity for a large time step scheme withboundary conditionscostly flux functions make multiresolution really rewarding
PerspectivesStability condition for the full model with driftAnalysis of the LTS algorithm