Biomedical Engineering Reference
In-Depth Information
p 1 (
r
,
z | ∂Σ a )
n
+
2
με r , z (
u 1 (
r
,
z | ∂Σ a ))
n
=
n d
1
2
p 2 (
r
, φ ,
z | Γ a )
n
+
2
με (
u 2 (
r
, φ ,
z | Γ a ))
φ
for r
∂Σ a .
(9.47b)
π
0
In this case the opposite situation occurs, as expressed by Eqs. (9.37) and (9.39),
respectively.
9.3.2 Coupling 3D and 1D models
In this case let us assume that over
Ω 1 a 1D model is set up, so we have that
Ω 1
I
( I an interval), and we keep a full 3D model in
Ω 2 (see Fig. 9.3). Specifically, for all
t we consider the following kinematics
u 1
u 1 , z (
z
) ∈Q 1
u 2
u 2 (
x
) ∈Q 2 ,
(9.48)
where u 1 , z is the component of u 1 in the axial direction e z , so it would be equivalent
to define u 1 as u 1 (
z
)=
u 1 , z (
z
)
e z . In this case, the spaces
Q i and
T Γ a ( Q i )
, i
=
1
,
2
are
H 1
Q
=
(
)
T Γ a ( Q
)=R ,
I
(9.49)
1
1
H 1
H 1 / 2
Q 2 =
( Ω 2 )
T Γ a ( Q 2 )=
( Γ a ) .
(9.50)
Here the coupling interface corresponds to
Γ a when viewed from
Ω
2 , and to the
point
ζ a when viewed from I , as seen in Fig. 9.3. As in the previous case, we write
x | Γ a (
,
,
z | Γ a )
γ =
x
y
. For this case, noting that n coincides with e z over
Γ
a ,
1 yields
1
| Γ a |
u 1 , z (
z | ζ a )=
u 2 (
x
,
y
,
z | Γ a ) ·
n d
Γ a
at
ζ a ,
(9.51a)
Γ a
μ
u 1 , z
p 1 (
z | ζ a )+
2
z (
z | ζ a )=
p 2 (
x
,
y
,
z | Γ a )+
2
με (
u 2 (
x
,
y
,
z | Γ a ))
n
·
n
for
(
x
,
y
) Γ a ,
(9.51b)
0
=
2
με (
u 2 (
x
,
y
,
z | Γ a ))
n
·
t
for
(
x
,
y
) Γ a ,
(9.51c)
where t is any vector tangent to the coupling interface
Γ a . In Eq. (9.51a) the de-
pendence on
is eliminated, establishing the continuity in a mean sense from
the point of view of the entire area
(
x
,
y
)
Γ a , whereas in Eq. (9.51b) an extension is made
such that the value of the quantity from the 1D model, a real value, is mapped to all
the points
(
x
,
y
)
over
Γ a from the 3D model. Finally, Eq. (9.51c) establishes that the
traction over
Γ a , when viewed from the 3D model, has no component tangent to the
coupling interface. Here the continuity of the velocity field is given in the sense of
R
(mean value) while the dual variable is continuous in a strong sense, that is the
sense of H 1 / 2
.
In the same manner, for
( Γ
)
a
γ =
0 we obtain the following coupling equations
u 1 (
z
| ζ a )
n
=
u 2 (
x
,
y
,
z | Γ a )
for
(
x
,
y
) Γ a ,
(9.52a)
Search WWH ::




Custom Search