15
pages

- approximation
- dimensional subspace
- then eliminated
- dependent problems
- can then
- mortar approximations
- sciences program
- applied mathematical
- feti

Voir plus
Voir moins

Vous aimerez aussi

Sud

A

is

FETI

t

PRECONDITIONER

u

F

of

OR

),

TW

curl

O

ectors

DIMENSIONAL

Comm

EDGE

Applied

ELEMENT

of

APPR

2

O

v

XIMA

curl

TIONS

space

OF

MAXWELL

h

EQUA

elations

TIONS

URL

ON

1

NONA

b

TCHING

of

GRIDS

:=

FRANCESCA

RAPETTI

pro

v

AND

2

ANDREA

u

TOSELLI

b

y

[17

Abstract

anishing

A

;

class

Orsa

of

as

FETI

and

metho

the

ds

ork

for

as

the

t

mortar

has

appro

tangen

ximation

w

of

the

a

space

v

H

ector

(

ld

2

problem

H

in

the

t

norm

w

Z

o

Z

dimensions

;

is

u

in

tial

tro

of

duced

;

and

b

analyzed

(

Edge

The

elemen

;

t

onen

discretizations

y

of

9029

lo

,

w

Eail

est

This

degree

y

are

Con

con

Direction

sidered

CNRS

The

y

metho

Mercer

d

E

prop

This

osed

part

can

of

b

Con

e

domain

emplo

diameter

y

the

ed

to

with

.

geometrically

form

conforming

(1)

and

tro

nononforming

Hilb

partitions

url

Our

b

n

;

umerical

v

results

2

sho

curl

w

2

that

The

its

;

condition

ed

n

wing

um

and

b

u

er

curl

increases

u

only

x

with

curl

the

d

n

u

um

:=

b

u

er

The

of

onen

unkno

t

wns

v

in

H

eac

on

h

@

sub

to

domains

1

and

);

is

8

indep

of

enden

H

t

with

of

tial

the

on

n

denoted

um

0

b

er

P

of

ersit

sub

506,

domains

Cedex

and

.

the

rancescaap

size

ork

of

orted

the

Europ

problem

y

Key

TMR

w

y

ords

R

Edge

ternationnelles

elemen

rance

ts

PICS

Maxw

t

ell

Sciences

equations

New

domain

10012,

decomp

toselliimsyudu

osition

ttpwwath

FETI

ork

precondition

orted

ers

y

nonatc

Sciences

hing

U

grids

Energy

AMS

DEF

sub

The

ject

classiations

unit

F

and

N

is

N

unit

N

t

1.

its

In

oundary

tro

The

duction

eak

In

ulation

this

problem

pap

requires

er

in

w

duction

e

the

consider

ert

the

H

b

;

oundary

deed

v

y

alue

url

problem

)

L

u

2

:=

L

curl

(

))

(

j

a

v

curl

L

u

(

)

)

:

+

space

A

url

u

)

=

equipp

f

with

in

follo

inner

;

duct

u

graph

(

t

;

=

)

0

:=

on

@

d

;

+

(1)

with

u

v

a

x

b

k

ounded

k

p

curl

olygonal

(

domain

;

in

)

R

:

2

tangen

.

comp

Here

t

curl

v

,

:=

a

2

ector

6

2

6

url

4

)

@

the

v

oundary

@

,

x

elongs

2

the

@

H

v

2

@

@

x

see

1

,

3

].

7

subspace

7

v

5

in

;

url

curl

)

u

v

:=

tangen

@

comp

u

t

2

@

@

is

x

b

1

H

@

url

u

).

1

ASCIPR

@

CNRS

x

aris

2

Univ

;

y

see

Building

e

91403

[17

y

].

FRANCE

The

rapettiscir

co

URL

eien

ttpwwscir

t

etti

matrix

w

A

w

is

supp

a

b

symmetric

the

uniformly

ean

p

unit

ositiv

under

e

tract

deite

ERBGT

matrix

b

alued

the

function

des

with

en

In

tries

du

A

F

ij

under

2

fund

L

478.

1

Couran

(

),

Institute

1

Mathematical

251

i

Street

j

Y

N

2,

USA

and

mail

a

.

2

h

L

yudutoselli

1

w

(

)

w

is

supp

a

in

p

b

ositiv

the

e

Mathematical

function

Program

b

the

ounded

Departmen

a

of

w

under

a

tract

y

GOR

from

zerofor

2

refer

F

Lagrange

RAPETTI

the

and

N

A

,

TOSELLI

the

F

of

or

in

an

Analogously

y

cal

D

ha

the

,

equations

w

static

e

metho

dee

st

the

non

bilinear

obtained

form

of

a

In

D

preconditioners

(

y

u

20

;

In

v

endency

)

e

:=

mo

Z

from

D

t

(

pap

a

ing

curl

of

u

conforming

curl

original

v

i

+

A

a

u

sides

v

v

b

)

ariables

d

the

x

and

;

,

u

of

;

e

v

7,

2

detailed

H

w

url

(3),

;

v

)

een

:

it

(2)

d

The

the

v

of

ariational

t

form

time

ulation

the

of

The

equation

build

(1)

Elemen

is

ETI

Find

elemen

u

metho

2

for

H

of

0

this

url

decomp

;

domains

)

:

suc

h

h

cal

that

elemen

a

of

(

con

u

onding

;

enforced

v

the

)

domain

=

algorithm

Z

b

and

f

is

b

v

solution

d

,

x

,

;

ds

v

from

2

een

H

[22

0

10

url

,

;

therein

)

5

:

pap

(3)

only

W

del

e

the

discretize

the

this

or

problem

has

using

and

edge

refer

elemen

Maxw

ts

go

also

for

kno

problem

wn

step

as

t

N

systems

edge

ed

ximation

and

elec

frequencyep

elemen

arising

ts

appro

see

ell

[24

of

].

is

These

iterativ

are

of

v

T

ector

In

alued

yp

ite

mortar

elemen

appro

ts

(1).

that

w

only

tro

ensure

solution

the

pro

con

y

tin

].

uit

h

y

of

in

the

erlapping

tangen

i

tial

1

comp

:

onen

On

t

do

across

a

the

matrix

common

the

side

discretization

of

i

adjacen

t

a

mesh

t

triangles

built

as

uit

is

solution

ph

the

ys

is

ically

y

required

ultipliers

for

terface

the

the

electric

b

and

original

magnetic

primal

lds

then

whic

solving

h

prob

are

equation

solutions

m

of

Sev

Maxw

v

ell

prop

equations

for

In

e

this

,

pap

,

er

,

w

].

e

e

consider

the

a

systems

mortar

appro

appro

e

ximation

osed

of

in

this

1

problem

2

The

11

com

,

putational

]

domain

the

is

a

partitioned

[6,

in

].

to

this

a

er

family

e

of

consider

non

mo

v

problem

erlapping

where

sub

dep

domains

on

and

time

indep

ariable

enden

on

t

frequency

triangulations

b

are

eliminated

in

w

tro

generically

duced

to

in

as

eac

ell

h

A

sub

o

domain

preconditioner

The

this

w

del

eak

is

con

st

tin

for

uit

eien

y

solution

of

linear

the

arising

tangen

the

tial

elemen

comp

appro

onen

of

t

problems

of

of

the

or

solution

enden

is

problems

then

from

enforced

quasitatic

b

ximation

y

Maxw

using

equations

suitable

aim

in

this

tegral

er

conditions

to

that

an

require

e

that

d

the

Finite

jumps

t

across

ear

the

and

sub

terconnecting

domain

t

inner

e

b

a

oundaries

edge

are

t

p

ximation

erp

problem

endicular

FETI

to

ds

suitable

ere

ite

in

elemen

duced

t

the

spaces

of

deed

ap

on

ximations

the

elasticit

edges

problems

of

[15

the

In

partition

approac

W

the

e

domain

note

is

that

osed

the

to

mortar

v

metho

sub

d

w

,

as

=

originally

;

in

:

tro

;

duced

.

in

eac

[9

sub

]

main

for

i

ite

lo

elemen

stiess

t

is

appro

from

ximations

ite

in

t

H

of

1

.

(

Mortar

;

appro

).

ximations

,

for

set

edge

righ

elemen

hand

t

is

appro

The

ximations

tin

in

y

t

the

w

corresp

o

to

and

primal

three

ariables

dimensions

then

ha

b

v

using

e

m

b

across

een

in

studied

deed

in

y

[3

sub

]

inner

and

oundaries

[4],

the

resp

FETI

ectiv

the

ely

v

.

are

There

eliminated

has

y

also

lo

b

Neumann

een

lems

additional

an

recen

in

t

Lagrange

w

ultipliers

ork

obtained

for

eral

the

ha

case

e

of

een

sliding

osed

meshes

studied

for

its

the

see

study

[14

of

16

elec

23

tromagnetic

13

lds

29

in

26

electrical

19

engines

32

see

Man

e

iterativ

[25

metho

].

for

The

solution

applications

linear

that

arising

w

mortar

e

ximations

ha

v

v

b

e

prop

in

W

mind

cite

are

particular

mainly

,

problems

,

arising

,

from

,

static

,

and

,

quasitatic

12

Maxw

21

ell

33

equations

and

ddy

to

curren

references

t

for

problems

more

see

discussion

el

FETI

preconditioner

i

f

some

or

um

Maxwell

edges

's

a

equa

domain

tions

2

3

also

T

y

o

:

our

condition

kno

pressions

wledge

2.

the

;

application

eac

of

conforming

FETI

preconditioners

to

dee

mortar

terface

appro

@

xi

in

mations

In

w

t

as

subspace

st

are

explored

nononforming

in

consider

[18]

;

and

i

then

l

tested

connected

more

do

systematically

i

in

b

[28

:

].

b

The

wing

idea

ery

is

that

fairly

g

simple

the

and

on

relies

[

on

section

the

partition

observ

ite

ation

e

that

4,

mortar

d

appro

to

ximations

cal

with

jection

Lagrange

in

m

for

ultipliers

ted

and

spaces

FETI

erlapping

form

=

ulations

;

where

[

the

p

1

oin

;

t

wise

en

con

remark

tin

to

uit

the

y

i

across

of

the

H

substructures

f

is

of

enforced

.

b

tangen

y

hosen

using

t

Lagrange

left

m

i

ultipliers

its

giv

on

e

j

rise

then

to

the

indeite

the

linear

that

systems

N

that

@

ha

j

v

ws

e

w

the

duce

same

the

form

lo

FETI

t

preconditioners

3,

can

the

then

in

b

e

e

FETI

devised

terms

for

jection

mortar

lo

appro

a

ximations

The

in

the

a

the

straigh

giv

tforw

5

ard

umerical

w

conforming

a

are

y

section

see

elemen

[28

e

].

non

In

of

this

F

pap

er

=

w

e

j

apply

=1

the

FETI

\

preconditioner

;

in

k

tro

N

duced

h

in

sub

[19

is

]

olygonal

for

in

the

W

case

F

of

not

nonedundan

e

t

e

Lagrange

of

m

y

ultipliers

the

to

the

the

sub

mortar

H

appro

max

ximation

originally

i

studied

elemen

in

H

[3].

substructur

Our

t

w

the

ork

to

gener

,

alizes

that

that

direction

for

,

FETI

on

preconditioners

or

for

domain

t

w

w

set

oimensional

en

conforming

not

edgelemen

t

i

appro

I

ximations

W

in

the

[32

also

].

eleton

As

osition

opp

of

osed

F

to

not

the

:

H

i

1

i

case

=

the

=1

generalization

i

of

follo

FETI

In

pre

2,

conditioners

e

to

tro

mortar

a

appro

of

ximations

domain

requires

and

some

cal

mo

elemen

diations

spaces

in

section

H

w

url

consider

).

mortar

More

and

precisely

section

,

w

the

presen

coarse

our

comp

metho

onen

in

ts

of

of

pro

the

on

preconditioners

a

need

wimensional

to

and

b

lo

e

preconditioner

mo

ex

did

for

here

pro

in

and

order

preconditioner

to

then

obtain

en

a

section

scalable

and

metho

n

d

results

and

geometrically

a

and

suitable

partitions

scaling

presen

matrix

in

Q

6.

has

Finite

to

t

b

W

e

st

in

a

tro

v

duced

partition

see

the

section

,

5.

H

As

(

sho

i

wn

i

in

1

[28

],

no

N

mo

N

diation

i

app

ears

=

to

;

b

k

e

necessary

=

for

;

no

dal

<

ite

elemen

)

ts

suc

in

that

H

h

1

domain

.

i

Finally

a

,

p

w

op

e

set

note

R

that

.

in

e

this

that

pap

H

er

es

w

need

e

b

only

geometrically

consider

W

problems

denote

without

diameter

jumps

of

b

the

H

co

and

eien

maxim

ts

of

F

diameters

or

the

conforming

domains

appro

y

ximations

:

FETI

:=

metho

1

ds

i

that

N

are

H

robust

g

with

The

resp

ts

ect

F

to

are

large

called

v

es

ariations

Let

of

i

jumps

e

of

unit

the

t

co

@

eien

i

ts

c

ha

so

v

follo

e

the

b

of

een

i

dev

elop

is

ed

the

and

F

studied

ev

see

sub

[26

,

,

19

e

],

the

but

of

the

op

case

edges

of

do

no

lie

dal

@

or

b

edge

f

elemen

j

t

2

appro

i

ximations

.

on

e

nonatc

dee

hing

in

grids

,

still

called

needs

k

to

of

b

decomp

e

as

explored

union

and

the

is

of

left

H

for

do

a

lie

forthcoming

@

pap

:=

er

[

The

=1

outline

of

n

the

remainder

N

of

i

this

[

pap

2I

er

i

is

asctors

4

space

F

ector

RAPETTI

i

and

subscript

A

of

TOSELLI

ectors

W

e

e

v

also

u

dee

b

the

e

lo

i

cal

spaces

:=

of

in

restrictions

the

of

denote

v

freedom

ectors

e

in

b

H

the

0

url

e

;

Y

)

)

to

=

@

i

pro

:

W

H

the

?

consist

url

Throughout

;

e

v

i

i

)

the

:=

with

f

degree

u

(

i

W

2

corresp

H

W

url

i

;

k

j

i

i

)

duce

j

=

u

i

i

url

of

t

@

=

f

0

t

on

@

@

g

h

\

=1

@

w

diren

i

that

g

the

:

t

F

i

or

w

ev

v

ery

same

substructure

i

i

i

ector

,

old

w

u

e

of

consider

using

a

(

triangulation

.

T

onding

i

in

,

.

made

duct

of

denoted

triangles

u

or

of

rectangles

letter

Let

the

E

X

i

onding

b

i

e

u

the

e

set

k

of

2

edges

(4)

of

in

T

pro

i

h

.

:=

F

=1

or

N

ev

H

ery

edge

the

e

tial

2

h

E

i

i

i

,

u

w

i

e

restricted

x

i

a

j

direction

X

giv

and

en

space

b

=

y

Y

a

i

unit

note

v

ha

ector

hosen

t

deition

e

spaces

.

y

The

].

length

W

of

piecewise

the

tial

edge

@

e

@

is

pap

denoted

will

b

wing

y

tions

j

use

e

for

j

in

.

tangen

The

in

lo

W

cal

generic

triangulations

in

are

a

assumed

with

shap

,

eegular

,

and

v

quasiniform

degrees

and

in

they

same

do

sup

not

),

need

i

to

k

matc

freedom

h

the

across

,

the

is

inner

)

b

generic

oundaries

the

of

X

the

is

sub

y

domains

letter

W

and

e

v

dee

of

h

the

as

u

the

will

maxim

notation

um

of

of

and

the

the

meshizes

of

of

(

the

)

triangulations

:=

W

i

e

t

next

k

consider

e

the

;

lo

k

w

E

estrder

:

N

W

next

ed

tro

the

elec

duct

ite

X

elemen

(

)

t

X

E

N

spaces

i

originally

X

in

tro

Y

duced

=1

in

?

[24

;

],

i

deed

;

on

spaces

eac

tangen

h

v

sub

W

domain

(

i

)

as

W

X

:=

h

(

(

i

i

t

)

)

=

i

X

to

i

:=

n

f

u

u

i

2

2

i

H

;

?

the

url

duct

;

W

()

i

W

)

N

j

i

u

W

i

:

j

e

t

that

2

e

R

v

(

c

t

a

)

t

;

of

t

trace

2

than

T

emplo

i

ed

g

[32

;

Here

where

spaces

in

i

the

of

case

constan

of

tangen

triangular

ve

meshes

on

w

e

n

ha

.

v

this

e

er

R

e

(

use

t

follo

)

con

:=

en

W

will

1

the

+

notation

the

3

ectors

x

X

2

and

tial

2

ectors

W

3

.

x

e

1

a

v

j

function

X

k

using

2

b

R

letter

the

:

i

W

e

e

i

recall

and

that

column

the

ector

tangen

its

tial

of

comp

deed

onen

(4),

t

the

of

letter

a

the

v

erscript

ector

i

u

e

i

(

2

)

X

Its

i

h

is

of

constan

corresp

t

to

on

edge

the

k

edges

deed

of

(4),

the

u

triangulation

i

T

k

i

A

,

v

and

in

that

pro

the

space

degrees

r

of

)

freedom

also

can

b

b

a

e

old

c

e

hosen

,

as

the

the

onding

v

ector

alues

degrees

of

freedom

the

y

tangen

same

tial

e

comp

.

onen

e

t

use

on

same

the

for

edges

spaces

functions

e

i

k

W

(

and

u

corresp

i

spaces

)

degrees

=

freedom

u-1 1

-11Ω Ω

3-11 -1

1

-1 1

1 -1

Ω -1-1

1 11 Ω

4 1 -1 -1

1 -1

1-1 -1 1 1

-1-1

1 11-1

1

-1 Ω 2

1

-1

of

5

m

's

eleton

Maxwell

(

or

section

f

f

preconditioner

k

tions

u

equa

t

i

The

t

studied

:

as

of

W

space

I

mortar

u

;j

FETI

in

call

;j

partition

p

2

recall

A

energy

this

i

b

all

(

tangen

k

i

o

3.

mortars

presen

I

dev

e

e

m

ose

6

union

M

,

I

e

Fig

general

1

example

.

indices

Example

w

of

:=

p

)

artition

mortar

of

w

the

m

domain

W

i

.

;

We

n

show

an

the

e

dir

W

e

;

ctions

en

of

(

the

k

sub

2

domain

that

b

are

oundaries

W

given

tro

by

to

the

m

unit

W

ve

H

ctors

minimizes

f

t

i

i

H

g

)

,

v

those

i

of

comp

the

to

e

@

e

@

dges

mortar

on

metho

the

in

interfac

as

e

ed

,

[3].

and

the

the

c

c

splitting

orr

disjoin

esp

some

onding

;j

values

w

of

.

the

that

de

not

gr

see

e

for

es

decomp

of

set

fr

onds

e

hoice

e

denote

dom

I

t

m

(

;

i

h

)

is

.

:

Giv

the

en

denote

the

y

unit

m

v

g

ectors

ha

t

M

i

m

,

n

the

if

column

n

v

2

ectors

F

t

m

(

,

i

b

)

the

are

;j

deed

(

b

)

y

)

t

y

(

:=

i

k

)

)

k

to

:=

u

t

k

i

e

v

t

k

e

the

k

ectors

;

to

e

Before

k

the

e

@

a

t

i

e

n

:

@

e

that

;

i

e

i

k

the

2

a

E

i

i

H

:

u

W

;

e

i

will

i

need

among

these

the

tangen

ectors

tial

X

v

with

ectors

tial

in

onen

the

equal

deition

u

of

on

our

FETI

n

metho

.

d

A

see

condition

section

mortar

5.

d

W

ted

e

this

remark

w

that

originally

in

elop

case

and

that

in

all

W

the

consider

edges

sk

e

and

k

ho

on

a

@

of

the

i

t

ha

of

v

edges

e

k

the

g

same

that

direction

e

of

mortars

the

W

b

note

oundary

this

@

is

in

i

unique

,

Figure

the

eft

en

an

tries

of

of

osition

the

unique

v

of

ector

corresp

t

to

(

c

i

and

)

e

are

it

equal

y

to

M

one

f

Figure

=

1

k

sho

j

ws

suc

an

that

example

;j

of

a

a

g

partition

T

with

simplify

the

notation

directions

e

of

the

the

b

sub

f

domain

j

b

2

oundaries

M

and

.

of

e

the

v

e

:=

edges

[

on

=1

the

;

in

\

terface

=

,

;

and

m

the

=

corresp

and

onding

m

degrees

I

of

.

freedom

or

t

y

(

in

i

M

)

w

.

denote

Finally

y

,

m

for

space

i

k

=

(

1

=

;

k

j

2

M

;

giv

N

b

,

W

w

;j

e

f

dee

u

the

discrete

k

harmonic

t

extensions

restricted

with

k

resp

j

ect

k

to

X

the

g

bilinear

W

forms

note

a

the

ectors

i

W

(

;j

also

;

restrictions

v

)

in

in

k

to

k

the

.

in

in

terior

ducing

of

mortar

w

i

need

H

x

i

last

:

oin

W

Let

i

b

!

a

X

i1

Ω 0.911,2

0.8Γ

0.7

1,1

0.6Γ

0.5Ω3

0.4

0.3

2,1

0.2Γ

0.1Ω2

0

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

dge

Here

shap

mortar

(

imp

u

ds

2

form

the

ecomes

these

deed

As

decomp

A

b

mortar

the

elements

[3

c

elements

TOSELLI

sp

RAPETTI

ac

essary

subsp

[9

b

of

the

dges

and

m

j

remark

m

F

tial

at

adjacen

6

m

required

Fig

2

;

.

transmission

A

n

k

example

stringen

wher

lds

e

and

the

in

domain

w

is

a

de

suitable

c

ultiplier

omp

onen

ose

the

d

The

into

for

thr

for

e

dal

e

L

(r

chosen

e

ac

ctangular

only

nonverlapping

it

sub

is

domains

e

The

ac

skeleton

ciate

has

no

two

e

p

is

artitions

2

gur

2

e

g

on

a

the

.

left

b

one

are

in

means

terms

u

of

constrain

mortars

t

olid

v

dark

v

line

k

and

(5)

the

e

other

wing

in

k

terms

of

m

nonortars

to

ong

condition

dashe

w

d

in

lines

diren

On

hing

the

usually

right

mortar

an

metho

example

imp

of

ts

discr

eak

etization

means

of

m

the

Lagrange

sub

consists

domains

tial

by

of

me

functions

ans

edges

of

Remark

triangular

of

grids

ac

that

e

do

simpler

not

dal

match

the

at

ase

the

e

interfac

ange

es

b

b

a

etwe

of

en

on

adjac

dges

ent

a

sub

e

domains

the

edge

ase

with

ne

m

de

=

the

(

multiplier

k

sinc

;

is

j

to

)

not

and

se

u

for

2

The

X

ortar

h

w

:

f

for

2

almost

j

ev

m

ery

I

x

W

2

this

m

of

,

ectors

there

transmission

exists

in

an

w

index

sub

l

w

(1

b

these

l

A

W

N

satisfy

),

Z

l

k

6

u

=

k

k

t

suc

0

h

W

that

=

x

j

2

M

m

set

\

can

@

in

the

l

a

.

t

A

and

t

k

this

t

p

at

oin

b

t

a

x

o

w

t

e

since

ha

t

v

o

e

are

t

general

w

on

o

t

lds

nonatc

namely

grids

u

is

k

done

and

nononforming

u

domain

l

osition

.

ds

Since

e

the

ose

domain

constrain

decomp

in

osition

w

is

form

in

y

general

of

nononforming

Lagrange

the

ultipliers

v

the

alue

m

of

space

l

of

dep

tangen

ends

comp

on

ts

x

the

and

e

w

at

e

mortar

denote

see

b

].

y

3.1.

J

deition

l

the

the

sp

set

e

of

the

indices

dge

l

is

(1

than

no

l

ite

In

N

no

)

c

suc

the

h

ac

that

of

m

agr

\

multipliers

@

annot

e

l

as

6

sp

=

e

;

tr

.

es

W

the

e

e

then

but

dee

as

u

suitable

k

ac

(

of

x

In

)

e

:=

c

u

it

l

not

(

c

x

to

)

cr

;

ase

x

dimension

2

the

m

sp

\

e

@

e

information

l

asso

;

d

l

e

2

and

J

to

l

des

:

e

The

]

function

mor

u

details

k

Lagrange

is

ultiplier

deed

space

for

no

almost

deed

all

y

x

v

2

L

m

()

.

v

In

m

general

W

it

;

is

2

not

M

the

:

tangen

e

tial

that

comp

is

onen

space

t

tangen

at

v

m

on

of

The

a

conditions

ld

the

u

terface

2

et

H

een

?

t

url

domains

;

then

eakly

i

osed

):

y

it

of

can

Lagrange

indeed

ultipliers

corresp

solution

ond

2

to

is

tangen

to

tial

the

comp

ts

onen

m

ts

u

from

diren

k

t

k

sub

t

domains

)

whic

h

k

share

=

a

;

subset

2

of

m

m

m

and

(

liv

;

e

)

on

I

diren

:

t

The

grids

of

The

conditions

equalit

b

y

expressed

b

matrix

et

in

w

follo

een

w

u

y

ku

FETI

and

preconditioner

e

f

S

or

tin

Maxwell

b

's

freedom

equa

wing

tions

2

7

appro

Let

cal

w

k

k

(

i

h

b

i

e

B

the

:=

basis

The

function

system

asso

assem

ciated

(

to

that

the

e

k

b

h

the

mesh

i

edge

W

of

prop

@

=

see

i

write

n

0

@

)

.

W

f

e

ultiplier

in

0.

tro

mortar

duce

(3).

t

cal

w

forms

o

),

matrices

The

C

the

and

substructure

D

b

b

f

y

righ

C

i

k

complemen

l

W

:=

to

Z

i

m

the

(

satisfy

w

(

k

u

i

(

u

t

31

i

w

)(

as

w

f

l

:=

i

.

2

t

(1)

i

g

)

(1)

ds

3

=

is

the

k

t

l

the

j

from

e

elemen

k

of

j

e

;

the

e

matrices

k

the

;

e

;

l

the

v

m

of

;

not

D

terface

k

to

n

can

:=

in

Z

blo

m

elimination

(

i

w

the

k

hand

i

S

b

t

h

i

matrices

)(

)

w

!

n

relativ

r

degrees

@

t

@

i

recall

)

cal

ds

complemen

e

follo

k

y

)

m

i

;

i

e

n

i

H

@

);

[27

r

F

;

32

r

can

2

mortar

J

u

m

;

u

where

where

m

6

=

.

i

(

.

7

Then

;

the

f

matc

hing

S

conditions

f

(5)

4

ha

.

v

N

e

5

the

ector

form

Lagrange

B

e

u

eak

=

y

0

u

where

of

B

linear

=

arising

C

the

D

edge

:

t

W

ximation

e

problem

remark

W

that

st

the

ble

en

lo

tries

stiess

of

relativ

C

to

and

bilinear

D

a

dep

i

end

on

the

and

particular

lo

c

load

hoice

ectors

of

degrees

degrees

freedom

of

are

freedom

on

deed

in

in

only

(4).

elong

The

one

matrix

and

B

b

can

eliminated

also

parallel

b

y

e

c

written

Gaussian

as

Let

B

(

=

)

h

e

B

resulting

(1)

t

B

sides

(2)

let

(

)

e

B

Sc

(

ur

N

t

)

S

i

i

;

:

where

i

the

W

lo

;

cal

e

matrices

the

B

of

(

on

i

)

n

act

.

on

e

v

that

ectors

lo

in

Sc

W

ur

i

ts

.

the

The

wing

en

ert

tries

u

of

i

B

t

do

(

not

)

b

(

elong

)

in

a

general

i

to

H

f

u

0

;

;

i

1

i

;

(6)

1

e

g

,

as

].

in

ollo

the

[19,

conforming

],

case

e

describ

then

ed

our

in

problem

[32

S

]

+

but

t

since

=

w

B

e

=

are

(7)

w

u

orking

2

with

4

nonatc

(1)

hing

.

grids

u

they

N

tak

3

e

5

in

W

to

S

accoun

diag

t

S

the

;

edge

in

;

tersections

(N)

at

;

the

:=

in

6

terfaces

f

4.

.

A

.

FETI

(

metho

)

d

7

In

:

this

v

section

w

a

e

m

in

relativ

tro

to

duce

w

a

con

FETI

uit

metho

constrain

d

B

for

=

the

solution