The influence of quadrature formula in 2D and 3D mortar methods

Yvon Maday

Description

Niveau: Secondaire, Lycée, Terminale
The in uen e of quadrature formulas in 2D and 3D mortar element methods Yvon Maday 1 , Fran es a Rapetti 1 , and Barbara I. Wohlmuth 2 1 Laboratoire d'Analyse Numerique, Paris 6 University, Bo^te ourrier 187, 75252 Paris edex 05, Fran e 2 Mathematis hes Institut A, Universitat Stuttgart, Pfaenwaldring 57, 70 569 Stuttgart, Germany Abstra t. The paper is on erned with the mortar nite element dis retization of s alar ellipti equations in three dimensions. The attention is fo used on the in uen e of quadrature formulas on the dis retization error. We show numeri ally that the optimality of the method is preserved if suitable quadrature formulas are used. 1 Introdu tion The mortar element method, rstly proposed in [3?, is a non- onforming non- overlapping domain de omposition method. The oupling of dierent phys- i al models, dis retization s hemes or non-mat hing triangulations at the interfa es an be eÆ iently realized in terms of mortar element methods. To preserve the global optimality of the lo ally adapted dis retizations, the in- terfa es between the dierent regions have to be handled appropriately. Due to its high exibility, this approa h has been analysed and implemented in many situations.

• all ontinuous pie


• ewise linear


• mortar side


• fun tion


• ontinuity ondition


• pie ewise


• linear nite


• ondition prevents

Informations

• Publié par: profil-zyan-2012
• Nombre de lectures: 15
• Langue: English
• Poids de l'ouvrage: 1 Mo

Extrait

The

of
quadrature
the
error
m
form
standard
ulas
An
in

2D
b
and
(see
3D
uit
mortar
eak
elemen
olv
t
sions,
metho
In
ds
tegrals
Yv
based
on
error
Mada
ds
y
the
1
implemen
,
b
F
traces.

non-matc
Rap
the
etti
orts
1
that
,
supp
and
quadrature
Barbara
men
I.
in
W
master
ohlm
w
uth
the
2
mortar
1
the
Lab
the
oratoire
a
d'Analyse
asp
Num
realization

in
erique,
terms
P
the
aris
in
6
dened
Univ
as
ersit
the
y
t
,
t
Bo
easily
^
due
te
one

in
187,
trian-
75252
task.
P
ev
aris


ev
05,
and
F
a
rance
sla
2
not

umerical
hes
b
Institut
form
A,
side
Univ
master
ersit
t
at
to
Stuttgart,

Pfaen

w
eleton
aldring
osition
57,
eak
70
ortan
569
for
Stuttgart,
is
German
the
y
across

It
The
written
pap
a
er
tiplier
is
of

asso
with
in
the
discrete
mortar
dieren
nite
grids
elemen

t
es
discretization
tersection
of
orts.
scalar
o
elliptic
nite
equations
the
in
b
three
tersected
dimensions.
in
The
the
atten
latter
tion

is
dimensions,
fo
of

(union
on
or
the
a

use
of
ulas
quadrature
these
form

ulas
the
on
Ho
the
using
discretization
appro
error.
the
W
b
e
form
sho
on
w
e
n
do
umerically
optimal
that
for
the
in
optimalit
dimensions).
y
appro
of
a
the
based
metho
v
d
the
is
on
preserv
As
ed
elemen
if
metho
suitable
is
quadrature
replace
form
exact
ulas
tin
are
y
used.
at
1
sk
In
of
tro


with
The
w
mortar
one.
elemen
imp
t
t
metho
ect
d,
the
rstly
tation
prop
the
osed
of
in
w
℄

is
the
a
terfaces.


non-
e
o
in
v
of
erlapping
Lagrange
domain
ul-

and
osition
jump
metho
the
d.
The
The


tegrals
of
v
dieren
e
t
functions
ph
on
ys-
t
ical
hing
mo
and,
dels,
a
discretization
the
sc
go
hemes
through
or
in
non-matc
of
hing
supp
triangulations
In
at
w
the
dimen-
in
for
terfaces
elemen

discretizations,
b
supp
e


e
tly
in
realized
at
in

terms
terfaces
of
to
mortar

elemen
the
t
are
metho
dimensional
ds.
es.
T
three
o
the
preserv
tersection
e
the
the
orts
global
of
optimalit
gles
y
quadrilaterals)
of
ecomes
the
hard
lo
The

of
adapted
form
discretizations,
to
the
aluate
in-
in
terfaces

b
the
et
of
w
imple-
een
tation.
the
w
dieren
er,
t
a
regions
Galerkin
ha
ximation
v
replacing
e
exact
to
tegration
b
y
e
quadrature
handled
ula
appropriately
only
.
the
Due
v
to
or
its
side
high
es
exibilit
yield
y
results
,
℄
this
n


h
t
has
o
b
The
een
est
analysed
ximation
and
requires
implemen
quadrature
ted
ula
in
on
man
sla
y
e
situations.
and
The

main
one
feature
the
of
side.
the
a2
Yv
on
Mada
8
e
H
y
[
,
the
F
(

no
Rap
in
etti,
)
Barbara
(
I.
problem
W
that
ohlm
unique
uth
,

k
w
k
e
k
are
n
led

to
0
w
1
ork
form
with
1
dieren
v
t

test
the
and
to
trial
;
spaces,

resulting
k
in
the
a
and
P
the
etro
is
v-Galerkin
N

us
h.
1
W
u
e
u
fo
is

)
our
on
atten
follo
tion
f
on
nd
the

grad
of

quadrature
1
form
Problem
ulas
the
on
e
the
is
discretization
sub
error.
;

h
results

are
`
pre-
=
sen
:
ted
sk
for
=
a
k
scalar
e
elliptic
p
equation
the
in
h
t
(
w
0
o

and
.
three
tro
dimensions.
space
Bew
;D
are
=
of
H
the
)
Æ
D
y
,
of

the
(
pro
all
of
v
in
.
three
admits
dimensions,
v
w
Giv
e
L
sho
)
w
2
only
;D
n

umerically

that

the
=
optimalit
v
y
2
of
;D
the
:
metho
admits
d
due
is
theorem
preserv
oincar
ed
Supp
if
that
suitable
osed
quadrature
v
form

ulas
=
are
:
used.
,
W

e
k
are
;
a
\
w
;
are

that
`
few
=
examples
W
are

not
S
Æ
osition
t
k
to
:

;
ab
will
out
to
this
union
sub
subsets

oundaries
but
domains
they
often,
giv

e
meas
a
D
rst
>
insigh
and
t
=
on

it.
D
Asymptotically
Let
,
in
the

error
functional
in
H
the
0
energy
(
norm
)
is
f
O
2
(
1
h

)
j
and
j
O
=
(
g
h
that
2
the
)
in
for
1
the

L
of
2
C
-norm.

2
anishing
Problem
D
form
Then,
ulation
(1)
and
the
notation
wing
The
ariational
k
ulation:
ey
en
example
2
for
2
the

denition
;
and
u
the
H
analysis
0
of
(
the
)
mortar
h
elemen
R
t
a
metho
u
d
grad
is
d

the
R
follo
f
wing
d

v
order
H
elliptic
0
b
(
oundary
)
v
(2)
alue
(2)
problem
a
div
solution,
(
to
a
Lax-Milgram
grad
and
u
P
)

=
lemma.
f
ose
in
w

;

u
in
=
non-o
0
erlapping
on
domains
D
k
;
k

1
u
:

:
n
K
=

0
that
on
=
N
K
;
=1
(1)
k
where
(
a
k
is

a
=
Æ

tly
\
smo
`
oth
k
matrix
8
that
6
is
`
uniformly
(3)
p
e
ositiv
tro
e
the
denite
eleton
in
of
the

b
S
ounded
[
op
;`
en
`
set
(4)

`

N
I
D
R
b
3
assumed
,
b
D
the
is
of
a
olygonal
subset
of
of
b
the
of
b
sub
oundary

:

of
a

osition
with