pgd in computational mechanics f. chinesta, a. leygue eads chair, ecole centrale of nantes reduced...
TRANSCRIPT
![Page 1: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/1.jpg)
PGD in computational mechanics
F. Chinesta, A. Leygue
EADS Chair, Ecole Centrale of Nantes
Reduced Bases Workshop
![Page 2: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/2.jpg)
The past
The dream
Computational mechanics
![Page 3: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/3.jpg)
1 – Models defined in high dimensional spaces;
2 – Separating the physical space;
3 – Parametric models: optimization, inverse methods and simulation in real time;
4 – Miscellaneous.
PGD & …
![Page 4: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/4.jpg)
1 - Multidimensional models
Ψ x,t,q1,L ,qN( )
1
11
, , ,..,, , ,..,
( )( )
N
jj
NN
j
tt
t
d
d
Ψ Ψ x q q
x q qq
q
(3 1 3 )N D
![Page 5: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/5.jpg)
Ψ q1,L ,qN( ) ≈ Qi
1 q1( )×L ×i=1
i=M
∑ QiN qN( )
JNNFM, 2006
![Page 6: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/6.jpg)
1
* * *
( , ) ( ) ( ) ( ) ( )
( , ) ( ) ( ) ( ) ( )
n
j jj
u x y X x Y y x yR
u x y R x S y R x y
S
S
*( ) ( ) 0x y
u u u dx dy
M A M A
Successive enrichment:
R S
R⋅S* M (u(R,S,X j ,Yj )) −A( )x×y
dx dy=0•
x
⏐ →⏐ ⏐ S
S R
R* ⋅S M (u(R,S,X j ,Yj )) −A( )x×y
dx dy=0•
y
⏐ →⏐ ⏐ R
![Page 7: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/7.jpg)
Variants:
- Newton linearization;
- Residual minimization in the case of non-symmetric operators;
- Other more efficient and optimal constructors are in progress.
![Page 8: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/8.jpg)
Transient: INNFM 2007
Ψ t,q1,L ,qN( ) ≈ Ti t( )×Qi
1 q1( )×L ×i=1
i=M
∑ QiN qN( )
Non linear models: AR 2007
M Ψn−1( )Ψn =A
Extension of the radial approximation, Ladeveze 1985
Non linear models: CMS 2010
Newton and advanced fixed point M Ψn(s)( )Ψn(s+1) =A
![Page 9: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/9.jpg)
BCF descriptions: MSMSE 2007
Quantum chemistry & Master Equation: IJMCE 2008, EJCM 2010
Complex flows:
Ψ t,x,q1,L ,qN( ) ≈ Ti t( )×Xi x( )×Qi
1 q1( )×L ×i=1
i=M
∑ QiN qN( )
ACME 2009, MCS 2010
![Page 10: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/10.jpg)
2 - Separating the space …
Ψ x, y,z( ) ≈ Xi x( )×Yi y( )×Zi z( )
i=1
i=M
∑
CMAME 2008
Boundary conditions and separating non hexahedral domains
IJNME 2010
![Page 11: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/11.jpg)
u(x, y, z,t)≈ F j ,1(x, y) ⋅F j ,2(z)
j=1
n
∑ ⋅F j ,3(t)
3D modeling in plate domains
CMAME, In press
3D solution with 2D complexity
![Page 12: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/12.jpg)
1
( , , ) ( , ) ( )n
j jj i i
i
u x y z X x y Z z
![Page 13: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/13.jpg)
Shells
xx
yy
zz
FEM solution:
PGD solution
610 dof
500 times steps, implicit
14s using matlab
![Page 14: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/14.jpg)
∂T∂t
−k∂2T∂x2
=S
T (t, x,k)
2
2
TS
tk
x
T
,( )T t x
3 - Parametric modeling
3D
2D
![Page 15: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/15.jpg)
allows the solution of a single problem instead the many solutions required by
usual strategies in:
Optimization & Inverse identification
( , , )T t x k
( , , )dT
t x kdk
is explicitly available
![Page 16: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/16.jpg)
accounting for variability
( , , )g gT t x kk
![Page 17: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/17.jpg)
Parametric solution computed off-line with all the sources of variability considered as extra-coordinates
On-line processing of such solution in real time and using light computing platforms: e.g. smartphones, …
![Page 18: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/18.jpg)
Parametric modeling: MCS 2010, ACME 2010, JNNFM 2010
Process optimization: Composites Part A, 2011
Inverse identification: JNNFM 2011
Control & DDDAS: CMAME, MCS - Submitted
Accounting for variability: CMAME In press
Using smarthphones: CMAME In press, CMAME Submitted
u(x,t, p
1,L , p
Q) ≈ Xi (x) ⋅T(t)
i=1
n
∑ ⋅Pi1 p1( )L Pi
Q pQ( )
Shape optimization: ACME 2010
![Page 19: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/19.jpg)
4 – Miscellaneous Error estimator: CMAME 2010
Coupling PGD & FEM: IJMCE 2011
Coupled models & multi-scale: IJNME 2010, IJMF 2010
![Page 20: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/20.jpg)
Time multi-scale & parallel time integration: IJNME 2011, In press
u =u(τ )
u =u(τ ,t)
![Page 21: PGD in computational mechanics F. Chinesta, A. Leygue EADS Chair, Ecole Centrale of Nantes Reduced Bases Workshop](https://reader035.vdocument.in/reader035/viewer/2022062422/56649f385503460f94c55306/html5/thumbnails/21.jpg)
PGD based transient boundary element discretizations: EABE 2011
Convective stabilization:- Separating and then stabilizing- Stabilizing and then separating
Computational Mechanics, Submitted