Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
dumux-repositories
dumux-course
Commits
f7c2a629
Commit
f7c2a629
authored
Mar 31, 2020
by
Ned Coltman
Browse files
[CoupledExercise] Set the Beavers Joeseph condition using setBeaversJoeseph rather than set BJS
parent
275a56c4
Changes
6
Hide whitespace changes
Inline
Side-by-side
exercises/exercise-coupling-ff-pm/README.md
View file @
f7c2a629
...
...
@@ -131,7 +131,7 @@ values.setDirichlet(Indices::velocityXIdx); // assume no slip on interface
```
with a Beavers-Joseph-Saffman (BJS) boundary condition for the respective momentum balance:
```
cpp
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
```
at the position where the coupling boundary conditions are set in
`interface/freeflowsubproblem.hh`
.
...
...
@@ -187,7 +187,7 @@ in `main.cc` in the subfolder `interface`.
The same is true for the BJS condition, however, here we need to consider the tangential direction:
```
cpp
values
.
setB
JS
(
1
-
scvf
.
directionIndex
());
values
.
setB
eaversJoseph
(
1
-
scvf
.
directionIndex
());
```
The final result should look something like this:
...
...
exercises/exercise-coupling-ff-pm/models/freeflowsubproblem.hh
View file @
f7c2a629
...
...
@@ -184,7 +184,7 @@ public:
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
);
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
+
1
);
values
.
setCouplingNeumann
(
Indices
::
momentumYBalanceIdx
);
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
}
return
values
;
...
...
exercises/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh
View file @
f7c2a629
...
...
@@ -209,7 +209,7 @@ public:
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
+
1
);
values
.
setCouplingNeumann
(
Indices
::
energyEqIdx
);
values
.
setCouplingNeumann
(
Indices
::
momentumYBalanceIdx
);
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
}
return
values
;
}
...
...
exercises/solution/exercise-coupling-ff-pm/interface/freeflowsubproblem.hh
View file @
f7c2a629
...
...
@@ -196,11 +196,11 @@ public:
values
.
setDirichlet
(
Indices
::
velocityXIdx
);
// assume no slip on interface
#elif EXNUMBER == 2
// set the Beaver-Joseph-Saffman slip condition for the tangential momentum balance equation
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
#else
// set the Beaver-Joseph-Saffman slip condition for the tangential momentum balance equation,
// consider orientation of face automatically
values
.
setB
JS
(
1
-
scvf
.
directionIndex
());
values
.
setB
eaversJoseph
(
1
-
scvf
.
directionIndex
());
#endif
}
...
...
exercises/solution/exercise-coupling-ff-pm/models/freeflowsubproblem.hh
View file @
f7c2a629
...
...
@@ -187,7 +187,7 @@ public:
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
);
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
+
1
);
values
.
setCouplingNeumann
(
Indices
::
momentumYBalanceIdx
);
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
}
return
values
;
...
...
exercises/solution/exercise-coupling-ff-pm/turbulence/freeflowsubproblem.hh
View file @
f7c2a629
...
...
@@ -220,7 +220,7 @@ public:
values
.
setCouplingNeumann
(
Indices
::
conti0EqIdx
+
1
);
values
.
setCouplingNeumann
(
Indices
::
energyEqIdx
);
values
.
setCouplingNeumann
(
Indices
::
momentumYBalanceIdx
);
values
.
setB
JS
(
Indices
::
momentumXBalanceIdx
);
values
.
setB
eaversJoseph
(
Indices
::
momentumXBalanceIdx
);
}
return
values
;
}
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment