The mortar edge element method on non matching grids for eddy current calculations in moving structures

Description

The mortar edge element method on non-matching grids for eddy current calculations in moving structures F. Rapetti y Abstract The subject presented in this paper concerns the approximation of the eddy current problem in non-stationary geometries with sliding interfaces. The physical system is supposed to be composed of two solid parts: a xed one (stator) and a moving one (rotor) which slides in contact with the stator. We consider a two dimensional mathematical model based on the transverse electric formulation of the eddy currents problem in the time domain and the primary unknown is the electric eld vector. The rst order approximation of the problem that we propose here is based on the mortar element method combined with the edge element discretization in space and an implicit Euler scheme in time. Numerical results illustrate the accuracy of the method and allow to understand the in uence of the rotor movement on the currents distribution. Key words. domain decomposition method, eddy current problem, edge element approximation, electric eld as primary variable, moving structures, non-matching grids 1 Introduction The computation of the space and time distribution of induced currents in electromagnetic systems is of great importance for performance predictions and devices design. One important aspect to take into account is the presence of moving structures. The subject of our research activity is the development of simulation tools to eectively predict the induced current distribution in non-stationary geometries with sliding interfaces.

• sliding-mesh mortar


• dimensional congurations


• rotor domain


• rotor


• electric eld


• mesh elements


• variable has


• go into

Informations

• Publié par: pefav
• Nombre de lectures: 15
• Langue: English

Extrait

e
The
ld
mortar
our
edge
eddy
elemen
implicit
t
con
metho
purp
d
mo
on
the
nonatc
CNRS
hing
osed
grids
v
for
the
eddy
elemen
curren
ctric
t
time
calculations
are
in
on
mo
ersit
ving
1
structures
eddy
F
otor
Rapetti
ers

e
y
t
Abstract
analysed
The
t
sub
a
ject
the
presen
aim
ted
distribution
in
oid
this
ork
pap
that
er
prop
concerns
t
the
-
appro
cedex
ximation
de
of
e
the
the
eddy
in
curren
o
t
to
problem
other
in
[2
nontationary
ysical
geometries
ed
with
slidingesh
sliding
as
in
t
terfaces
precision
The
article
ph
d
ysical
Still
system
dimensional
is
e
supp
phenomena
osed
.
to
ect
b
on
e
w
comp
er
osed
a
of
equations
t
ariables
w
their
o
with
solid
the
parts
st
a
of
ed
space
one
in
tator
P
and
506,
a
LGEP
mo
ersit
ving
ette
one
particular
otor
fo
whic
tion
h
of
slides
t
in
domain
con
t
tact
parts
with
w
the
in
stator
with
W
In
e
]
consider
the
a
ed
t
ha
w
een
o
the
dimensional
ector
mathematical
ite
mo
d
del
ted
based
de
on
ork
the
terms
transv
robustness
erse
of
electric
to
form
same
ulation
the
of
framew
the
orking
eddy
w
curren
del
ts
ts
problem
describ
in
ph
the
the
time
ve
domain
w
and
studying
the
the
primary
emen
unkno
curren
wn
e
is
in
the
Moreo
ele
a
ctric
presence
ld
v
ve
in
ctor
e
.
Lagrangian
The
problem
st
ed
order
wn
appro
one
ximation
stator
of
rotating
the
W
problem
then
that
discretization
w
com
e
edge
prop
d
ose
of
here
sc
is

based
9029
on
Sud
the
,
mortar
Orsa
element
rance
metho
UMR
d
elec
com
,
bined
91192
with
F
the
In
edge
w
elemen
will
t
cus
discretization
atten
in
on
space
description
and
the
an
curren
implicit
x
Euler
a
sc
comp
heme
of
in
w
time
solid
Numerical
one
results
allo
illustrate
ed
the
rotate
accuracy
electric
of
tact
the
the
metho
tator
d
pap
and
[1
allo
and
w
]
to
in
understand
olv
the
ph
inence
phenomena
of
v
the
b
rotor
describ
mo
through
v
magnetic
emen
v
t
A
on
mortar
the
elemen
curren
metho
ts
w
distribution
presen
Key
in
w
no
ords
elemen
domain
framew
decomp
and
osition
in
metho
of
d
and
eddy
The
curren
ose
t
this
problem
is
edge
presen
elemen
the
t
metho
appro
in
ximation
edge
electric
t
ld
ork
as
w
primary
with
v
t
ariable
o
mo
mo
ving
of
structures
curren
nonatc
w
hing
will
grids
e
1
related
In
ysical
tro
through
duction
ele
The
ld
computation
ctor
of
As
the
e
space
at
and
the
time
of
distribution
rotor
of
v
induced
t
curren
the
ts
t
in
w
electromagnetic
will
systems
ork
is
the
of
domain
great
v
imp
to
ortance
v
for
the
p
of
erformance
con
predictions
ectiv
and
term
devices
the
design
w
One
w
imp
in
ortan
v
t
the
asp
equations
ect
solv
to
in
tak
o
e
frames
in
are
to
ed
accoun
the
t
and
is
other
the
with
presence
rotor
of
e
mo
ose
ving
a
structures
order
The
based
sub
the
ject
bination
of
the
our
elemen
researc
metho
h
in
activit
and
y
an
is
Euler
the
heme
dev
time
elopmen
ASCI
t
UPR
of
CNRS
sim
aris
ulation
Univ
to
y
ols
Bat
to
91403
ectiv
y
ely
F
predict
y
the
-
induced
8507
curren
Sup
t
Univ
distribution
y
in
Plateau
nontationary
Moulon
geometries
Gifurv
with
cedex
sliding
rance
in
terfacesd
2
spatial
W
t
e
del
remark
the
that
et
in
sho
a
consider
Lagrangian
in
approac
is
h
h
the
w
mesh
ortan
no
as
des
N
alw
ulation
a
domain
ys
idea
coincide
e
with
o
the
functions
same
of
material
angle
particles
write
throughout
e
the
n
mo
The
v
pap
emen
in
t
w
Due
the
to
on
this
metho
feature
ws
ev
meshes
en
w
if
Lagrangian
the
[7
stator
m
and
In
the
ultipliers
rotor
terface
meshes
w
ha
without
v
h
e
Concerning
b
system
een
problem
generated
mortar
in
3
suc
metho
h
the
a
v
w
the
a
estimates
y
[12]
that
emen
the
electromagnetic
t
ite
w
tac
o
space
sets
dimensional
of
curren
no
mortar
des
in
lying
nononforming
on
hnique
the
t
sliding
at
in
and
terface
metho
coincide
the
at
b
the
turn
initial
those
conuration
k
it
of
could
order
not
on
b
of
e
space
the
of
case
at
when
y
the
allo
rotor
whole
part
and
has
vy
mo
mesh
v
terface
ed
eac
If
organization
the
mo
appro
equations
ximation
form
metho
w
d
in
do
metho
es
time
not
e
allo
asp
w
Section
for
ulations
nonatc
d
hing
of
grids
the
at
w
the
problem
in
ergence
terface
discussed
either
is
a
o
o
of
c
that
ktep
mortar
tec
three
hnique
elec
is
b
adopted
tly
i
conurations
the
nononforming
time
ximation
step
w
is
ld
deed
the
dep
problem
ending
so
on
t
the
stly
rotation
].
sp
is
eed
v
and
osition
elemen
h
t
indep
size
th
at
nonatc
this
in
slip
the
in
domains
terface
the
as
is
in
imp
[3],
conditions
or
in
the
means
distortion
ultipliers
of
to
the
t
elemen
osed
ts
8].
has
argumen
to
explicit
b
particular
e
space
someho
ensure
w
prop
faced
discrete
Classical
con
remedies
edge
in
d
v
Lagrangian
olv
the
e
shap
a
tial
partial
sliding
or
p
total
using
reeshing
metho
of
us
the
with
computational
comp
domain
ed
this
rotating
pro
osing
cess
ts
and
een
most
t
of
the
all
the
the
ciated
in
time
terp
pap
olation
Section
of
deduce
the
from
unkno
Maxw
wns
w
at
v
the
of
new
Section
mesh
discuss
elemen
b
ts
ducing
increase
elemen
the
and
computational
t
cost
In