v
ABSTRACT

Nowadays, numerical techniques become the effective tools to solve the
problems in science and engineering. Eventhough the impressive recent progresses
attained in computer technologies and computational simulation techniques,
numerous models intractable when the usual and well-experienced discretization
techniques are applied for their numerical simulation due to their high complexity
and requirements. One of the typical difficulties is highly multi-dimensional models
arising from quantum mechanics or kinetic theory descriptions of solids and

the number of degrees of freedom involved scales exponentially with the dimension
of the space concerned.
In order to overcome the drawbacks above, one lastest technique in recent
years proposed to support, activate in using the mesh-based discretization
techniques -FEM -is called Proper Generalized Decomposition (PGD). This is a
powerful model reduction technique by means of successive enrichment a separated
representation of the unknown field, so the computational complexity of the PGD
scales linearly with the dimension of the space.
And a coupling Proper Generalized Decomposition and Finite Element
Method  PGD-FEM briefly  will open a new approach in searching a powerful
kind of simulation techinique in both terms of computing time and accuracy.
ized Decomposition and Finite

Eventhough the topic just started to invest PGD-FEM for fluid problem in a
small term of the pressure Poisson equation from 2D unsteady imcompressive
Navier-Stokes flow, the comparative results speaked out the outstanding innovative
property of PGD-FEM in both computing time and accuracy from the traditional
vi
discretization technique (FEM). Moreover, in order to overcome its remaining
drawbacks and enlarge, develop further research trends, I also provided to solve

sut 2D cho bài toán Navier-Stokes ca dòng chy nht không nén ph thuc vào
thng hng nhu kin biên
viii
hn hp ( Dirchlet-ng kt qu c y s t
khi gii quyt b-FEM v mt th chính
xác so vi rc truyn thng thi vi mong mun to
mt s thun li trong vic hoàn thi rng, phát trin  tài
i, tác gi  c n vic gi-Stokes
cho dòng chy nht không nén ph thuc thi gian v u kin biên lid-driven
cavity ba trên k thut tham chiu Chorin-Temam.

ix
MC LC
TRANG
Trang ta
Quy tài
Lý lch cá nhân i
L iii
Cm t iv
Tóm tt v
Mc lc ix
Danh mc kí hiu-t vit tt xi
Danh mc hình v xii
NG QUAN 1
1.1 Tng quan v ng nghiên cu 1
1.2 Mu, khách th ng nghiên cu 2
1.3 nh nhim v và phm vi nghiên cu c tài 3
1.4 u 3
 LÝ THUYT 4
2.1 Generalized Decomposition (PGD) 4
|| . ||
2

Res
n


x


y




D


N

H
1n


Min kho sát ca bài toán

Min kh

Min kh

Biên ca mi

Biên Dirichlet ca 

Biên Neumann ca 

Không gian hàm Sobolev mà có giá tr trit tiêu trên 
DVec-n ng ra ngoài) 
N

Vec-ng ti các nút trên mi

Vec-ng ti các nút trên mi
x

Vec-ng ti các nút trên mi
y

Tng s nút trên min 


FEM

PGD

PGD-FEM MBS

ROM

LATIN

POD

SVD

PDE

Vec- hàm Y,S,G ti các nút trên mi
y

Toán t gradient
,
pp
xy



coupling Proper Generalized Decomposition and Finite
Element Method

Multi-Bead Spring

Reduced-Order Model

LArger Time INcremential

Proper Orthogonal Decomposition

Singular Value Decomposition

Partial Differential Equations xiii
DANH MC HÌNH V

Hình 2.1 Min khu kin biên ca 
Hình 3.1 Min khng nht c
Poisson
Hình 3.2  th c
gii tích, FEM, PGD-ng kiu i 30 x 30 và 80 x 80.
Hình 3.3  th cson vi f(x,y)=x
2
-y
2


c lc cho các bài toán trong khoa hc và k thut. Mng tin b, phát trin
t bc trong công ngh máy tính và nhng k thut tính toán s u bài toán vn
còn b hn ch c bit v mt thi gian tính toán khi  ri rc khó có
th gii quyt bi tính phc tp và m yêu cu ngày càng cao ca bài toán.
Có th nêu ra mt s v p phi là:
(i)bài toán có s chiu không gian kho sát lng gp  ng t,
thuyng hc ct phc tp [8], hay sinh hc, hóa hc [16]. V
khi áp di r phc tp ca bài toán  t l i s
chiu không gian cn kho sát min thi gian th
kho sát giàn khoan ph thuc vào thi gian thc [13].(iii) bài toán có min kho sát suy
bin xut hin trong thanh, tm, vn nhng thông s, tham s khác
(ngoài yu t không gian-thi gian vt lý ng) ví d kho sát h s truyn nhit
ca vt liu trong bài toán truyn nhit.
 khc phc nhng v trên, mi i n
n b try trong vic phi hp vi rc truyn th
pháp Proper Generalized Decomposition (PGD), mt mô hình gim bc bài toán
da trên c s tách bin giúp gim mt cách hiu qu  phc tp ca bài toán vi t l
tuyn tính theo s chiu không gian so vi t l i rc truyn
thng nên nó mang lt tri trong thi gian x m b chính
xác so vi rc truyn thng.
Mt s công trình nghiên cu quc t ni bt quan trng liên quan trc tip n quá
trình nghiên c tài là:

2
[4]: s dng  gii quyt v nhiu chiu cho bài toán Poisson
vu king nhng gp trong thuyng hc ct phc t
m rng cho bài toán MBS (multi-bead-spring) ca không gian hai, ba chiu.
[13]: ng d   gii quyt bài toán Navier-ng
hp lid-driven cavity vi các h s Reynolds khác nhau và so sánh kt qu v
gii tích da trên hai tiêu chí: thi gian tính toán  chính xác.

t.
 Lp trình tính toán và mô phng kt qu bng ngôn ng lâp trình k thut Matlab
trên máy tính.
 So sánh kt qu gi-FEM vi u khác.
Phm vi nghiên cu:
Do còn nhng hn ch nhnh trong qúa trình nghiên cu t phía tác gi
s xut hin rt mi c l s gii hn
nghiên cu mt phn nh c   -Stokes áp dng cho dòng chy nht
không nén ph thuc thi gian trong không gian 2D bng vic gii quy  
Poisson áp sut vng hp khác nhau cu kin biên b-
ng thi ver-Stokes cho dòng chy nht không nén ph thuc
thi gian trong không gian 2D, tác gi  tip cn mô hình gii b
FEM du Chorin-Temam, tu ki hoàn thin, phát
tri u 
u.
S dng phn mm Matlab h tr vic lp trình tính toán và mô phng trên máy tính
laptop có cu hình trung bình.
Thc hin phép so sánh kt qu gi-FEM v
chiu khác theo tiêu chí th chính xác. 
4
2
T

Decomposition (PGD)[5,6,7,9]

Ngày nay vic mô phng s hóa cho nhng h tht phc t
i ngày càng cao trong khi khó có th gii quyt mt cách d dàng bi các


Cùng vi s phát trin ca trên ROM, ta có th nhn
mt trong nh bin nht da trên mô hình gim bc này c
 (Proper Orthogonal Decomposition)- mô hình tách bin d
trc giao- c hình thành da trên vic thit lp mt ma trng vi mi
m thi gian ri tin hành tìm tr riêng và các vec-
i

n
cho
bài toán (Q Q
T
 

5
Gi s u(x,t) là hàm tha mãn yêu cu bài toán vi x

IR
D
(D=1,2,3) và t


IR
+
ti mi x
i
, ta có bin thng t
p
i i=1: N
n

Q
i
i
i i N
x ,x , ,x x xx



(2.2)

i
là bin không gian, thi gian hoc tham s bt k ca bài toán
thuc min kh

IR
d
ng v
Công thc (2.2t tri trong thi gian x lý,
c th nu mi bin x
i
ri rc thành M bc t do thì tng s bin cn tìm là Q x N x
1
12
1 1 1
2
22
2
12
x
x 
y
=IR
2
vu ki
ca bài toán. (2.3)
Tìm U(x,y)
 vi phân, g là thành phn th hai ca bài toán.
t, PGD là mgii lp m c nh (fixed interation
method)dùng  tìm nghim xp x U(x,y) vi:
2
U( ) IR , IR , IRx,y X Y x X y Y       

Gi s, tc lp th n, hàm F
i
và G
i
t.
Bây gi ta mun tìm hàm F
n
, G
n
.
Gc biu din tc lp th n 
1
1
U ( , ) F ( )G ( ) F ( )G ( )
n
n i i n n

hàm F
n
, G
n
, ta thc hin phép tham chiu cho tng bin F
n
,G
n
vào (2.5), ta có:

7
2
2
2
1
L (X)
L (X)
1
L (X)
L F ( )G ( ) F ( )G ( ) ,F g,F Res ,F
n
i i n n n n n n
i
x y x y



  



2
(Y)
là tích trong trên L
2
ng theo mi 
ng thi, t phn sai s Res
n
phi vuông góc vi hàm F
n
, G
n
nên:
2
2
1
L (X)
1
L (X)
L F ( )G ( ) F ( )G ( ) ,F g,F
n
i i n n n n
i
x y x y








n
,G
n
có th
c s dng bi hoc      c phn t hu hn
(FEM).Và  tài nghiên cu này, tác gi s s dng cho
vic tính toán F
n
,G
n
.   m mu ch th hin tính kt hp gia
n t hu hn (
tm gi tt PGD-FEM)  u.
 gii thut tng quát c


8



9
2.1.4 Tiêu chu
n cht c-FEM là mt quá trình
gii lp vm c nh (Fixed Point Iteration) nên nhng công thc tiêu chut
ra trong vi sai s t vai trò quan trng và quy chính
xác c-FEM so vu khác.
Có th phân loi mt s nhóm công thc da trên phn gii thut c
pháp PGD 
Dùng để kiểm tra điều kiện hội tụ của F
(k)
và G
(k)
.
( ) ( 1) ( ) ( 1)
( ) ( ) ( 1) ( 1)
F F G G
F G F G
k k k k
k k k k




   


(2.8)
Dùng để kiểm tra điều kiện hội tụ của U
n

xy
ror x y x y dxdy




(2.9)
Trong ph tai nghiên cu, u
ref
có th là nghim ci tích
ho

10
Dùng để kiểm tra độ chính xác của phương pháp PGD-FEM mà không có phương
pháp tham chiếu.
( ) ( )
1
L F G g
n
ii
i










  


11
[1, 12,15]

i quá mi m, xa l i vi
bt kì mt bài toán k thut nào  trong nhic vt
rn, t, truyn nhin t hng sinh hóa hi
i bt i vi kh  x lý
dng hình hc phc tu kin biên phc hp; phân tích các bài toán tuyn tính,
phi tuyng lc h
 tìm thông tin v 
n hin nay.
 ni dung lý thuyn
t hu hc tích hp khi tính toán gii F
k
,G
k
  gii thut ca
 Vì th n lý thuyt c tp
c gi làm rõ vic tính toán F
k
, G
k
.
 n, ta s  Poisson ca hàm u trong không
gian hai chiu:

22

12
2.2.1. Thit lp dng yu
Nhân hàm trng s u
*
(x,y)

H
1
2.10) và thc hin tích phân trên toàn
mi
22
**
22
uu
u u f( , )d x y d
xy



    





S dng tích phân tng phn, ta có:
**
**
u u u u
u u u f( , )d d d x y d

xy

 
n

Và u
*

D
do u
*
(x,y)

H
1

**
**
**
**
u u u u
u g( , )u u f( , )
u u u u
u f( , ) u g( , )u
N
N
d d x y d x y d
x x y y
d d x y d x y d
x x y y

nod
N
ii
i
N
T
i
i
x y x y
x y x y






Vu
V


V
Nnod
] là vec-a hàm dng ti các nút trên toàn mi.
u
Nnod
]
T
là vec-a giá tr hàm u(x,y) ti các nút trên toàn mi.
Xét phn t 
i

Thay (2.13) vào (2.12):
**
**
u u u u
u f ( , ) u g( , )
f ( , ) g( , )
e e N
e e N
e e e e
ee
ee
TT
TT
ee
d d x y d x y d
x x y y
d d x y d x y d
x x y y
   
   
   
      
   

   
       


   


 
 
 
-1
U = K F
(2.16)
15
Sơ đồ giải thuật tổng quát của phương pháp FEM

Khai báo d liu
bài toán
Tính toán ma tr cng, vec-ti
k,f  (2.14)
Lp ghép vào ma trn toàn cc
K, F  (2.15)
u kin biên
Nghim ca bài toán 
(2.16)


16

NG D và FEM CHO BÀI
T



u
: s i ca vn tc theo thi gian.
 
.uu
: s c hình thành bng vec-vn tc.